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INTRODUCTION 


y,' 


Canopy  sun  glint  is  a major  cue  in  the  visual  detection  of  aircraft.  There- 
fore, it  would  be  desirable  to  have  a method  to  calculate  the  magnitude  of 
these  sun  glints  for  a proposed  or  experimental  canopy.  The  computer  pro- 
gram described  in  this  report  is  an  attempt  to  provide  an  analytical 
method  to  solve  this  problem  for  smooth  reflecting  surfaces.  However,  cer- 
tain limitations  were  necessary  to  maintain  a degree  of  simplicity.  These 
limitations,  which  include  the  restriction  that  emy  section  of  a reflecting 
surface  have  a quadratic  representation,  are  discussed. 
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DESCRIPTION  OF  PROBLEM 


An  aircraft  canopy  reflects  the  rays  of  the  sun  at  azimuth  and  elevation 
angles  that  are  a function  of  the  geometry  of  the  canopy  and  the  orienta- 
tion of  the  aircraft  with  respect  to  the  sun.  As  the  aircraft  rotates,  as 
shown  in  Figure  1,  the  sun  glint  also  moves;  and,  as  the  orientation  of  the 
aircraft  changes  with  respect  to  the  sun,  more  or  different  canopy  panels 
can  beccme  the  reflecting  surfaces.  At  any  one  instant,  canopy  sun  glints 
may  be  visible  at  more  than  one  set  of  azimuth  and  elevation  coordinates . 

The  objective  then  in  designing  a canopy  to  minimize  visual  detection  is  to 
limit  the  size  cind  number  of  sun  glints  occurring  within  the  band  of  co- 
ordinates which  are  probable  positions  for  hostile  observers.  To  evaluate 
the  effectiveness  of  such  a design,  there  is  a choice  between  experimental 
testing  and  cuialysis.  There  are  at  least  two  problems  with  experimentally 
testing  new  canopy  designs.  First,  if  a scaled  model  is  used,  getting  the 
proper  curvature  on  the  panels  could  be  time-consuming  and  expensive. 

Second,  if  a mockup  is  used,  greater  distances  are  required  to  simulate  a 
real  situation.  In  either  case,  the  evaluation  of  a large  family  of  curves, 
or  Vciriations  of  a basic  design,  would  be  difficult.  This  is  not  to  say 
that  an  aned.ytical  approach  is  without  problems;  there  are  certain  limita- 
tions that  must  be  imposed.  But  the  simple  fact  that  a large  number  of 
different  designs  can  be  evaluated  rapidly  makes  the  analytical  scheme 
attractive.  The  problem  then  is  to  develop  an  analytical  method  to  deter- 
mine the  coordinates  of  the  sun  glint. 
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DEFINITI(»S 


Before  proceeding  any  further  with  the  discussion,  some  useful  definitions 

are  presented: 

• Reflecting  surfaces  (panels)  - Refers  to  the  transparent  surfaces  of 
interest;  i.e.,  all  see-through  surfaces  of  the  canopy. 

• Fences  - Opaque  surfaces  which  may  either  prevent  the  rays  of  the  sun 
frcm  striking  a reflecting  surface  or  obscure  the  reflection  from  the 
observer.  The  program  assumes  that  all  fences  can  be  described  by  flat 
surfaces.  This  seems  permissible,  since  it  is  not  the  curvature  of  the 
fence  surface  which  is  of  interest  but  merely  where  in  space  the  fence 
is  located. 

• Sun  vector  - A sliding  vector  which  represents  the  rays  of  the  sun  as 
they  strike  a reflecting  surface.  The  orientation  of  this  vector  as 
referred  to  inertial  axes  (earth-fixed)  is  defined  by  the  elevation  of 
the  sun.  The  vector  has  unit  magnitude  and  terminates  at  all  boundary 
points  of  all  reflecting  surfaces.  The  vector  can  be  viewed  as  simply 
a sliding  cir.Tow  which  remains  parallel,  but  opposite  in  direction,  to 
^ul  ^u:row  whose  tail  is  at  the  origin  of  the  inertial  axes  and  points 
toweurd  the  sun,  as  shown  in  Figure  2. 

• Reflection  vectors  - Each  boundeiry  point  on  a reflecting  surface  has  a 
reflection  vector  associated  with  it,  as  shown  in  Figure  3.  Mathe- 
matically, it  is  formed  by  taking  the  sun  vector  and  retaining  the 
portion  which  is  tangent  to  the  reflecting  surface  while  reversing  the 
portion  which  points  toward  the  surface  along  the  normal.  As  the  nor- 
mal varies  from  point  to  point,  so  does  the  reflection  vector. 

• Inertiad.  auces  - Earth-fixed  axis  system  with  origin  at  the  center  of 
rotation  of  the  aircraft.  The  positive  X axis  points  in  the  direction 
of  the  sun  while  the  positive  Z axis  points  away  from  the  center  of 
earth.  The  Y eucis  forms  right  angles  with  the  X and  Z axes  euid  has 
positive  sense  to  the  right  when  facing  the  sun  (see  Figure  4) . 
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Figure  2.  Sun  Vector. 


Figure  3.  Basic  Relationship  Between  Striking  Sun  Ray  and  Reflection  Vector. 
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Figure  4.  Axis  Systems. 
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Aircraft  headings  are  measured  relative  to  this  axis  system  and  are 
positive  as  shown  in  Figure  1.  Sun  glint  position  is  also  measured 
relative  to  these  axes  as  shown  in  Figure  5. 


• Body  axes  - Body-fixed  axis  system  with  origin  at  the  center  of  rotation 
of  the  aircraft.  The  positive  X axis  points  toward  the  tail  of  the  air- 
craft while  the  positive  Z axis  points  toward  the  top  of  the  aircraft. 

The  Y axis  is  positive  out  the  right  side  of  the  aircraft  and,  along 
with  the  X and  Z axes,  forms  an  orthogonal  system.  These  axes  are  as- 
sumed parallel  to  the  station  line,  buttline,  and  waterline  reference 
system  as  shown  in  Figure  4. 

• Sun  glints  - In  this  report,  a sun  glint  refers  to  the  single  image 
produced  by  one  reflective  surface  for  one  set  of  aircraft  and  sun  po- 
sitions. If  the  reflection  vectors  are  extended,  as  shown  in  Figure  5, 
until  they  intersect  an  imaginary  cylinder  in  space,  the  image  produced 
is  a sun  glint.  The  whole  thing  can  be  visualized  as  if  the  aircraft 
were  sitting  in  the  center  of  a circular  stadium  and  the  stands  were 
filled  with  hostile  observers.  These  stands  are  strictly  vertical  and 
form  the  rim  of  the  stadium.  The  reflection  vector  from  each  boundary 
point,  when  extended,  will  strike  someone  in  the  stands.  If  each  person 
that  is  struck  stands  up,  the  outline  of  an  image  is  formed.  The  image 
formed  from  the  reflection  vectors  of  one  reflective  surface  will  be 
called  the  sun  glint  of  that  surface  at  the  specified  aircraft  and  sun 
conditions.  Since  the  radius  (distance  to  observers)  of  the  stadium  is 
a fixed  input  value,  the  location  of  any  person  in  the  stands  is  defined 
completely  by  specifying  two  angles . The  two  angles  are  measured  rela- 
tive to  the  inertial  axes.  To  locate  a person,  fir^t  rotate  the  X 
inertial  axis  about  the  Z axis  until  it  intersects  the  stands  directly 
below  the  person  (azimuth  angle) ; then  rotate  the  X axis  a)x)ut  the  Y 
axis  until  the  X axis  points  directly  at  the  person  (elevation  angle) . 

In  the  same  way  the  boundary  points  of  the  sun  glint  are  defined  by  an 
azimuth  angle  and  an  elevation  angle.  These  coordinates  (angles)  form  the 
results  of  this  program. 
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SUN  GLINT 


Figure  5.  Definition  of  Sun  Glint. 
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• Sun  glint  signature  - This  is  the  plot  of  all  the  coordinates  of  all 
the  sun  glints  from  all  the  reflective  surfaces  for  one  sun  elevation. 
A sample  plot  for  an  imaginary  canopy  is  shown  in  Figure  fo. 


• Probability  - The  probability  of  being  observed  is  defined  to  be  the 

ratio  of  two  areas . The  numerator  is  the  area  within  the  glint  signature 
plot  which  is  swept  out  by  the  reflected  canopy  sun  glints  eis  the  air- 
craft is  rotated  360  degrees  in  heading.  The  denominator  is  the  total 
area  enclosed  by  the  rectangle  formed  by  the  maximum  and  minimum  allow- 
aJsle  elevation  euid  azimuth  cuigles  of  the  observer.  As  an  example,  con- 
sider Figure  6.  The  curea  swept  out  is  approximately  (135-  (-150)  ) X 
(5-  (-5)  ) = 2,850  degrees^.  The  total  allowable  area  is  (180-  (-180)  ) 

X (10-  (-10)  ) = 7,200  degrees^.  The  probability  would  then  be 


2,850  ^ 
7,200 
that  is 


0.396.  In  the  calculation  of  the  swept  area,  overlap  is  ignored; 
if  a portion  of  the  sun  glint  signature  plot  is  covered  by  the 


sun  glints  from  more  than  one  surface,  the  area  is  only  counted  once. 


AZIMUTHAL  POSITION  OF  OBSERVER 
RELATIVE  TO  THE  SUN  - DEG 


Figure  6.  Sun  Glint  Signature. 
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METHOD  OF  SOLUTION 


I 

The  rays  of  the  sun,  treated  as  parallel  vectors,  strike  the  canopy  and 
are  reflected  by  the  canopy.  The  method  of  solution  is  basically  simple; 
the  trimmings  make  it  con^licated.  The  reflecting  surface  is  assumed  to  be 
smooth,  and  therefore,  the  reflection  vector  becomes  a function  of  the 
angle  between  the  striking  sun  rays  and  the  normal  to  the  surface,  as  was 
shewn  in  Figure  3.  This  angle  may  vary  as  the  relative  position  between 
the  aircraft  and  the  sun  varies,  for  instance,  due  to  aircraft  heading 
changes.  It  could  also  vary  from  point  to  point  on  the  surface  as  the 
normal  to  the  surface  varies,  as  in  the  case  of  a curved  panel. 


It  is  necessary  then  to  determine  at  all  the  boundary  points  of  all  the 
reflective  surfaces  the  normals  to  those  surfaces.  One  method  considered 
was  to  input  the  normals  for  all  the  reflective  surfaces. 


However,  these  normals  are  not  readily  available  in  most  cases.  Therefore, 
an  alternate  approach  was  taken.  The  location  of  the  boundary  points  in 
terms  of  aircraft  station  lines,  buttlines,  and  water lines  was  used. 

From  these  points  the  program  fits  a surface  through  the  points  in  a least- 
squares  error  sense.  This  method  does  have  a drawback.  There  are  an  un- 
limited number  of  types  of  surfaces  that  could  be  used.  To  keep  things 
relatively  simple,  and  still  handle  the  majority  of  cases,  a general  second- 
order  surface  was  chosen.  More  complex  reflective  panels  can  be  handled 
by  breciking  them  into  smaller  panels  (sections) , thereby  allowing  each 
section  to  fit  a different  second-order  equation.  For  instance,  if  a panel 
looked  like  that  shown  in  Figure  7,  breaking  it  up  into  three  panels  as 
shown  should  yield  an  acceptable  fit. 

Once  the  equation  of  the  reflective  surface  is  known,  the  normal  at  any 
point  on  the  surface  is  known.  Then,  selecting  an  aircraft  orientation 
(pitch,  roll,  and  heading)  and  sun  elevation,  the  cuigle  between  the  sun 
vector  and  any  surface  normal  can  be  calculated.  The  reflection  vector  at 
each  boundary  point  is  then  known;  and  from  the  reflection  vectors,  the 
sun  glints  can  be  formed,  as  was  shown  in  Figure  5. 
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It  has  been  assumed  so  far  that  the  sun  rays  are  free  to  strike  every 
boundeury  point  on  the  surface.  The  program,  however,  is  cap^tble  of  de- 
termining whether  an  opaque  surface  (fence) , defined  by  boundary  points, 
ccmes  between  the  sun  and  a reflective  surface  boundary  point.  The  pro- 
gram is  also  capckble  of  determining  whether  a reflection  vector,  if  ex- 
tended, will  intersect  cui  opaque  surface  before  reaching  a possible 
observer.  Finally,  the  program  is  capcd>le  of  determining  whether  the  re- 
flection is  internal  to  the  canopy. 

Basically,  this  is  the  method  used  by  the  program  to  solve  the  sun  glint 
problem.  The  progreun  is  dressed  up  a bit  by  providing  print  plots  and 
graphic  plots  of  the  resulting  sun  glint  signatures.  The  total  scheme 
is  as  follows: 

• Input  desired  aircraft  orientations 

• Input  desired  sun  elevations 

• Input  desired  observer  restraints 

• Input  the  boundary  points  for  all  fences  and  all  reflective  surfaces 

• Curve-fit  boundary  points 

• Calculate  relative  position  of  sun  vector  and  reflective  surface 

• Determine  if  any  fences  lie  between  sun  and  reflective  surface 

• Calculate  normals  to  reflective  surface  at  boundary  points 

• Calculate  reflection  vectors 

• Determine  if  reflection  vectors  when  extended  intersect  any  fences 

• Determine  if  reflection  is  internal  to  canopy 

• Calculate  observer  positions  (sun  glints) 

• Calculate  probabilities  from  sun  glint  signatures. 
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PROGRAM  FLOW  CHARTS 


Presented  in  this  section  are  the  flow  cheu:ts  for  the  executive  (main) 
routine  and  each  subroutine  which  contains  at  least  one  CALL  statement. 

The  bloc)cs  that  are  shown  within  a particular  flow  ch2u:t  refer  only  to  the 
CALL  statements  in  that  routine.  If  a subroutine  is  called  which  calls 
^mother  subroutine,  only  the  first  CALL  statement  appears.  However,  in 
the  flow  cheurt  for  the  first  subroutine  called  will  appear  the  CALL  state- 
ment to  the  second  subroutine.  In  other  words,  only  one  level  of  calls 
appear  in  any  flow  chart.  For  example,  in  the  executive  routine  there  is 
a call  to  subroutine  CFITF,  which  would  result  in  cauLls  to  subroutines 
NORM,  CURFIT,  SIGNF,  AND  EQNSOL.  This  ceui  be  seen  by  looking  at  the  flow 
charts  for  CFITF  and  CURFIT. 


17 


CALL  INPUT 

CALL  CFITS 

READ 

INPUT 

DATA 

CURVE  FIT 

REFLECTIVE 

SURFACE 

DATA 

CALL 

GLINT 

CALCl 

REFL( 

VECTC 

FENCE  SH 

ANDINTEF 

AND  SUN ( 

JLATE 

CTION 

)RS. 

hDOWS 

(FERENCE. 

jLINTS 

CALL  MINMAK 


CALCULATE 
MINIMUM  AND 
MAXIMUM 
COORDINATES 
OF  EACH 
GLINT 


Figure  8.  Executive  Flow  Chart 
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Figure  9.  Flow  Charts  for  Subroutines  (Sheet  1 of  3). 
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Figure  9.  Flow  Charts  for  Subroutines  (Sheet  2 of  3). 
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DESCRIPTION  OF  SUBROUTINES 

< 

A general  description  of  all  the  subroutines  used  in  the  program  follows:  j 

I 

• INPUT  - Reads  data  into  progreun  from  cards.  Generates  boundary  points  j 

for  symmetrical  surfaces.  Sets  certain  variables  to  default  values,  j 

when  variables  are  left  undefined.  ' j 

i 

• OUTPTl  - Prints  out  input  data  with  suitable  headings. 

• 0UTPT2  - Prints  out  results  from  curve-fitting.  Coefficients  of  the 

curve  best  approximating,  in  the  least-squares  sense,  the  boundary  points  ; 

of  all  fences  and  reflective  surfaces  are  printed  and  labeled. 

• 0UTPT3  - Tabular  printout  of  boundary  points  for  each  reflective  surface  '■ 

with  condition  flags  and  required  observer  coordinates  (azimuth  and 

elevation).  These  observer  coordinates  are  those  shown  in  Figure  5. 

The  condition  flags  refer  to  whether  the  reflection  vector  from  each 
particular  boundary  point  is  an  internal  reflection,  whether  the  reflec- 
tion vector  is  blocked  by  a fence,  or  whether  the  boundary  point  lies  in 

1 

the  shadow  of  a fence. 

• PLOTOL  - Uses  line  printer  to  plot  sun  glint  signatures.  Points  plotted 
represent  minimum  and  maximum  values  of  elevation  and  azimuth  for  each 
glint. 

• DRPLOT  - Prepares  and  sorts  tabulated  glint  data  to  be  plotted  by 
CALCOMP  plotter. 

' 

• GPLOT  - Gives  commands  to  CALCOMP  software  routines  to  draw  glint  signa- 
tures. Three  signatures,  corresponding  to  three  sun  elevations,  are 
drawn  per  page  with  a maximum  of  three  pages  (nine  elevations)  per 

computer  run.  . 
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• PLOTS,  FACTOR,  PLOT,  AXIS,  SYMBOL,  NUMBER,  and  LINE  - CALCOMP  software 
functions  which  are  supplied  by  CALCOMP  BASIC  SOFTOARE  package. 

• GRID  - Normally  a CALCOMP  software  routine;  this  subroutine  was  written 
and  included  as  part  of  the  program  since  it  is  not  pcirt  of  the  basic 
packages  supplied  to  some  CALCOMP  users.  If  system  already  has  GRID, 
modify  subroutine  name  and  call  statement  to  make  it  unique,  i.e.,  GRIB 
instead  of  GRID. 

• CFITF  - Along  with  the  subroutines  NORM,  CURFIT,  SIGNF,  and  EQNSOL, 
curve- fits  the  boundary  points  of  a fence  to  a plane  (linear)  surface 
in  space. 

• NORM  - In  an  attempt  to  minimize  numerical  errors  during  curve- fitting, 
boundary  points  are  normalized.  This  routine  calculates  a norm  which 

is  the  root-mean- square  value  of  the  set  of  coordinates  representing  the 
boundary  points. 

• CURFIT  - Sets  up  the  normal  coefficient  matrix  and  constant  terms  for 
the  least-squares  solution.  Calculates  the  root-mean- square  error  of 
the  fit  by  plugging  boundary  points  into  the  equation  of  the  fitted 
surface. 

• EQNSOL  - Solves  a system  of  linear,  homogeneous  or  inhomogeneous, 
equations  using  the  technique  of  Gaussian  elimination.  It  is  used  to 
solve  for  coefficients  of  the  fitted  surface. 

• SIGNF  - A surface  which  is  represented  by  an  equation  in  terms  of  its 
coordinates  can  also  be  represented  by  the  negative  of  the  equation.  In 
other  words, 

Ax+By+Cz+D=0=  -Ax  - By  - Cz  - D. 

In  one  case,  the  progrcun  will  generate  an  outward  normal  to  the  canopy 
while  in  the  other  case  it  will  generate  an  inward  normal.  In  subsequent 
calculations,  the  program  requires  that  the  normal  be  outward.  Therefore, 
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SIGNF  is  used  to  determine  correct  sign  on  coefficients.  It  does  this 
using  tests  on  the  two  possible  normal  vectors. 


• CFITS  - A set  of  boundary  points  representing  a reflective  surface  is 
used  to  produce,  in  general,  a quadratic  equation  (a  plane  surface  is 
considered  a degenerate  quadratic  surface) . This  equation  is  an  analyti- 
cal representation  of  the  points  in  a least-squares  error  sense.  There 
are  three  options  available  when  trying  to  match  an  equation  to  a set  of 
points.  The  first  is  to  specify  the  coefficients,  either  all  or  some. 

In  this  case  the  calculated  values  are  overridden  by  the  input  values. 

The  second  option  is  to  specify  which  coefficients  are  to  be  included  in 
the  matching  equation.  For  instance,  if  the  program  user  wished  the  sur- 
face through  the  boundary  points  to  be  quadratic  in  x while  linear  in  z, 
he  would  specify  through  his  inputs  that  the  other  terms  in  the  general 
quadratic  equation  be  ignored.  The  third  option,  which  is  probably  the 
most  useful,  starts  out  by  trying  to  match  the  points  to  a linear 
equation  (plane  surface) . If  this  equation  fails  to  satisfy  the  error 
criteria,  the  program  continues  to  form  quadratic  combinations.  The 
combination  with  the  minimum  least-squares  error  is  selected. 

• GLINT  - Calculates  the  sun  glints  for  each  reflective  surface  as  air- 
craft varies  heading.  First  the  components  of  the  sun  vector  are  trans- 
formed from  the  inertial  axes  to  the  body  axes  by  the  transformation 

in  Figure  4.  Then  each  boundary  point  of  the  reflective  surfaces  is 
checked  to  see  if  it  lies  in  the  shadow  of  a fence.  Next  the  normal 
to  the  reflective  surface  at  each  boundary  point  is  calculated  followed 
by  the  calculation  of  the  reflection  vectors.  Each  reflection  vector 
is  checked  to  see  if  it  intersects  any  fences.  Then  the  components 
of  the  reflection  vectors  are  transformed  from  the  body  axes  to  the 
inertial  euces.  Finally,  these  reflection  vectors  are  extended  as  in 
Figure  5 to  find  the  coordinates  of  each  sun  glint. 

• SHADOW  - Answers  the  question,  "Is  a boundary  point  of  a reflective 
surface  in  the  shadow  of  a fence?"  To  answer  this  question,  first  it 
is  determined  whether  any  of  the  planes  containing  the  fences  are 
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located  between  the  boundary  point  and  the  sun.  If  a fence  plane  lies 
between  the  sun  and  the  reflective  surface  boundary  point,  the  inter- 
section of  the  line  formed  by  projecting  the  boundary  point  along  the 
sun  vector  with  the  plane  of  the  fence  is  saved.  Finally,  it  is  deter- 
mined whether  any  of  the  saved  intersections  lie  within  the  boundaries 
of  their  respective  fences.  If  any  of  them  do,  then  the  reflective  sur- 
face boundary  point  is  indeed  in  a shadow  and  the  ISHAD  flag  is  set  to 
one  for  this  point.  The  coordinates  of  the  sun  glint  for  this  point 
are  still  calculated  and  tabulated  but  not  plotted  and  not  used  in  the 
probability  calculation. 

• PLANPT  - Given  the  equation  of  a plane  surface  and  the  representation  of 
a line  by  a point  and  direction  cosines,  the  intersection  of  the  line 
with  the  plane  is  calculated. 

• SORNOT  - Given  a point  and  the  coordinates  of  a unit  vector,  a line 
through  the  point  with  the  same  direction  cosines  as  the  vector  is  formed. 
It  is  then  determined  whether  a second  point  on  the  line  is  in  the  posi- 
tive or  negative  sense  of  the  vector  (see  Figure  10) . 

• INTFER  - Answers  the  question,  "Does  a reflection  vector  intersect  a 
fence  before  reaching  a possible  observer?"  To  answer  this  question, 
first  it  is  determined  whether  any  of  the  planes  containing  the  fences 
are  located  between  the  reflective  surface  boundary  point  from  which 
the  reflection  vector  eminates  and  any  possible  observer.  If  a fence 
plane  lies  between  the  boundary  point  and  the  observer,  the  intersection 
of  the  line  formed  by  extending  the  reflection  vector  with  the  plane  of 
the  fence  is  saved.  Finally,  it  is  determined  whether  any  of  the  saved 
intersections  lie  within  the  boundaries  of  their  respective  fences.  If 
any  of  them  do,  then  the  reflection  vector  does  indeed  strike  a fence 
and  the  INTFER  flag  is  set  to  one  for  the  reflective  boundary  point 
corresponding  to  the  reflection  vector.  The  coordinates  of  the  sun 
glint  for  this  point  are  still  calculated  and  tabulated  but  not  plotted 
and  not  used  in  the  probeibility  calculation. 
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• QUADPT  - Calculates  the  intersection  or  intersections  of  a line  with  a 
quadratic  surface.  Inputs  to  the  subroutine  are  the  coordinates  of  a 
point  on  the  line,  the  direction  cosines  of  the  line,  and  the  coefficients 
of  the  equation  describing  the  quadratic  surface. 


• MINMAX  - Determines  minimum  and  maximum  values  of  the  elevations  and 

azimuths  for  each  sun  glint.  At  this  point,  0.25  degree  is  added  to  the 
maximums  and  subtracted  from  the  minimums  to  represent  sun  dispersion 
angle.  These  values  are  plotted  and  used  to  determine  probcibility . 


• PROBL  - As  the  aircraft  changes  heading , with  fixed  sun  elevation  and 
fixed  aircraft  pitch  and  roll  attitudes,  the  sun  glint  from  a reflective 
surface  may  change  position  on  the  sun  glint  signature  plot,  as  was 
shown  in  Figure  6.  If  all  intermediate  positions  between  the  discrete  sun 
glint  positions  calculated  are  assiamed  possible,  then  there  is  a continu- 
ous area  swept  out- during  the  heading  changes.  The  total  of  all  the  areas 
swept  by  all  the  reflective  surfaces,  but  counting  any  overlap  only  once, 
divided  by  the  total  area  possible  will  yield  the  probability  for  this 
sun  elevation.  In  programing  a method  to  calculate  the  swept  area,  cer- 
tain assumptions  and  simplifications  were  found  necessary.  Figure  11  shows 
a hypothetical  example  of  the  motion  of  a sun  glint  from  one  reflective 
surface.  The  four  points  shown  for  each  sun  glint  are  the  results  from 
MINMAX,  which  represent  the  minimum  and  maximum  observer  elevations  and 
azimuths  for  this  glint.  This  is  simplification  number  one,  to  use  only 
the  minimum  and  maximums  to  represent  each  glint.  This  allows  for  the 
generation  of  the  polygon  shown  in  Figure  11  by  the  dotted  lines.  The 
single  point  at  zero  aircraft  heading  is  another  simplification  or 
assumption  and  is  caused  when  all  the  boundary  points  of  the  reflective 
surface  at  this  heading  are  in  the  shadow  of  one  or  more  fences  or  when 
the  generated  reflection  vectors  are  blocked  by  one  or  more  fences.  In 
other  words,  there  is  no  glint  for  this  surface  at  zero  heading.  However, 
at  headings  of  ±80  degrees  there  are  glints;  therefore,  the  approach 
chosen  was  to  place  a single  point  at  the  average  azimuth  and  elevation 
of  the  observers  had  the  fences  not  been  there  and  thus  a glint  created. 
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By  doing  this, the  area  is  tapered  to  zero  on  both  sides  of  the  ooint  (c) 
As  shovm,  the  points  (f)0@(c)@@(N)(M)@(g) 

® ® when  connected  by  straight  lines  in  this  order  form  a 

polygonal  shape. 

A polygon  will  be  generated  for  the  glint  motion  of  each  reflective  sur- 
face. The  overlap  must  now  be  considered.  Figure  12  shows  two  polygons 
representing  the  areas  swept  by  two  different  surfaces.  The  shaded 
area  represents  the  common  area  to  both  polygons.  The  program  treats 
polygon  0 as  if  polygon  0 does  not  exist;  thus  the  shaded  area  is 
included  in  0.  Then  polygon  0 is  broken  up  into  two  polygons  formed 
by  the  points  © 0 @ ® ® and  (0  ©00®© 

When  finished  there  are  three  polygons,  instead  of  two,  but  with 
no  common  area. 

• OVRLAP  - Inputs  to  the  routine  are  two  sets  of  points  representing  two 
polygons.  One  set  of  points  is  designated  the  shadow,  while  the  other 
set  is  the  surface.  It  is  the  purpose  of  this  routine  to  generate , if 
necessary,  a new  set  or  sets  of  points  which  represent  the  portion  or 
portions  of  the  surface  which  lie  outside  of  the  shadow.  This  routine 
is  used  in  probability  calculation  to  account  for  overlap.  Referring 
to  Figure  12,  the  polygon  representing  the  swept  area  from  surface 
number  one  would  be  the  shadow  while  the  polygon  representing  the 
swept  cirea  from  surface  number  two  would  be  the  surface.  In  this  case, 
the  output  from  OVRLAP  would  be  two  sets  of  points,  representing  the 
two  portions  of  surface  lying  outside  of  shadow. 

• INTERC  - Finds  intersections,  if  any,  of  a line  between  two  points  in 
a plane  with  the  set  of  lines  given  by  a set  of  connected  points  in 
the  same  plane. 

• CLOSE  - Determines  which  point  in  a set  of  points  in  a plane  is  closest 
to  a given  point  in  the  same  plane. 

• CHECK  1 - Determines  whether  a given  point  is  already  contained  in  a 

set  of  points. 
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• CHECK  2 - Similar  to  CHECK  1,  except  it  automatically  checks  a number 
of  sets  for  common  point. 

• CHECK  3 - Checks  to  see  whether  the  X or  Y values  of  a point  in  a plane 
are  outside  of  a given  range. 

• CHECK  4 - During  the  overlap  process,  it  becomes  necessary  to  check  to  see 
if  cuiy  lines  from  one  polygon  lie  on  top  of  any  lines  from  the  other 
polygon.  If  so,  one  of  the  lines  is  moved  slightly.  It  became  necessary 
to  add  this  routine  to  avoid  certain  numerical  problems  in  the  subroutine 
OVRLAP. 

• CIRCI£  - Used  to  determine  if  a given  point  in  a plane  is  enclosed  by 

a polygon  created  by  connecting  a given  set  of  points.  First,  a trans- 
lation is  performed  to  make  the  origin  of  the  axis  sytem  lie  on  top  of 
the  given  pnaint.  With  this  done,  the  number  of  positive  or  negative 
X-axis  crossings  is  counted.  If  the  number  is  odd,  the  point  is 
enclosed;  if  the  number  is  even,  the  point  is  not  enclosed.  This  is 
shown  by  two  examples  in  Figure  13. 

• AREA  - Determines  the  area  enclosed  by  a polygon.  Figure  14  shows  the 
method  used. 

• SORT  - If  part  of  the  polygon  representing  glint  sweep,  or  part  of  the 
glint  itself,  is  outside  of  observer  limits,  new  points  at  the  inter- 
sections with  the  boundaries  will  replace  those  outside  limits.  This 
is  done  for  the  probcibility  calculation  and  for  the  CALCOMP  plots. 

• S0RT2  - There  is  a natural  discontinuity  at  the  +180-degree  azimuths 

on  the  sun  glint  signature  plot.  When  a glint  is  wrapped  around  the 
180-degree  azimuth,  or  -180-degree  azimuth,  as  shown  in  Figure  15,  the 
points  ® ® (E)  of  the  glint  will  be  plotted  close  to  the  -180- 

degree  azimuth  line,  while  the  points  ® ® will  be  plotted  close 

to  the  +180-degree  line.  The  program  will  assume  that  the  glint  ® 
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POINT  ENCLOSED 


NUMBER  OF  POSITIVE  X CROSSINGS  IS  1 (ODD) 
NUMBER  OF  NEGATIVE  X CROSSINGS  IS  3 (ODD) 


Y 


NUMBER  OF  POSITIVE  X CROSSINGS  2 (EVEN) 
NUMBER  OF  NEGATIVE  X CROSSINGS  2 (EVEN) 


Figure  1 3.  Determination  of  Point  Enclosure. 
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Figure  14.  Area  Enclosed  by  Polygon. 


OBSERVER 

ELEVATION 


® © ® © is  connected  as  shown  with  the  solid  lines.  To  h^mdle 

this  situation  properly,  the  glint  ® (g)  @ @ (e)  is  replaced 

with  two  glints,  @ ® © ® © and  © ® © ® . 

The  new  glints  are  shown  shaded  in  Figure  15.  The  method  employed  is 
to  first  add  360  degrees  to  points  ® ® (e)  , moving  them  to  the 
right  as  shown.  Subroutine  SORT  is  then  used  to  find  points  ® © . 

Next,  360  degrees  are  subtracted  from  points  © © , moving  them  to 
the  left  as  shown.  Subroutine  SORT  is  then  used  to  find  points  © ® . 


EQUATIONS  USED  IN  THE  PROGRAM 

The  mathematical  e;Q)ressions  used  in  calculating  the  sun  glint  coordinates 
are  presented  in  this  section. 

The  general  form  of  the  equation  used  to  represent  all  fences  is  a plane, 

Ax  + By  + Cz  + D = 0. 

The  general  form  of  the  equation  used  to  represent  all  reflective  surfaces 
is  a quadratic  surface, 

2 2 2 

Ax  + By  +Cz  +Dxy  + Exz  + Fyz  +Gx+Hy+Kz+L=0. 

The  sun  vector  referred  to  inertial  coordinates  is  of  the  form 

S = ai  + bj  + ck  , 

where  (i,  j,  k)  are  unit  vectors  in  the  direction  of  the  inertial  x,  y,  z 
cixes,  respectively.  The  coefficients  a,  b,  and  c are  defined  using  y,  the 

sun  elevation,  as 

a = -cos  Y 
b = 0 

c = -sin  Y • 

The  normal  to  the  reflective  surface  in  body  axes  is 

d(^  3^  3^ 

N=N  I+N  3+N  k=  — i+~j+7~k, 

B B B B B B B 3x  B 8y  B 3z  B 
X y z 

where  (i  , j k ) are  unit  vectors  in  the  direction  of  aircraft  body  x, 

® ® 3(4  3(J>  ^ 3(t  ^ 

y,  z axes,  respectively,  and  the  coefficients  — , — , and  — are  the  first 

3x  3y  3z 

partial  derivatives  of  the  equation  representing  the  surface. 
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9(J) 

~ = 2Ax  + Dy  + Ez  + G 
34) 

— = 2By  + Dx  + Fz  + H 
3y 

34) 

— = 2Cz  + Ex  + Fy  + K. 

3z 

Before  the  reflection  vector  can  be  calculated,  the  sun  vector  must  be 
represented  in  terms  of  body  axes. 


T T 

(S  , S , S ) = D (-cos  Y,  0,  -sin  y) 

B ^B  B 

where  (S  , S , S ) are  the  components  of  the  sun  vector  in 
*B  ^B  ^B 

D is  the  transformation  matrix  shown  in  Figure  4. 


body  axes  and 


The  magnitude  of  the  sun  vector  along  the  normal  to  the  surface  is  the 
projection  of  the  sun  vector  onto  the  normal  vector.  Mathematically,  it 
is  the  dot  product  of  the  sun  vector  and  the  normal  vector,  divided  by  the 
magnitude  of  the  normal  vector. 


where  S=Si+S  j+S  k and 
B Xg  B yg  B Zg  B 


2 2 2 

+ N + N 
B B B 

X y z 


Using  the  value  of  S , the  reflection  vector  R is 

N B 


j+R  k=S  (2S)N  /In  |, 

B ’'b  B ^ ^B  ® ® ~ N B ' B ' 


Transforming  the  components  of  the  reflection  vector  from  body  to  inertial 
axes  is  done  by 


T 


D 


B 


T 
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where  D is  transpose  of  the  matrix  of  Figure  4. 

Each  boundary  point  of  the  reflective  surface  is  transformed  from  body 
coordinates  to  inertial  coordinates. 

T T T 

' ''s  % ’ 

III  B B,  B 

The  intersections  caused  by  the  line  passing  through  the  point  (x  , y , 

z ) with  direction  numbers  (R  , R , R ) and  the  cylinder 

2 2 2 
X + y = R , 


which  is  the  imaginary  surface  containing  the  possible  observers,  are 

‘V'  ^0/  ‘*0  ' "^0  ' 


1 1 


2 2 


X = (-b  + ]!  b^-4ac)  /2a 

°1 


X = (-b  - ||b^-4ac)/2a 


-X3^) 


°2  °1 


Z = Z + (R  /R  ) (Xq  -X  ) 
1 I I I 1 I 


°2  ’‘l  °1 


where  a = 1 + (R  /R  ) 

Yl  Xj 


b = 2y^  (R  /R  ) -2x  (R  /R  )‘ 
’'l  ®I  ^1  ^ 


The  equations  are  modified  slightly  if  R =0. 
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The  point  which  lies  in  the  positive  direction  of  the  reflection  vector  is 


designated  (x  , y , z ) , and  the  polar  coordinates  of  this  point  are 
o o o 


6 = tan 


where  a is  the  elevation  angle  and  B is  the  azimuth  angle  as  shown  in 
Figure  5. 


INPUT  VARIABLES 


The  input  variables  to  the  con5)uter  program  are  presented  iij  this  section. 

• TITLE  - Title  printed  at  top  of  each  output  page. 

• DEFLT  - Specifies  whether  default  option  is  to  be  used  during  inputting 
of  reflective  surface  data.  The  default  option  reduces  the  amount  of 
necessary  input. 

DEFLT  = 0 Default  option  not  used 
DEFLT  = 1 Default  option  used 

• PRINT  - Controls  how  much  and  what  type  of  data  will  be  printed.  Any 
plots,  whether  printer  plots  or  CALCOMP  plots,  will  be  controlled  by 
another  input.  Figure  16  shows  the  results  of  each  option. 

• PLOT  - Plot  options.  When  plots  are  requested,  a sun  glint  signature 
plot  for  each  sun  elevation  is  produced  either  on  CALCOMP  plotter  or  on 
printer,  or  both. 

PLOT  = 0 No  plots 

PLOT  = 1 CALCOMP  plots  only 

PLOT  = 2 Print  plots  only 

PLOT  = 3 Both  CALCOMP  and  print  plots 

• PROBL  - Specifies  whether  to  calculate  probability. 

PROBL  = 0 Probability  calculation  bypassed 
PROBL  = 1 Calculate  probability 

• THETA  - Aircraft  pitch  attitude  in  degrees.  Positi^?^'^for  nose  of  air- 

} 

craft  up.  Used  in  forming  the  inertial-to-body  axis  transformation 
shown  in  Figure  4. 
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PRINT  TABLE 
OF  GLINT 
COORDINATES 


Figure  16.  Computer  Output  From  Option  PRINT. 
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• PHI  - Aircraft  roll  attitude  in  degrees.  Positive  for  right  wing 
down.  Used  in  forming  the  inertial-to-body  axis  transformation  shown 
in  Figure  4. 

• GAMN  - Number  of  different  sun  elevations  to  be  analyzed.  Value  should 
be  less  than  or  equal  to  9.  If  not  specified,  the  program  assumes  a 
value  of  3. 

• GAMI  - Value  of  initial  sun  elevation  in  degrees.  Positive  for  sun 
ckbove  aircraft.  If  not  specified,  the  program  assumes  a value  of  zero. 

• DGAM  - Increment  by  which  each  successive  sun  elevation  is  increased. 
Value  is  in  degrees.  If  not  specified,  the  program  will  assume  a value 
of  30. 

Example  of  use  of  GAMN,  GAMI,  DGAM:  Suppose  it  is  desired  to  analyze 
three  sun  elevations,  0,  15,  and  30  degrees.  The  inputs  would  be 

GAMN  = 3.  GAMI  = 0.  DGAM  = 15. 

• PS IN  - Number  of  different  aircraft  headings  to  be  analyzed.  Value 
should  be  less  than  or  equal  to  9.  If  not  specified,  a value  of  9 
is  assumed  by  the  program. 

• PSII  - Value  of  initial  aircraft  heading  (yaw)  in  degrees.  Heciding  is 
positive  for  nose  right  as  shown  in  Figure  1 . A value  of  -180  degrees 
is  assumed  by  the  program  when  no  value  is  specified. 

• DPS I - Increment  by  which  each  successive  heading  is  increased.  Value 
is  in  degrees.  If  no  value  is  specified,  a value  of  45  degrees  is  used. 

• XROT  - Station  that  locates  the  center  of  rotation  of  the  aircraft  (see 
Figure  4).  Usually  the  aircraft  center  of  gravity.  If  not  specified, 
a value  of  zero  is  used. 

• YROT  - Buttline  location  of  the  center  of  rotation.  If  not  specified,  a 
value  of  zero  is  used. 
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• ZROT  - Waterline  location  of  the  center  of  rotation.  If  not  specified, 
a value  of  zero  is  used. 

• DISTG  - Radial  distance  (feet)  in  inertial  x-y  plane  from  origin  of 
inertial  axes  to  projected  glint.  This  distance  is  specified  as  R in 
Figure  5.  If  not  specified,  a value  of  1,000  feet  is  used  by  the  progrcun. 

• ELVMIN  - Lower  limit  on  elevation  angle  of  possible  observers.  Value 
is  in  degrees.  If  not  specified,  a value  of  -10  degrees  is  used. 

ELVJ'.T?a  is  used  in  probability  calculation  and  lower  boundary  on  sun 
glint  signature  plots.  Value  Ccinnot  be  less  than  -90  degrees. 

• ELVMAX  - Upper  limit  on  elevation  angle  of  possible  observers.  Value  is 
in  degrees.  If  not  specified,  a value  of  10  degrees  is  used.  ELVMAX 

is  used  in  probability  calculation  and  upper  boundary  on  sun  glint 
signature  plot.  Value  cannot  be  greater  than  90  degrees. 

• AZMMIN  - Minimum  allowable  azimuth  angle  to  observer,  in  degrees.  If 
not  specified,  a value  of  -180  degrees  is  used.  AZMMIN  is  the  boundary 
on  the  left  of  the  sun  glint  signature  plot.  Value  cannot  be  less  than 
-180  degrees. 

• AZMMAX  - Maximum  allowable  azimuth  angle  to  observer,  in  degrees.  If 
not  specified,  a value  of  180  degrees  is  used.  AZMMAX  is  the  boundary 
on  the  right  of  the  sun  glint  signature  plot.  Value  CeUinot  be  greater 
than  180  degrees. 

• FENCES  - Number  of  fences.  Meiximum  number  is  10. 

• PANELS  - Number  of  reflective  (transparency)  panels.  Maximum  numJaer 
is  12. 

• FPTS  - Number  of  boundary  points  on  a fence.  Maximum  value  is  10. 

• XFPT  - Station  for  a fence  boundary  point. 

• YPPT  - Buttline  for  a fence  boundary  point. 
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ZFPT  - Waterline  for  a fence  boundary  point. 


• TITLES  - Transparency  (reflective)  surface  identif icatioi). 

• SYMTRY  - Specifies  whether  the  transparency  surface  just  identified  is 
symmetric  to  the  surface,  whose  data  has  just  been  read  in.  This  option 
will  be  made  clearer  with  the  sample  cases  in  Appendix  A.  The  possible 
values  for  SYMTRY  are: 


SYMTRY  = 

0. 

Not  symmetric  to  previous  surface  and  therefore  data 

for  this  surface  is  necessary. 

SYMTRY  = 

1. 

Symmetry  is  about  x-z  plane. 

SYMTRY  = 

2. 

Syrmnetry  is  cibout  x-y  plane. 

SYMTRY  = 

3. 

Symmetry  is  about  y-z  plane. 

If  SYMTRY  ^0,  then  the  following  variables  are  not  necessary  for  this 
reflective  surface. 

The  variables  XFIT,  YFIT,  ZFIT,  XYFIT,  XZFIT,  YZFIT,  XSQFIT,  YSQFIT, 
ZSQFIT,  XCOF,  YCOF,  ZCOF,  XYCOF,  XZCOF,  YZCOF,  XSQCOF,  YSQCOF,  ZSQCOF, 
and  CTERM  are  used  only  if  DEFLT  = 0.  With  these  variables,  the 
program  user  may  either  specify  which  terms  in  the  general  quadratic 
equation  he  wishes  to  use  in  the  curve-fit  routine,  or  he  may  specify 
the  actual  values  for  these  coefficients. 

• XFIT  - Should  the  linear  x term  in  the  quadratic  equation  be  used  in 
the  curve-fit  analysis  of  this  surface? 

XFIT  = 0.  No 
XFIT  =1.  Yes 

• YFIT  - Should  the  linear  y term  in  the  quadratic  equation  be  used  in 
the  curve-fit  analysis  of  this  surface? 
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YFIT  =0.  No 
YFIT  = 1.  Yes 

• ZFIT  - Should  the  linear  z term  in  the  quadratic  equation  be  used  in 
the  curve-fit  analysis  of  this  surface? 

ZFIT  =0.  No 
ZFIT  = 1.  Yes 

• XYFIT  - Should  the  quadratic  xy  term  in  the  quadratic  equation  be 
used  in  the  curve-fit  analysis  of  this  surface? 

XYFIT  =0.  No 
XYFIT  = 1.  Yes 

• XZFIT  - Should  the  quadratic  xz  term  in  the  quadratic  equation  be  used 
in  the  curve-fit  analysis  of  this  surface? 

XZFIT  =0.  No 
XZFIT  = 1.  Yes 

• YZFIT  - Should  the  quadratic  yz  term  in  the  quadratic  equation  be  used 
in  the  curve-fit  analysis  of  this  surface? 

YZFIT  =0.  No 
YZFIT  = 1.  Yes 

2 

• XSQFIT  - Should  the  quadratic  x term  in  the  quadratic  equation  be  used 
in  the  curve-fit  analysis  of  this  surface? 

XSQFIT  =0.  No 
XSQFIT  = 1.  Yes 

2 

• YSQFIT  - Should  the  quadratic  y term  in  the  quadratic  equation  be  used 
in  the  curve-fit  analysis? 
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YSQFIT  =0.  No 
YSQFIT  = 1.  Yes 


• ZSQFIT  - Should  the  quadratic  z term  in  the  quadratic  equation  be  used 
in  the  curve-fit  analysis? 

ZSQFIT  =0.  No 
ZSQFIT  = 1.  Yes 

• XCOF  - If  a nonzero  value  is  read  in  for  XCOF,  then  the  value  for  XCOF 
will  be  used  for  the  coefficient  of  x in  the  quadratic  equation  repre- 
senting this  reflective  surface.  XFIT  must  be  set  to  one  when  XCOF  is 
used. 

• YCOF  - If  a nonzero  value  is  read  in  for  YCOF,  then  the  value  for  YCOF 
will  be  used  for  the  coefficient  of  y in  the  quadratic  equation  repre- 
senting this  reflective  surface.  YFIT  must  be  set  to  one  when  YCOF  is 
used. 

• ZCOF  - If  a nonzero  value  is  read  in  for  ZCOF,  then  the  value  for  ZCOF 
will  be  used  for  the  coefficient  of  z in  the  quadratic  equation  repre- 
senting this  reflective  surface.  ZFIT  must  be  set  to  one  when  ZCOF  is 
used. 

• XYCOF  - If  a nonzero  value  is  read  in  for  XYCOF,  then  the  value  for 
XYCOF  will  be  used  for  the  coefficient  of  xy  in  the  quadratic  equation 
representing  this  reflective  surface.  XYFIT  must  be  set  to  one  when 
XYCOF  is  used. 

• XZCOF  - If  a nonzero  value  is  read  in  for  XZCOF,  then  the  value  for 
XZCOF  will  be  used  for  the  coefficient  of  xz  in  the  quadratic  equation 
representing  this  reflective  surface.  XZFIT  must  be  set  to  one  when 
XZCOF  is  used. 
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XSQCOF  - If  a nonzero  value  is  read  in  for  XSQCOF,  the  value  for 

2 

XSQCOF  will  be  used  for  the  coefficient  of  x in  the  quadratic  equation 
representing  this  reflective  surface.  XSQFIT  must  be  set  to  one  when 
XSQCOF  is  used. 


• YSQCOF  - If  a nonzero  value  is  read  in  for  YSQCOF,  the  value  for 

2 

YSQCOF  will  be  used  for  the  coefficient  of  y in  the  quadratic  equation 
representing  this  reflective  surface.  YSQFIT  must  be  set  to  one  when 
YSQCOF  is  used. 


• ZSQCOF  - If  a nonzero  value  is  read  in  for  ZSQCOF,  the  value  for 

2 

ZSQCOF  will  be  used  for  the  coefficient  of  z in  the  quadratic  equation 
representing  the  reflective  surface.  ZSQFIT  must  be  set  to  one  when 
ZSQCOF  is  used. 

• CTERM  - If  a nonzero  value  is  read  in  for  CTERM,  the  value  for  CTERM 
will  be  used  for  the  constant  term  in  the  quadratic  equation  for  this 
reflective  surface. 


• RPTS  - Number  of  boundary  points  on  reflective  surface.  Maximum  value 
is  30. 


• XSPT  - Station  for  a reflective  surface  boundary  point. 

• YSPT  - Buttline  for  a reflective  surface  boundary  point. 


• ZSPT  - Waterline  for  a reflective  surface  boundary  point. 


Fifty-one  input  paraimeters  have  been  defined,  but  the  actual  number  of 
inputs  depends  on  the  particular  case.  The  arrangement  of  the  input 
parameters  on  data  cards  for  a configuration  consisting  of  two  fences  and 
two  reflective  surfaces  is  shown  in  Figure  17.  The  first  fence  has  three 
boundary  points  emd  the  second  fence  has  four.  The  two  reflective  sur- 
faces are  symmetric  about  the  x-z  plane  and  are  defined  by  four  boundeiry 


45 


COLUMN 


COLUMN 

•0 


FIRST 

REFLECTIVE 

SURFACE 


1 TITLE 

1 

1 

11 

21 

31 

41 

SI 

SQ 

60 

1 OEFLT* 

]FRINT 

1 THETA 

] PHI 

1 

1 

1 

1 

1 

11 

21 

31 

41 

61 

60 

60 

1 QAMN 

Igami 

1 OGAM 

1 PSIN 

|PSII 

1 OPS 

1 

1 

1 

11 

21 

31 

41 

61 

60 

60 

1 XROT 

1 VROT 

I ZROT 

1 

1 

1 

1 

1 

1 

11 

21 

31 

41 

61 

60 

60 

1 OISTG 

1 ELVMIN 

1 ELVMAX 

1 AZMMIN 

1 AZMMAX 

1 

1 

1 

1 

11 

21 

31 

41 

61 

60 

60 

1 FENCES 

1 PANELS 

|PLOT 

1 PflOBL 

1 

1 

1 

1 

1 

11 

21 

31 

41 

51 

60 

60 

T n □ 

1 

11 

21 

31 

41 

61 

60 

60 

FIRST  j 
FENCE 

1 XFFTI1I 

1 VFPTMI 

1 ZFPTdl 

1 XFPTI3I 

1 VFPTI2I 

( ZFPTI2I 

1 

1 

1 

11 

21 

31 

41 

61 

60 

60 

1 XFPTI3I 

1 VFPTI3I 

1 ZFPTOI 

1 

1 

1 

1 

1 

, 

11 

21 

31 

41 

61 

60 

60 

iFFTs  1 1 n r— n n 

SECOND 
FENCE  ] 

1 

11 

21 

31 

41 

61 

60 

60 

1 XFFTMI 

1 YFFTI1) 

1 ZFFT(I) 

1 XFPT(2I 

1 YFPTIZI 

1 ZFPT(2) 

1 

1 

1 

11 

21 

31 

41 

61 

60 

60 

1 XFFTO) 

1 YFPTOI 

1 ZFPTO) 

1 XFPTI4) 

1 YFPT(4) 

1 ZFPT(4) 

1 

1 

1 

Ifi 

1 TITLES 

1 

1 

1 1 

21 

31 

41 

61 

60 

60 

1 SVMTRY* 

1 

1 

1 

1 

1 

1 

1 

11 

21 

31 

41 

61 

60 

60 

1 XFIT 

1 YFIT 

1 ZFIT 

1 XYFIT 

1 XZFIT 

1 VZFIT 

1 

1 

•• 

1 

11 

21 

31 

41 

61 

60 

60 

TH  ESE 
CAROS 

1 XSOFIT 

1 YSQFIT 

(ZSQFIT 

1 

1 

1 

1 

1 

OMITTED 

IF  1 

1 

11 

21 

31 

41 

61 

60 

60 

DEFAULT 

orrioN 

USED 

1 XCOF 

1 VCOF 

1 ZCOF 

1 XYCOF 

[XZCOF 

1 YZCOF 

"1 

1 

1 

11 

21 

31 

41 

61 

60 

60 

1 XSQCOF 

1 YSQCOF 

1 ZSQCOF 

1 CTERM 

1 

1 

1 

t 

11 

21 

31 

41 

61 

60 

60 

1 RPTS 

1 

1 

1 

1 

1 

1 

1 

11 

21 

31 

41 

61 

60 

60 

1 XS^TIII 

1 Ysrrit) 

] ZSPT(l) 

1 XSPT(2) 

1 YSPT(2» 

1 ZSPT(2) 

1 

1 

1 

11 

21 

31 

41 

61 

60 

60 

1 XSFTOI 

1 YSFTO) 

1 ZSPTO) 

1 XSPT(4I 

1 VSPT(4) 

1 ZSPT<4) 

1 

1 

1 

11 

21 

31 

41 

61 

60 

60 

1 XSFTIS) 

1 YsrriB) 

1 ZSPTIS) 

I 

1 

1 

“I 

1 

1 

SECOND 

(titles  I 

reflective  • 

' 

_1J 



31 

41 

61 

60 



SURFACE 

^ 1 1 1 1 I 

FOR  THISSETUF 
OEFLT-0 

SVMTRY-0  (FIRST  reflective  SURFACE) 


SVMTRY  • 1 (SECOND  reflective  SURFACE) 

••  oeflt  - t 


Figure 


7.  Input  Data  Scheme. 


46 


points.  The  default  option  will  not  be  used  in  readinq  in  the  reflective 
surface  data;  however,  the  continents  show  which  cards  are  omitted  if  the 
default  option  is  used. 


GENERATED  SYMBOLS 


The  symbolic  names  for  the  variables  of  interest  which  are  assigned  values 

during  the  execution  of  the  progrcun  are  presented  in  this  section. 

• COEFFB  - Array  containing  the  coefficients  of  the  equations  represent- 
ing the  fences. 

• CX3EFSB  - Array  containing  the  coefficients  of  the  equations  represent- 
ing the  reflective  surfaces. 

• SIGF  - Array  containing  the  root-mean-square  estimate  of  the  error  for 
each  fence  curve-fit. 

• SIGS  - Array  containing  the  root-mean-square  estimate  of  the  error  for 
each  reflective  surface  curve-fit. 

• GAMMA  - Sun  elevation  angle.  Takes  on  values  specified  by  the  inputs 
GAMN,  GAMI,  and  DGAM. 

• ALPHA2  - Array  containing  the  elevation  angles  for  all  the  sun  glints 
generated . 

• BETA2  - Array  containing  the  azimuth  angles  for  all  the  sun  glints 
generated. 

• INTRFL  - Array  containing  a set  of  condition  flags.  Each  flag  records 
whether  a particular  set  of  sun  glint  coordinates  was  generated  by  an 
internal  reflection. 

INTRFL  =■  0 External  reflection 
INTRFL  “ 1 Internal  reflection 


• ISHADW  - Array  containing  a set  of  condition  flags.  Each  flag  records 
whether  a particular  set  of  sun  glint  coordinates  was  generated  from  a 
boundary  point  in  the  shadow  of  a fence. 

ISHADW  = 0 Boundary  point  not  in  shadow 
ISHADW  = 1 Boundary  point  in  shadow 

• INTERF  - Array  containing  a set  of  condition  flags.  Each  flag  records 
whether  a particular  set  of  sun  glint  coordinates  was  generated  by  a 
reflection  vector  which  intersected  a fence. 

INTERF  =0  No  intersection  with  fence 

INTERF  = 1 Reflection  vector  blocked  by  a fence 

• APLOT  - Array  containing  the  minimum  and  maximum  elevation  angles  which 
were  formed  at  each  aircraft  heading. 

• BPLOT  - Array  containing  the  minimum  and  maximum  azimuth  angles  which 
were  formed  at  each  aircraft  heading. 

• SUN  - Three-dimensional  array  containing  the  coordinates  of  the  unit 
sun  vector  in  inertial  auces. 

• SUNB  - Three-dimensional  array  containing  the  coordinates  of  the  unit 
sun  vector  in  aircraft  body  axes. 

• DIRCOS  - Three-by- three  matrix  containing  the  orthonormal  tremsformation 
shown  in  Figure  4 used  to  transform  coordinates  of  a vector  frcxn  inertial 
to  body  axes. 

• 2®  ~ Coordinates  of  a reflective  surface  boundary  point  in  body  cixes . 
Array  is  redefined  for  each  boundary  point  being  analyzed. 

• - Coordinates  of  a reflective  surface  boundary  point  in  inertial  axes. 

• VNORMB  - Coordinates  of  normal  to  reflective  surface. 
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• RFLTN  - Coordinates  of  reflection  vector  referred  to  inertial  axes. 

• i^CPLOT  - Array  of  azimuthal  values  of  points  forming  the  polygon  repre- 
*^senting  the  area  swept  by  glints  from  a reflective  surface. 

• YPLOT  - Array  of  elevation  angles  for  points  forming  the  polygon  repre- 
senting the  area  swept  by  glints  from  a reflective  surface. 

• SWEEP  - Area  swept  by  the  sun  glints  from  all  the  reflective  surfaces, 
ignoring  overlap. 

•'  SWA2  - Area  swept  by  the  sun  glints  from  all  the  reflective  surfaces, 
counting  the  overlap  area  only  once. 

• TAREA  - Total  possible  area  of  sun  glint  signature  plot. 

• VPROBL  - Value  of  probability  calculation.  SWA2  divided  by  TAREA. 
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OUTPUT 


The  output  from  the  program  is  broken  up  into  five  sections  and  controlled 
by  the  input  variables  PRINT  and  PLOT  as  defined  previously.  The  five 
types  of  output  are: 

• Input  data 

• Curve-fit  results 

• Tabular  suh  glint  results 

t 

• Sun  glint  signature  plots  using  printer 

• Sun  glint  signature  plots  using  CALCOMP  plotter. 

The  printout  of  the  input  data  is  just  a feedout  of  the  data  read  in  from 
cards.  Sheet  1 of  B'igure  18  shows  the  first  page  of  this  type  output.  The 
first  group  of  input  data  is  printed  out  under  the  caption  **CONTROL 
OPTIONS.  The  input  parameters  presented  are  DEFLT,  PRINT,  PLOT,  and  PROBL, 
respectively.  The  second  group  of  input  data  comes  under  the  caption 
♦‘AIRCRAFT  INITIAL  EULER  ORIENTATION.  The  input  parameters  printed  are 
THETA,  PHI,  AND  PSII.  Under  the  third  and  fourth  captions,  “SELECTED  YAW 
ANGLE  ROTATIONS  and  **SELECTED  SUN  ELEVATIONS,  are  printed  the  values  for 
the  variables  PS  IN,  PSII,  DPSI,  and  GAMN,  GAMI,  DGAM,  respectively.  The 
fifth  set  of  information,  “REFERENCE  INFORMATION  FOR  MEASURING  SUN  GLINT 
ANGLES,  in  the  printout  of  XROT,  YROT,  ZROT,  and  DISTG.  The  sixth  and 
seventh  sets  of  information,  “PLOT  SCALING,  and  “BOUNDARIES  FOR  CALCU- 
LATING probability,  are  printouts  of  the  same  variables,  AZMMIN,  AZMMAX, 
ELVMIN,  and  ELVMAX.  They  are  printed  out  twice,  under  different  titles,  to 
show  that  the  same  values  are  used  in  both  places.  Sheet  2 of  Figure  18 
is  a printout  of  the  fence  input  data.  The  number  of  fences  (FENCES) , the 
number  of  boundeury  points  (FPTS) , and  the  location  of  each  boundary  point 
(XFPT,  YFPT,  ZFPT)  are  printed  out.  ■ Sheet  3 of  Figure  18  shows  how  the 
input  data  for  a reflection  surface  is  output.  Under  ♦♦REFLECTIVE  SURFACES 
DATA,  the  number  of  reflective  panels  (PANELS)  is  printed.  Then,  for  each 
reflective  surface,  cin  identification  title  (TITLES) , preselected  curve-fit 
terms  (XFIT,  YFIT,  ZFIT,  XYFIT,  XZFIT,  YZFIT,  XSQFIT,  YSQFIT,  ZSQFIT) , 
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Figure  18.  Typical  Printout  of  Input  Data  (Sheet  1 of  3), 


preselected  curve-fit  coefficients  (XCOF,  YCOF,  ZCOF,  XYCOF,  XZCOF,  YZCOF, 
XSQCOF,  YSQCOF,  ZSQCOF,  CTERM) , and  boundary  points  are  printed. 

The  second  type  of  output  is  the  result  of  curve-fitting  the  boundary  points 
of  each  fence  and  reflective  surface.  Sheet  1 of  Figure  19  shows  the  co- 
efficients of  the  plane  best  approximating  the  boundary  points  of  the  fence. 
For  the  case  shown,  the  equation  would  be  0.01476y  = 0 or  just  y = 0,  which 
is  the  equation  of  the  x-z  plane.  Sheet  2 of  Figure  19  shows  the  coeffi- 
cients of  several  quadratic  equations  approximating  several  reflective 
surfaces.  For  the  case  shown,  the  equation  for  the  first  surface  is  a 
plane,  -0.006923x  + 0.005206z  -1  = 0. 

The  third  type  of  output  is  the  tabular  sun  glint  results.  Each  table 
represents  the  results  for  one  reflective  surface  at  constant  sun  elevation 
as  aircraft  heading  changes.  Part  of  a typical  table  is  shown  in  Figure 
20.  At  the  top  of  the  page,  the  reflective  surface  is  identified.  Next, 
the  sun  elevation,  aircraft  pitch  and  roll  attitude,  location  of  aircraft 
center  of  rotation,  and  radial  distance  to  observers  are  printed.  Then, 
for  each  aircraft  heading,  each  boundary  point  is  listed  along  with  three 
condition  flags  and  the  point  on  the  sun  glint  resulting  from  the  reflec- 
tion from  this  boundary  point.  The  condition  flags  INRFL,  ISHAD,  and 
INTRF  specify  whether  the  reflection  from  this  boundary  point  is  internal, 
whether  the  boundary  point  lies  in  a shadow,  or  whether  the  reflection 
vector  from  this  boundary  point  intersects  a fence.  A value  of  zero  means 
no  while  a one  means  yes . 

The  fourth  type  of  output  is  a print  plot  of  a sun  glint  signature,  as 
shown  in  Figure  21.  Each  plot  is  for  a constant  sun  elevation.  The  values 
that  are  plotted  are  the  minimum  and  maximum  azimuths  and  elevations  for 
each  sun  glint  resulting  from  each  reflective  surface  at  each  aircraft 
heading  analyzed.  Each  letter,  as  the  legend  below  the  plot  shows,  is  for 
a constant  aircraft  heading  and  may  have  resulted  from  more  than  one  sur- 
face. The  probability  as  previously  defined  is  also  shown  in  the  right- 
hand  corner  below  the  plot. 
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Figure  19.  Typical  Printout  of  Curve-Fit  Results  (Sheet  1 of  2) 


Figure  21 . Typical  Printout  of  On-Line  Plot. 


The  fifth  and  final  type  of  output  comes  from  a CALCOMP  plotter.  It  is,  as 
can  be  seen  from  Figure  22,  a continuous  plot.  As  was  the  case  for  the 
print  plots,  each  CALCOMP  plot  is  the  total  sun  glint  signature  of  the  con- 
figuration at  a constant  sun  elevation.  The  aircraft  heetdings  are  shown 
on  the  plots  as  near  as  possible  to  the  sun  glint  they  represent.  However, 
when  more  than  one  surface  produces  a glint  in  the  same  area,  as  in  the 
bottom  plot,  the  values  are  sometimes  printed  on  top  of  each  other.  This 

nuisance  has  not  been  corrected  yet  but  should  cause  no  problems,  since  1 

the  numbers  can  still  be  distintinguished. 

The  plotter  should  be  started  with  the  pen  origin  at  0,0  or  Lower  Left.  1 


SON  CLCVATION  RNCLC  r SO- 
N/C  PITCH  BTT  = -4.0 
R/C  ROLL  BTT  : 0-0 


HfROlNO  OF  RIRCRRFT  (OEOJ 

UiO 


Uj'-l8C.00  -120-00  -60-00  0-00  60-00  120.00  180-00 

-i  OBSERVER  POS  fOEG) 


SON  CLCVRT13N  RNCLC  : 30. 

R/C  PITCH  RTT  r -4.0 
R/C  ROLL  RTT  : 0-0 

HCROINO  OF  RIRCRRFT  I0E5I 

tuo 


SON  CLCVRTION  RNCLC  = 0. 

R/C  PITCH  BTT  : -4.0 
R/C  ROLL  RTT  ; 0.0 

HCROINO  OF  RIRCRRFT  lOEOI 

UJO 


Figure  22.  Typical  CALCOMP  Plot 
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RUNN ING  INSTRUCTIONS 


The  program  is  written  in  FORTRAN  IV  computer  language.  It  has  been  suc- 
cessfully compiled  and  executed  on  the  Boeing  Computer  Service's  IBM  370-158 
at  Philadelphia,  Pennsylvania.  In  addition  to  the  program  itself,  the 
CALCOMP  Basic  Software  package  is  used  for  continuous  plots.  This  CALCOMP 
package  includes  six  routines  that  are  called  from  the  program  if  the 
input  variable  PLOT  = 1 or  3.  The  routines  are  PLOT,  SYMBOL,  NUMBER, 

SCA1£,  AXIS,  and  LINE.  These  are  the  only  CALCOMP  subroutines  required. 

The  following  is  a physical  setup  of  the  deck  required  to  compile  and 
execute  on  the  IBM  370-158  at  Boeing. 

//  JCL  cards  particular  to  IBM  370-158  at  Boeing  Computer  Services. 
Input  is  read  fran  unit  5 and  output  printed  on  unit  6. 

//EXEC  C0MP2G0 

//FTC. SYS IN  DD* 

(Program  goes  here.) 

//G0.FTO5FOO1  DD* 

(Data  deck  goes  here.) 

/* 

// 

At  Boeing  Vertol,  the  compiled  program  is  stored  on  a disk;  therefore,  once 
compiled,  any  subsequent  runs  require  only  the  proper  JCL  cards  and  the 
data  deck.  It  is  not  known  whether  the  user  has  such  a capability,  so  any 
reference  to  the  JCL's  used  at  Boeing  Vertol  could  be  meaningless. 

The  program  takes  about  35  seconds  to  compile  on  the  IBM  370.  Typical 
execution  times , depending  on  the  amount  of  output  requested  and  the  number 
of  conditions  analyzed,  are  approximately  8 to  12  seconds. 

As  for  the  form  of  the  input,  all  titles  are  read  using  an  alphanumeric 
format,  which  allows  for  the  mixing  of  letters  and  numbers  in  the  input. 

The  rest  of  the  inputs  are  read  using  a floating-point  format  with  10 
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columns  alloted  for  each  variable.  As  long  as  the  decimal  point  is  in- 
cluded, the  variable  may  appear  anywhere  in  those  10  columns.  The  specific 
location  of  each  variable  on  the  card  was  shown  in  Figure  17 . 


APPENDIX  A 


SAMPLE  CASES 


Two  sample  cases  will  be  discussed  to  demonstrate  the  use  of  the  program. 
The  first  example  deals  with  the  configuration  shown  in  Figure  A-1.  The 
panels  are  all  flat.  There  are  six  panels,  four  of  them  being  identified 
in  the  isometric.  The  right- forward  and  the  right-aft  panels  are  hidden 
from  view.  The  right-forward  panel  is  symmetric  to  the  left-forward  panel 
about  the  x-z  plane  of  the  aircraft.  Likewise  the  right-aft  panel  is  sym- 
metric to  the  left-aft  panel  about  the  x-z  plane  of  the  aircraft.  There 
are  two  fences  on  the  canopy,  one  over  the  windshield  while  the  other  is 
over  the  top  panel.  Since  the  surfaces  are  all  flat,  as  few  as  three 
points  can  be  used  to  locate  the  boundaries  of  the  surface.  The  circled 
letters  represent  the  boundary  points  selected.  In  practice,  more  boundary 
points  may  be  selected;  however,  in  most  cases  of  flat  panels  a large 
number  of  boundary  points  is  not  necessary.  Using  the  three-view  drawing 
shown,  the  center  of  rotation  and  each  circled  letter  can  be  located  refer- 
enced to  the  station,  buttline,  waterline  axes.  The  panels  and  fences 
can  then  be  described  by  a set  of  locations  as  tabulated  in  Table  A-1.  For 
instance,  the  windshield  is  defined  by  the  locations  of  @ @ @ @ . 

These  are  the  boundary  points  that  are  used  by  the  program. 

With  the  panels  and  fences  defined  and  the  location  of  the  center  of  rota- 
tion known,  the  next  step  is  to  decide  on  what  conditions  to  analyze. 

Using  Figure  17  as  a guide,  the  first  input  card  would  contain  the  title 
for  this  run.  The  title  chosen  is  FLAT  PANELS  WITH  FENCES  SAMPLE  CASE  I, 
which  will  appear  on  the  first  card.  The  second  card  contains  the  input 
variables,  DEFLT,  PRINT,  THETA,  and  PHI.  The  default  option  is  chosen 
(DEFLT  = 1.);  therefore  the  optional  cards  shown  in  Figure  17  will  be 
omitted.  Since  all  the  printout  possible  is  desired,  PRINT  will  be  set 
equal  to  2.  The  aircraft  pitch  and  roll  attitudes  will  be  assumed  to  be 
zero  (THETA  = 0.  and  PHI  = 0.) . This  takes  care  of  the  second  card.  The 
third  card  will  specify  the  desired  sun  elevations  and  aircraft  headings 
(GAMN,  GAMI,  DGAM,  PSIN,  PSII,  DPSl) . Three  sun  elevations  are  chosen, 

0,  30,  and  60  degrees.  Therefore,  GAMN  will  be  equal  to  3.,  GAMI  will  be 
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TABLE  A-1.  COOEUDINATES  OF  CASE  I HELICOPTER  CONFIGURATION 


Coordinates  of  Center  of  Rotation  (CG) 


Station 


Buttline 


Waterline 


Coo’^dinates  of  Windshield  Boundary  Points 


Station 


Buttline 


Waterline 


Coordinates  of  Top  Panel  Boundary  Points 


Station 


Buttline 


Waterline 


118.8 

118.8 


Coordinates  of  Left  Forward-Panel  Boundary  Points 


Station 


Buttline 


Waterline 


Coordinates  of  Left-Aft  Panel  Boundary  Points 


Station 


Buttline 


Waterline 


118.8 

118.8 


Coordinates  of  Fence  No.  1 (Over  Windshield)  Boundary  Points 


Station 


Buttline 


Waterline 


22.32 


TABLE  A-1  - Continued 


r 


Coordinates 

of 

Fence 

No.  1 - Continued 

Station 

Buttline 

Waterline 

© 

54. 

0 

72. 

27.6 

0 

42. 

Coordinates 

of 

Fence 

No.  2 (Over  Top  Panel) 

Boundary  Points 

Station 

Buttline 

Waterline 

© 

54. 

0 

78. 

© 

1.18.8 

0 

78. 

118.8 

0 

72. 

© 

54. 

0 

72. 

0.,  and  DGAM  will  be  30.  Aircraft  headings  shall  cover  the  whole  360-degree 
range;  therefore,  PSIN  will  be  equal  to  9.,  PSII  = -180.,  and  DPSI  = 45. 
These  six  values  will  take  care  of  the  third  card.  The  fourth  card  contains 
the  location  of  the  center  of  rotation  (XROT,  YROT,  ZROT) . These  values 
are  measured  on  Figure  A-1  and  recorded  at  the  top  of  Table  A-1.  The  values 
are  XROT  = 150.,  YROT  = 0.,  and  ZROT  = 36. 

This  is  all  the  information  needed  for  the  fourth  card.  The  fifth  card 
contains  the  radial  distance  to  all  observers  (DISTG) , the  minimum  (ELVMIN) 
and  maximum  (ELVMAX)  elevation  angles,  and  the  minimum  (AZMMIN)  and  maxi- 
mum (AZMMAX)  azimuth  angles.  The  distance  to  the  observer  will  be  assumed 
to  be  1,000  feet.  The  elevation  limits  on  the  observers  will  be  ;tl0 
degrees  and  the  azimuth  limits  will  be  j^l80  degrees.  Therefore  the  fift'. 
card  has  DISTG  = 1,000  feet,  ELVMIN  = -10  degrees,  ELVMAX  = 10  degrees, 
AZMMIN  = -180  degrees  and  AZMMAX  = 180  degrees.  The  sixth  card  contains 
the  number  of  fences  (FENCES) , the  number  of  reflective  panels  (PANELS) , 
the  plot  option  chosen  (PLOT) , and  whether  to  calculate  probability.  For 
this  configuration,  FENCES  = 2.  and  PANELS  =6.  To  get  both  printer-type 
and  CALCOMP-type  plots,  set  PLOT  = 3.,  and  to  have  probability  calculated, 
PROBL  = 1.  This  takes  care  of  the  control  cards.  The  next  six  cards  will 
contain  fence  data.  The  seventh  card  contains  the  number  of  boundary 
points  associated  with  the  first  fence  (over  the  windshield) . There  are 
four  points;  therefore,  on  the  seventh  card  FPTS  = 4.  The  eighth  and 
ninth  cards  contain  the  four  boundary  points  @ @ @ 

from  Table  A-1. 

The  tenth  card  has  the  value  for  the  number  of  boundary  points  defined  for 
the  second  fence  (over  the  top  panel) , FPTS  = 4.  The  eleventh  and  twelfth 
cards  contain  the  four  boundary  points  (2)  This  completes 

the  fence  data.  The  rest  of  the  data  cards  contain  information  about  the 
reflective  surfaces.  The  thirteenth  card  is  an  identification  card  for  the 
first  reflective  surface  read  in.  For  this  case  it  is  the  WINDSHIELD.  The 
fourteenth  card  is  the  symmetry  card  for  this  surface.  Since  this  is  the 
first  surface  and  therefore  data  must  be  read  in  for  it,  SYMTRY  = 0.  Since 
DEFLT  = 1.,  the  optional  cards  shown  in  Figure  17  pertaining  to  selecting 
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certain  terms  in  the  curve-fit  equation  will  be  omitted.  The  fifteenth 
card  contains  RPTS,  the  number  of  boundary  points  for  this  surface;  RPTS 
has  a value  of  4 . The  sixteenth  and  seventeenth  cards  contain  the  four 
boundary  points  for  this  surface,  @ ® @ @ • 


This  completes  the  data  for  the  first  reflective  surface.  The  next  surface 
to  be  read  in  will  be  the  TOP  PANEL,  so  identified  on  the  eighteenth  card. 
The  nineteenth  card  is  the  symmetry  card  for  this  panel.  This  panel  is 
not  symmetric  in  any  way  to  the  windshield;  therefore  SYMTRY  = 0.  for  this 
surface.  The  twentieth  card  contains  the  number  of  boundary  points  for 
this  panel  (RPTS  =4.).  The  twenty-first  and  twenty-second  cards  are  for 
the  four  boundary  points  (?)  (§)  (h)  (^  . This  completes  the  data  for 

the  second  reflective  surface.  The  next  surface  to  be  read  in  is  the  LEFT 
FWD  PANEL,  so  identified  on  the  twenty-third  card.  This  panel  is  symmetric 
to  the  right  forward  panel,  but  since  it  is  read  in  first,  SYMTRY  = 0 
for  this  surface.  The  twenty-fourth  card  is  for  the  number  of  boundary 
points  recorded  for  the  left-forward  surface;  this  number  is  three.  The 
twenty-fifth  and  twenty-sixth  cards  are  used  for  the  three  boundary  points 
for  this  surface,  (a)  . The  next  panel  to  be  read  in  should  be 

the  RIGHT  FWD  PANEL,  since  it  is  symmetric  to  the  panel  just  read  in,  and 
therefore  only  an  identification  card  and  symmetry  card  are  required.  The 
twenty-seventh  card  will  contain  the  identification  for  the  right-forward 
panel,  while  the  twenty-eighth  card  is  the  symmetry  card  which  has  a value 
of  one  since  the  right-and  left-forward  surfaces  are  symmetric  about  the 
x-z  plane.  The  next  panel  to  be  read  in  is  the  LEFT  AFT  PANEL,  so  identi- 
fied on  the  twenty-ninth  card.  Since  this  panel  is  not  symmetric  to  the 
right- forward  panel  just  read  in,  the  symmetry  card  (the  thirtieth  card) 
contains  a zero.  The  thirty-first  card  represents  the  number  of  boundary 
points  for  the  left-aft  panel  (RPTS  = 4.).  The  thirty-second  and  thirty- 
third  Ccurds  contain  the  four  boundary  points  for  this  surface,  (g)  @ 

. The  next  surface  to  be  read  in  is  the  RIGHT  AFT  PANEL,  so  identified 
on  the  thirty-fourth  card.  Since  this  surface  is  symmetric  to  the  surface 
just  read  in,  namely  the  left-aft  panel,  the  thirty-fifth  card  (the  sym- 
metry card)  will  have  a value  of  one,  representing  symmetry  edjout  the  x-z 
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1 


plane.  This  completes  the  input  data  for  this  case.  The  complete  input 
in  card  image  form  is  shown  in  Figure  A-2.  The  complete  results  from  this 
input  are  shown  on  pages  71  through  116. 
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Figure  A— 2.  Input  for  Case  I:  Card  Images. 
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OBSERVER  AZIMUTH 


ELEVATION  ANCLE 
PITCH  ATT  = 0.0 
ROLL  ATT  =0.0 


HEAD I NO  OF  AIRCRAFT  (DEC) 


52-1- 

Ijj  '-180 


^•0 

fo 

.00  -120.00  -60.00  0.00  80.00  120.00  U 

OBSERVER  AZIMUTH 

SUN  OLINT  SIGNATURE 

ELEVATION  ANCLE  = 0. 

PITCH  ATT  = 0.0 
ROLL  ATT  =0.0 


PROBABILITY0.004 


180.00 


PR0BAB1LITY0.064 


HEAOINC  OF  AIRCRAFT  (OECi 


uj  '-'leo.oo 

_/ 

UJ 


■120.00 


-60.00  0.00  60.00 
OBSERVER  AZIMUTH 


120.00 


180  pOw 
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The  second  case  represents  a configuration  with  curved  surfaces,  shown  in 
Figure  A- 3.  The  tctbulation  of  points  is  shown  in  Table  A-2.  Some  of  the 
output  data  has  been  suppressed  by  setting  PLOT  2.,  no  CALCOMP  plots. 

Also,  only  one  sun  elevation,  zero  degrees,  was  chosen  by  setting  GAMN  = 1. 
and  GAMI  » 0.  And  since  the  glints  from  the  windshield  alone  will  vary 
from  -180  to  180  degrees  in  azimuth  at  an  aircraft  heading  of  zero,  only 
two  headings  were  chosen  to  avoid  much  overlapping.  The  headings  chosen 
were  -45  and  45  degrees,  by  setting  PSIN  = 2.,  PSII  = -45.,  and  DPSI  = 90. 
The  card  images  for  this  case  are  shown  in  Figure  A-4.  The  default  option 
was  not  chosen  for  this  case,  in  order  to  demonstrate  the  additional  inputs 
required.  These  additional  inputs,  which  specify  which  terms  to  use  in 
the  curve-fit  analysis,  come  after  the  symmetry  card  for  each  reflective 
surface  input  except  those  that  are  symmetrical  to  the  preceding  surface. 
Additional  engineering  judgment  must  be  exercised  when  using  this  option. 
For  instemce,  the  first  and  second  surfaces  appear  to  be  a function  of  x, 
z»  y , and  z^.  Therefore,  a value  of  one  is  placed  on  the  data  Ceirds 
corresponding  to  these  inputs.  The  side  panels  seem  to  be  independent  of 
X and  X . Therefore,  a value  of  one  is  placed  or  the  data  cards  corres- 
ponding  to  z,  y^,  and  z locations  only.  The  results  from  the  inputs  of 
Figure  A-4  are  shown  on  pages  123  through  138.  It  was  foiind  using  these 
terms  that  the  error  for  the  curve-fit  of  the  windshield  was  substantially 
higher  than  that  for  the  other  surfaces.  This  brings  up  a point.  By  in- 
specting the  results  of  the  curve-fit  analysis,  the  program  user  can  decide 
if  one  or  more  of  the  surfaces  needs  to  be  broken  up  into  smaller  pieces. 

A closer  look  at  Figure  A- 3 reveals  that  part  of  the  windshield  actually 

2 

looks  like  a side  panel,  thus  independent  of  x and  x . If  cm  imaginary 
line  is  drawn  on  the  windshield  connecting  points  (5)  , two 

surfaces  are  formed.  Surface  number  one,  which  will  be  identified  as 
WINDSHIELD  FORWARD  SECTION  LEFT  SIDE,  will  be  formed  from  points  @ (f) 

© © ® ® @ @ . The  second  surface,  WINDSHIELD  AFT 

SECTION  LEFT  SIDE,  will  be  formed  from  points  ® ® © ® 

. The  modified  inputs  for  case  II  with  this  decomposition  of  the 
windshield  are  shown  in  Figure  A-5.  The  curve-fit  result  for  the  windshield 
as  shown  on  page  145  is  greatly  improved;  the  complete  results  are  on  pages 
139  through  158. 
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TABLE  A- 2.  COORDINATES  OF  CASE  II 

HELICOPTER 

CCMJFIGURATION 

Coordinates  of 

Center  of  Rotation 

(CG) 

Station 

Buttline 

Waterline 

114. 

0 

18. 

Coordinates  of 

windshield  Boundary 

Points 

Station 

Buttline 

Waterline 

72. 

0 

72. 

66. 

-33.75 

60. 

© 

60. 

-41.57 

48. 

© 

48. 

-48. 

24. 

© 

16.25 

-36. 

24. 

6.43 

-24. 

24. 

© 

0 

0 

24. 

6.43 

0 

48. 

® 

16.25 

0 

60. 

48. 

0 

72. 

48. 

-37.47 

54. 

Coordinates  of 

Lower  Panel  (Left  Side)  Boundary  Points 

Station 

Buttline 

Waterline 

© 

48. 

-41.57 

0 

30. 

-37.47 

0 

8.2 

-12. 

0 

© 

18.6 

-12. 

-12. 

48. 

-12. 

-22.48 

© 

48. 

-37.47 

- 6. 

Coordinates  of 

Forward  Side  Panel 

(Left  Side) 

Boundary  Points 

Station 

Buttline 

Waterline 

® 

72 

-48. 

24. 

© 

96. 

-48. 

24. 

0 

96. 

-44.5 

42. 

© 

96. 

-31.75 

60. 

© 

84. 

-31.75 

60. 

© 

72. 

-46.48 

36. 

119 

> 


TABLE  A-2  - Continued 


s of  Aft  Side  Panel  (Left  Side)  Boundary  Points 
Station  Buttline  Waterline 


108. 

144. 

144. 

144. 

18. 

)8. 


-48. 

-48. 

-44.5 

-31.75 

-31.75 

-44.5 


24. 

24. 

42. 

60. 

60. 

42. 


I 


L2 


r 


CURVED  PANELS  WITH  NO  FENCES  SAMPLE  CASE  TWO 


i^m 

d. 

0. 

<0. 

1. 

0. 

0. 

2. 

-45. 

90. 

114. 

0. 

18. 

isee. 

0. 

8. 

2. 

1. 

LEFT  SIDE  OF  WINDSHIELD 

rO. 

- 

1. 

0. 

1. 

0. 

0. 

0. 

ADDITIONAL 

1. 

1. 

1. 

INPUTS 

0. 

0. 

0. 

0. 

0. 

0. 

WHEN 

0. 

0. 

0. 

0. 

DEFLT-0 

11. 

72. 

0. 

72. 

66. 

-31.75 

60. 

60. 

-41.57 

48. 

48. 

-48. 

24. 

16.25 

-36. 

24. 

6.43 

-24. 

24. 

0. 

0. 

24. 

6.43 

0. 

48. 

16.25 

0. 

60. 

48. 

0. 

72. 

48. 

-37.47 

54. 

RIGHT 

SIDE  OF  WINSHIELD 

‘lower 

FRONT  PflNEL 

LEn  SIDE 

r0. 

- 

1. 

0. 

1. 

0. 

0. 

0. 

ADDITIONAL 

1. 

1. 

1. 

INPUTS 

0. 

0. 

0. 

0. 

0. 

0. 

WHEN 

0. 

0. 

0. 

0. 

DEFLT-0 

^6. 

■ 

48. 

-41.57 

0. 

30. 

-37.47 

0. 

8.2 

-12. 

0. 

18.6 

-12. 

-12. 

48. 

-12. 

-22.48 

48. 

-37.47 

—6. 

LOWER 

FRONT  PfiNEL 

RIGHT  SIDE 

FORURRD 

SIDE  PfiNEL 

LEFT  SIDE 

r0. 

.1 

0. 

0. 

1. 

0. 

0. 

0. 

ADDITIONAL 

0. 

1. 

1. 

INPUTS 

0. 

0. 

0. 

0. 

0. 

0. 

WHEN 

0. 

0. 

0. 

0. 

DEFLT-0 

6. 

- 

72. 

-48. 

24. 

96. 

-48. 

24. 

96. 

-44.5 

42. 

96. 

-31.75 

60. 

84. 

-31.75 

60. 

72. 

-46.48 

36. 

FORWRRD 

SIDE  PfiNEL 

RIGHT  SIDE 

fin’  SIDE 

PfiNEL  LEFT 

SIDE 

r0. 

• 

0. 

0. 

1. 

0. 

0. 

0. 

ADDITIONAL 

0. 

1. 

1. 

INPUTS 

0. 

0. 

0. 

0. 

0. 

0. 

WHEN 

0. 

0. 

0. 

0. 

DEFLT-0 

6. 

- 

108. 

-48. 

24. 

144. 

-48. 

24. 

144. 

-44.5 

42. 

144. 

-31.75 

60. 

108. 

-31.75 

60. 

108. 

-44.5 

42. 

PFT  SIDE  PflrCL  RIGHT  SIDE 
1. 

Figure  A- 4.  Input  for  Case  II;  Card  Images. 
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CURVED  PANELS  WITH  NO  FENCES  SAMPLE  CASE  TWO 
- MODIFIED  WINDSHIELD 


0. 

C • 

0. 

0. 

1. 

0. 

0. 

Ca 

-45. 

114. 

0. 

13. 

1500. 

-10. 

10. 

-130. 

130. 

0. 

10. 

Ca 

1. 

WINDSHIELD  FORWfiRD  SECTION 

1% 

LEFT  SIDE 

1. 

0. 

1. 

0. 

0. 

1. 

1. 

1. 

l7i 

LT4 

0. 

U« 

0. 

Vm 

0. 

Ua 

0. 

Va 

43. 

— 4y . 

£4. 

16.25 

-36. 

6.43 

— y4. 

£4. 

0. 

0. 

6. 43 

0. 

4S. 

16.25 

0. 

43. 

0. 

1 L.a 

48. 

-37.47 

WINDSHIELD  FORWRRlD  SECTION 

RIGHT  SIDE 

1. 

WINDSHIELD  RFT  : 

SECTION  LEFT 

SIDE 

0. 

0. 

0. 

1. 

0. 

0. 

0. 

1. 

i; 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

6. 

?£. 

0. 

72. 

66. 

—31 . 7^ 

S0. 

-41.57 

48a 

48. 

-43. 

48. 

-37.47 

54. 

48. 

0. 

WINDSHIELD  AFT 

SECTION  RIGHT  SIDE 

1. 

LOWER 

FR:0NT  PFirCL  LEFT  SIDE 

0. 

1. 

0. 

1. 

0. 

0. 

1. 

1. 

1. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

6. 

43. 

-41.57 

0. 

30. 

-37.47 

3.2 

-12. 

0. 

18.6 

-12. 

48. 

-12. 

-22.48 

48. 

-37.  47 

LOWER  FRONT  PANEL  RIGHT  SIDE 

1. 


0 

0 


0. 


y. 


Figure  A-S.  Input  for  Case  II  (Windshield  Modified);  Card  Images  (Sheet  1 of  2). 
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Tj-  Tj-  o tT  • • O OJ 

OJ  CO  Vjji  lO  O O 'D  OJ  r-. 


FORWARD  SIDE  PANEL  LEFT  SIDE 


Q. 

0. 

1. 

0. 

0. 

0. 

0. 

1. 

1. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

6. 

1 

-48. 

c!4. 

96. 

-48. 

£4. 

'r'S. 

-44.5 

42. 

96. 

-31.75 

60. 

S!4. 

-31.75 

60. 

■'£. 

-46.48 

36. 

FCiRWARD 

SIDE  PArCL 

RIGHT 

SIDE 

1. 

AFT  SIDE 

PANEL  LEFT 

SIDE 

y . 

0. 

0. 

1. 

0. 

0. 

0. 

0. 

1. 

1. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

6. 

6. 

0. 

y . 

108. 

-48. 

£4. 

144. 

— 4i:L 

i:!4. 

144. 

-44.5 

4£. 

144. 

— 3i . 75 

60. 

108. 

-51.75 

60. 

108. 

-44.5 

4^. 

l=fT  SIDE  PANEL  RHilT  SIDE 

1. 


IT 

I 

E 

i- 

i 

I 

[ 

Figure  A— 5.  Input  for  Case  II  (Windshield  Modified);  Card  Images  (Sheet  2 of  2). 
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CUliVCO  PANCLS  VITH  NO  FCNCCS  SANPtt  CASC  TWO 


CUKVCD  PANELS  WITH  NO  FENCES  SAMPLE  CASE  TWO 


j ui  in  « 


w o o e I 

o • • • 

<j  cs  e o < 


w o z 

< ►MO 

> » K K O 


O ► • • 

O m r<t 

O 9 

I I 


& S o o o 

cc  at  • • • 

w uj  e o o 


► M tk  tfe  Ife 

M Ui  Ui  W 

o o e 
u u u 

* 0^ 


« e cj  o 

3 X • • K 

» e o 4k 

4a.  n ^ 


c z a: 

K OC  S ^ 

uj  w M in  « 


z o o 

4aJ  Z M 

C ••«  ► X 


uj  o o o • 

Ck  • • t • 

w o o o o 


^ o z 

« ► lA  O 

L ► > K K K U 


« ► M OC  C 

4a.  D K 4iJ  4Ai 

Z O.  Z ► ► 

S Z 4J 


CURVED  PRNCIS  UlIH  NO  FENCES  S4RPLE  CRSF  TWO 


CURVtD  PANELS  yiTN  NO  FENCES  SAHPLC  CASE  TUO 


REFLCCTTVC  SURFRCC  HO.  1 1 OCNTIFIC *TI ON:  LEFT  SIPE  OF  WINDSHIELD 


RCrtCCTtVC  SOKfACC  NO*  2 I DCNTIF ICAT ION:  RICHT  StOr  OF  WINSHltLO 


CUftVCO  PANELS  WITH  NO  FENCES  SA-PLE  CASE  TWf 
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CURVCO  PANCLS  WITH  NO  FENCES  SAMPLE  CASE  TWO 


cunyco  PUHtLS  with  ho  FCHces  s*hple  c»sc  two 


•acFLCCTivr  sunftcc  no,  s lorNTincATicN:  yiNcsHirLo  »ft  stction  left  sirt 


ctsiurr  Tturs  Ff/R  cuRve-» ittinc 

'll  » TFR*  1.  tfi  T TFRR  C.  «’>  7 TFRR  1.  Ml  »T  TF**  '•>  «Z  TrRR 

f‘)  TtRI*  a.  (7)  XSn  TER*'  1.  (F)  YSC  TER"  1.  «<»)  2Sft  TfR*  1. 


I 


Z •••  Z •••  Z •••  X ••• 

€<  o es  Vi  CAC.O  uco  v^  o 

C o o o 

* * » » 

a ifcifeife  a u.iku.  a a u.u.14 

UWU  '••U'lfc  WWW  Www 

••  oqc  m**  oco  000  ooo 

a wvu  c/  X oou  c.a  c.wu  ciz  woo 
o I a *0  10 
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HOOT  rttl  SOUtltr  rRROK  Of  fITi  a.9lt«9t-(i9 


POOT  rUN  SOU*Rf  ERROR  OF  FITi  (!.73?2r*(iR 
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1 
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APPENDIX  B 
PROGRAM  LISTING 


6 LEVCL  21 


r'AIN 


DATE  = 78103 


08/5A/4A 


570 


Co*******  PROGRAHMED  by  FRED  klHITE  ON  BOEING  IBH 

C***********  BOEING  VERTOL*  PHILADELPHIA,  PENNA.  ************ 

- * *-»««ji »J.*  AREA^  COr.E  215  B.22 -22'>8  -».o«*-««o«-««-«.«*-*« • *•  oj 

C 


RFAl  *H  lirilRFt«»f;l.FrURFt».Q1.5;lf. 
REAL*8  CFALPO) 

DIMENSION  ALP(5P) ,PALP»5C> 
DIMENSION  SlGFt.->11.SIGS(12> 


OOOOOOAl 
OQOOOOA2 
non  non  A 3 


DIMENSION  COEFFP* 80 > ,CCEFSE < 120 > 

DIMENSION  OATAC(I0),DATAFe2C0),CATAS(1500),TlTLES(96O> 

dimension  ALPHA2( 80C ),8ETA2(e  DC) * INTRFL(F00>,ISHAOU(80 
1 INTERFCROO) 

Ill«ASLDN,^L.tQU.l32J. 


OOOOOCAA 
C00000A5 
jmo  aium. 
00000080 
00000070 
00000080 
oonoooBO 
OOOOCICO 
nnn  n 01 1 p 


0)  < 


0000  0120 
000001 30 
00000140 


DIMENSION  IBUF<1  000  )*XPLOTl63>,YPLOT<63),APLCT(  432),ePLOT<  432) 
COMMCN/70LRS/TOLFGN,TOLSNT,TOLQPT,TOHNT.TOLCK1,TOLCK2, 

L-  I01.CK4,DEL.T*TDL£lR*ANGSI2.TaL£LJ 

T0LEQN=0. 000001 
TOLSNTrO.OOOl 

1 0L£PX=O.fl  aOG  ORCL 

TCL1NT=0.02 
TOLCK1=0.02 

TOLCm2=0.02..-  - 

70LC!<4r';.02 
OEL7=0.ia 

TOLClRiG.aoaL  

ANGST2=9C . 

TOLPLT=0. 0001 

OD-5>-l—L,.30 

5 OAYACimO.O 
DO  8 1=1,200 

lAF-LXlsO.. 

DO  7 1=1,1500 
7 OATAS«I>=0. 

CAU INPUT 

KPRINT=DATAC12)*0.5 
lPLOT=DATAC12‘t)*0.5 
LP-RQBU=JJA  TAC.125  ).♦  0*  5_ 


000  0 0150 
0000  0160 

aODOJllZD- 

00000171 
00000172 
00000173 
000  0 01  74 
00000175 
._.?.Q£,QC17_6 
00000177 
000  0 0178 
-JU1J1IU1L79 
00000180 
00000181 
^00000182 
OOOOOIBO 
00000200 
oonn'-?!  0 


CALL  0UTPT1CDATAC,CATAF,0ATAS,TITLES,MF,MS) 
CALL  INPUT(OATAC,OATAF .DATAS,TITLES,MF ,MS, 1 ) 

4*FE*.Ct=OATAt«22>*L,E- — 

I FINFENCE.EQ.O)  GO  TO  14 
CO  12  I=1,NFENCE 

-LsilF  t I-A 


12 


NPTS=DATAF(L)*0.5 

ISTART=L*1 

CAU CfCTFCOAT  Af_,ICIAlU  ,NeTS..CQEE£a *1-,S1G.XCUR-F ,f  CURE  1 

SIGF( I )=SIG 
CONTINUE 

-LtdlPIUNI  .Li*M.C  Q_tC_i3 

CALL  0UTPT2(C0EFFB,SIGF,NFENCE,TITLES,1  > 


00000220 

00000230 

!LaaiQ2Aa, 

00000250 

00000260 

OOiUlZJ.O 

000  0 0280 
00''  0 02°0 

aoo-ooiciL 

oa"oo3i 0 
00000320 

aocaaasa 

0CC0C34C 

00000350 

oor.acstc. 

00000370 
0000  0380 
000003°0 


00000400 


160 


i 


G LfVEL  21 


MAIN 


CATE 


TEICJ 


OP/b»/AA 


lA  CONTINUE 

NPANEL=DATAC<23)«0.5 

CO.,  li  .1=1  iNEAilU. 

L-ME( 1) 

NPT,S  = OATAS(L)  ♦O.b 

XSTARIiL.‘l — 

CALL  CFITS(DATAC.0ATAS.1START«NPTS«C0EFSB. l»SIGtXCUPF»FCURF) 
SI6S«1)=SIC 

-15 — COMTlItUf 


COOOOAIP 

00000*20 

00000*30 


00000**0 
0000  0*50 

00CCCA70 

00000*P0 

00000490 


15 


IF<KPRINT.LT.2)  go  to  16 

CALL  OUTPTZlCOCFsBtS IGS»NPANEL*TI TLES,2) 
CQNtlNUE 


00000500 

0000051  0 
n nn  n 


NGAM=CATAC<6)»0.b 
DG  *!*  = D A TA  C ( B) 
GAMMArOATACCTl-DGAP 


00000550 

000005*0 


IF  ( IPLOT  .EQ.O.OR  . IPLCT.  EC  .2)  CO  TO  17 
JC0NT=C 

f LllI  GllHUF.iaOD.H 


17  CONTINUE 

00  30  I=ltNGAM 
GAMMA  =GAf1HAffljiAt!.. 


00000560 

00000570 

OOOOOSBQ 

00000590 

00000600 


CALL  &LINT(OATAr,DATAS.MS.COEFSE,GAMMA»ALPHA2«eETA2.INTRFL» 

1 I$PADW.INTFRFfOwOBKF,DATAF,MF) 

IF  tKrRINI.LI.ll  GQ  TU-  la_ 

CALL  0UTPTitDATAt,DATAS.MS.TITLES,GAMMA,ALPHA2.6ETA2,INTRFL» 

1 ISHACu.INTERF) 

. IB-  .CUNTINUB  - — 

call  MINMAX<0ATACfCATAS.MS»ALPHA2,PETA2,APL0T,BPL07»INTRFL. 

1 ISHADW.INTERF.IPLTOL) 

IF  llF(lQaL.ta..Gl  GD— TC...15 


00000620 
0000  06  30 
.aiG.ac6Ao. 
00000650 
00000660 
-!m.ac6.7,a 


OOQ006SO 
00000690 
nnonnTPn 

CALL  PROPL(DATAC.CAT*S,MS.APLOTfHPLOTfXPLOT,YPLOT,DUORKS.OWKSP2,  00000710 
1 VPHOBL t IPLTOL > 00000720 

-la coal  1 dui goaocTio 


20 


IF ( IPLOT.EC.O .OR . IPLOT.EQ. 1 > GO  TO  20 

CALL  PLOT OL (OATAC.OATAS. MS. TITLES. gamma. IP  LOT, B PLOT, 

1..KPSQEL.1PL.LOL) 

CONTINUE 

IF(IPlOT.EO.O.OP.1PLOT.EG.2>  GO  TO  21 
XAU. 


000007*0 
OOC  C 0750 
nor.ac76a 


1 ISHADW,INTERF,XPL0T,YPL0T,IPUF,JC0NT.0W0RKS,DWKSP2,VPR0PL) 
21  CONTINUE 

.za crntTutur 


00000770 

0000C7P0 

-gjacozs-g 

OOCO'SCC 
0000081  0 
onoonppn 


-U- 


IF(IPLOT.EQ.C.OR.IPLOT.E0.2)  go  to  31 
CALL  FLOT(10.0,0.999) 

rONTl NUF 


0CC0CB30 
0 00  0 06*  0 
nnnr  ormr 


STOP 

end 


00000860 

OOOOOfiTO 


G LEVCL  21 


INPUT 


0*TE 


78102 


10/47/48 


SUBROUTINE  INPUT«0*T*Ct0«T*F«0*T4StTITLEStNF«NS<I0TR«N> 

01  HENS  I ON  0*T*CC1>,0*T*F(1),0«TAS(1>«TITLES(I><HF<1>«PS(1> 
data  BLANK/IH  / 

00000880 

00000890 

00000900 

IF ( IOTRAN.EQ.1 > GO  TO  80 

00000930 

REA0(S,101)  ( TITLES!  I > «I  = 1,72) 

00000940 

00  2 1=73,80 

00000950 

2 

TITL£S<I>=BLANK 

00000960 

LT=80 

00000970 

REAOI5.100I  OEFlT. PRINT. theta. PHI  . 

00001000 

1FIDEFLT.LT. 0.5)  OEFLT=0. 

0000 1010 

IF!PrINT.LT.0.5)  PRINT=0. 

0000  1020 

IF(APS(THETA).LT. 0.001)  THCTA=0. 

00001030 

IF(ABS(PHI>.LT.0.0Q1)  PHI=0. 

00001040 

OATAC! 1)=0EFLT 

00001050 

0ATAC(2)=PRINT 

0000106  0 

OATAC  <3  > = 0 . 

00001070 

OATAC  «4>  = THETA 

OOOOIOBO 

0ATAC(5>=PHI 

00001090 

REAOtS.lOO)  GANN. GAHI.OGAH.PS IN. PS  II. OP  SI 

00001100 

ifigahn.gt.o.s)  go  To  3 

00001110 

GAHN=3. 

GAHI=0. 

OG001130 

□ GAN=3a . 

00001140 

J 

IFCPSIN.GT.0.5)  GO  TO  5 

00001150 

PSIN=9. 

00001160 

PSII=-lbO. 

D0001I70 

0PSIs<»5. 

00001180 

5 

OATAC (8) =GAHN 

00001190 

DATAC(7)  = G4.»I 

00001200 

0ATAC<R)=0GAM 

oooofzic 

0ATAC<9).=PSIN 

00001220 

OATAC(10)=PSI1 

00001230 

0ATAC(11)=0PSI 

00001240 

READCS.IOO)  XROT.YROT, 2R0T 

00001250 

I F CAPS(XROT)  .LT  .0  . Cfll  > XROTrC. 

00001260 

IF(ABS<YROT).LT.O.OCI)  yrot=o. 

00001270 

IF(ABS(ZROT)  .LT. 0.001)  ZROT=0. 

00001280 

OATAC<12)=XROT 

00001290 

OATAC<13)=YROT 

0ATAC«14)=ZR0T 

00001310 

READ! 5, 100)  0 1 S TC , ELv« IN .ELVHA X , AZMMI N, AZN M AX 

00001320 

IF (AeS(OISTG).LT. 0.001)  OISTG=1000. 

00001330 

IF(ABS(ELVNIN).GT.r.001  ) GO  TO  6 

00C01340 

IF(ABS<ELVHAX).GT. 0.001)  GO  TO  6 

00001350 

ELVMIN=-10. 

00001360 

ELVHAX=10. 

oono 1370 

6 

1F(ABS(AZHHIN).GT.0.001>  GO  TO  7 

00001380 

IF<ABS(AZPHAX).GT. 0.001)  GO  TO  7 

00C013'’0 

AZHH1N  = -180. 

00001400 

AZHNAXSIBO. 

00001410 

7 

OATAC! 15)=DISTG 

vnTRTn 

OATAC! 16)=ELVniN 

00001430 

162 


G LEVEL 


INPUT 


DATE  = 78102 


10/47/ A8 


21 


0ATAC(17)=ELVNAX 

00001440 

OATACC 18)=AZMNIN 

00001450 

REAOIStlOO)  FENCES,PANELStPLOT*PROBL 

00001470 

DATAC  «22) =FCNCES 

00001480 

QATACt231=PANFLS 

^1— ill  ill I I 

DATAC<24JrPL0T 

00001500 

0ATAC<25)=PR0BL 

00001510 

LOCsl 

—^—11  null 

IF(DATAC(22>.LT.0.5)  GO  TO  30 

00001530 

NrENCE=0ATAC(22)-»0.5 

00001540 

00  20  I=1«NFENCE 

00001550 

REAO(5ilOO)  OATAF(LOC) 

00001560 

NPTS=DaT4F(LOC)*0.5 

00001570 

■(FCtlsLQC 

N3=3«NPTS 

0C001S90 

READ* 5t 100 > lOATAF  <L0C*J) »J  = 1 ,N3> 

oooouoo 

2Q 

LQCiLOC^Nl-vl 

30 

NPANEL  = DA  TAC<23>»0.5 

00001620 

LOC  = l 

OOOO 1630 

00  40  1=1.NPANEL 

00001640 

NSd  ) =LOC»20 

00001650 

RE4D(5*101)  (TITLES(LT*J),J=1,72) 

00CC1660 

00  301  J = 73.80 

00001670 

301 

T ITLES(LT*J)=BLANK 

0000  1680 

LTrLT*80 

00001690 

REAO(5.100)  0ATA5<L0C> 

00001700 

IF<OATAS(LOC).CT.0.5»  GO  TO  32 

00001710 

IF(DATAC(  D.GT.O.S)  GO  TO  310 

00001720 

00001730 

RE  APIS?  100)  <OAT4S(LOC»J)»J=10.19> 

0000  1740 

310 

REAO(5tlOG)  0ATAS(L0C*20) 

0000 1790 

nr^2ni  ♦n.s 

finfiflififin 

N3P20=20»3*SPTS 

00001810 

READCStlOO)  <OATASCLOC*J>«vl=21»N3P20) 

00001820 

60  TO  38 

00001830 

32 

IFU.EQ.D  GO  TO  70 

00001840 

ISYMrOilTAS<LOC)^0.5 

00001850 

00  34  J=1.20 

00001860 

34 

DATAS( LOC* J)=0aT4  S(LOCP»U) 

00001870 

NPTS=DATAS(LOC*20)*0.5 

00001880 

DO  35  J=1.NPTS 

^■^—1  III!  Ill 

L={J-1 ) *3 

00001900 

OATAS (L0C*21*L)=DATASCL0CP*21 ♦L) 

00001910 

DATAS(L0C»22*L):DATASIL0CP»22*L) 

00001920 

DATAS  <L0C*23*L)=0 ATAS(L0CP*23*L) 

00001930 

IFdSYN.EQ.l ) OATAS(LOC*22*L)  r-0  AT  AS  ( LO  C»2  2*L  ) 

00001900 

1F11ST8.EQ.2)  DATAS«L0C*23*L>  = rDATAS(L0C»2  3*L) 

0000195  0 

IF(ISTM.E0.3)  0ATAS<L0C*21»L)  =-D  A TAS  ( LOC  *2  1 »L  > 

00001960 

35 

CONTINUE 

00001970 

163 


— ^ 


G LEVEL 

21  INPUT  date  = 78102 

10/47/88 

NT  = NPANEL»1 

00002000 

DO  68  1=1, NT 

00002010 

L=«I-1)*8Q 

00002020 

NB  = 0 

00002030 

00  62  Jsl«80 

0000  208  0 

1F( TI TLES (L*d) .NE .PLANK ) GO  TO  63 

00002050 

G2 

Ne=NB*l 

00002060 

63 

1F«NB.EG.80>  GO  TO  68 

(70002070 

IFINB.EQ.Q)  GO  TO  6h 

——mill 

DO  66  J=1.80 

00002090 

K=j*Ne 

00002100 

IFIK.GT.BO)  GO  TO  6A 

OOQ0211C 

TITLES<L*J)  = TITLFS(L*K) 

00002120 

GO  TO  66 

00002130 

64 

TITLES(L«d>=eLANK 

00002180 

66 

CONTINUE 

'.'002150 

68 

CONTI NUE 

000  0 2160 

RCTURN 

00002170 

70 

URITC(6,500) 

00002180 

STOP 

00002190 

500 

FORMAT!! HI , ’ERROR  IN  INPUT,  STMMETRICAL  PANEL  CANNOT 

BE  FIRST*)  00002200 

lOO 

FORMA  T(faE ) 0. 0 ) 

00002210 

ICl 

FORr*AT  (80  A1  > 

00002220 

80 

CONT INUE 

00002230 

IF«DATAC<22>.LT.n.5)  GO  TO  90 

00002280 

NFENCE=DATAC(22)*0.5 

00002250 

no  P6  I=1,NFENCL 

00002260 

LCC=MF(I) 

000  02270 

NP  TSrCAT  AF  ( LOG  }*C  .5 

000  02280 

DO  86  J=1,NPTS 

00002290 

00  8b  Jd=l,3 

00002300 

LJ=< J“1 )»3^JJ 

00002310 

66 

OATAF(LOC*Ld)=DATAF (L OC ’L J ) -0  A T AC  < J 1 1 ) 

00002320 

S8 

continue 

00002330 

“iO 

CONTINUE 

00002380 

NPANEL=OATACC23)*0.5 

00002350 

DO  96  t=l,NPANEL 

00002360 

LOCF20=MS(I) 

00002370 

NPTS=DATAS(LOCP20)«0 .5 

00CG2380 

DO  9»  J = 1,NPTS 

00002390 

DC  98  JJ=1,3 

000024CO 

LJ=(  j-n*3**jvj 

00002801 

94 

DATAS(LOCP20*LJ)=OATAS(LOCP2C,Ld)-OATAC(Jd»ll> 

00002810 

96 

CONTI NUE 

00002820 

RETURN 

00002830 

END 

00002440 

L6 


G LEVEL  2 


OUTPTl 


0*TE  = 78102 


io/*7/4e 


2 IH  ,5X,*PL0T  OPTION* , lOX tF2. 0, ex, •PFOBABILITY  0PT*,6X,  000025^0 

3 F2.0)  00002550 

yPITE  t ^ Oe  ) D A2  * C < AJ  I 0 ATAC(5).CATAC(  10) 00002560 

IG*  FORMAT( 1H0,2X,»**AIRCRAFT  INITIAL  EULER  OR IENTATION »/ IH  » 5X , "P I TCMOOOO 2570 
l*,13X,F6.1,7X,*R0LL*,14X,F6.1,7x,*TAy*,l5x,F6.1)  00002580 

yRITEtb,105>  IDATACII ) .1=9,11) 00002590 

105  F ORMA T ( IHO ,2X , •••selected  TAU  ANGLE  fl OT ATI ONS •/ 1 H ,5X,*NO.  OF  ANGL00002600 
2ES*,6X,F4.0,ex,*IMTIAL  ANGLE  *,  5X  ,F6.  1, 7X  ,*  ANGLE  INCRE  PE  NT  • ,3X  .F6 . 0000261  0 

3 IJ 00002620 

URITE(6,f06)  (OATACt  I)  ,1  = 6,8)  00002630 


106  F0RMAT(1H0,2X, •••selected  SUN  ELEVATIONS*/ IH  ,5X,*N0.  OF  ANGLES*,  00002640 

1 6X.F4.0,8X ,*INIT I AL  ANGLE  * ,5X .F6 .1 .7X . *ANGLC  1 NCR EMENT* .3X , F6 .1 ) 00002650 

URITE«6,107)  OATACt 15>,0ATAC( 12 > ,OATAC( 13) ,OATaC ( 14)  00002660 

107  FORHAT(1HO,2X, •••REFERENCE  INFORPATION  FOR  MEASURING  SUN  GLINT  ANG00002670 

1LES*/1H  ,5X,*RE F£R  ENCE  POINT  FOR  MEA SPRING  ANGLES  OF  GLINT  S I GN ATUQ 00 02680 
2RE*,5X,*0ISTANCE  AT  WHICH  GLINT  SIGNATURE  FORMED  IF  T . ) *,  F 1 0 .3/ IH  ,00002690 
3 1HX,*X*,9X,*Y*,9X,*Z*/1H  , 14X,3< F7.2, 3X> ) OOOC2700 

WRITE!  6,108)_  DATAC(18>,DA  T Ai  ( J 6 ) ,0  A_T  A Q ( 1 97  ,_0  AT  A t<  J 7_) 00002710 

ira  FORMAT! IHO, 2X, •••PLOT  SCALING*/iH  ,5X,*GLINT  A2 I MUTH * , 18X , *6L INT  E00002720 
1LEVaTION*/1H  , 7X, •“ in IMUM* ,8X ,F6, 1 , 1 ox, *M I NIMUM * ,8X,F6 .1 /I H ,7x,  00002730 

2 *MAX1MUM  * ,6X  , F 6. 1,10  X , • MA  X 1“  UM_*_J  8 X^F  6^,  17 0CS02740 

WRITE  16,10  9)  DA T AC! 16), DA T AC! 16), DAT  AC! 19), DA TAC!  17)  00002750 

109  FORMAT! IHO, 2X , •••eouNDAR lES  FOR  CALCULATING  FR 06 AB I L IT Y* / IH  ,5X,  00002760 

1 *GLINT  AZIMUTH*. I8X .*GLINT  ELF V AT  I ON */ IH  . 7X . * MIN  I MUM* , 8X . F6 . 1 . OOC0277G 

2 1 OX , *MINIMUM* ,8X ,Fb.l/ IH  , 7X , *M AX  I MUM  * ,8 X , F6 . 1 , 1 0 X , * MAX  I MUM  * , 8X , 00002780 

3F6.1)  00002790 

IF  !OA TAC! 22) .LT.O .5)  GO  TO  25  00002ROC 


WRITE!6,102)  ! T IT lES ! MMM) ,MMM=1 , 80 ) 00002810 

WRITE(6,101)  00002820 

WRlTEl6,110)  DATAC(22) 00002830 

F0RMATUH0.2X, •••FENCE  0ATA*/1M  ,5X,*NO.  OF  F ENCE S * , 7X  ,F 3 • 0 ) 0000284  0 

LINES=5  00002850 

NFENCE=D  A T_A  C!22)*0.S 0000286 1 

DO  20  I=1,NFENCE  00002870 

M2=MF!I)  00002880 

MRJIE-!4iLll> 016.0^890 

L1NES=LINES*2  00002900 

IF!  !L1NES*6).LE.50)  GO  TO  4 0000291  0 

WRlTE!b.l  02)  !T  1TLES!MMM)  tMMMil  ,80) 30CQ2920 

00002930 
00002940 
00002950 


WRITE!6,10t> 
WRITC!6,111) 
NCS  = 5 


WRITE!6,112)  I,0ATAF!M2) 


00002960 


G LCVEL  21 


DATE  = 78102 


10/47/A8 


111  FORHATdHO)  00002970 

112  FORHAKIH  t5X»*«FENCE  NO.'flS/lH  ,7X,»NUNaER  OF  POINTS* , 4X  «F3 . 0/  00002980 

NOART  POINTS  — BOOT  AXES  (STATION! INE.BUTTL INE. RATE  ROOD  0299 


XLINEl*  00003000 

1 /IH  *3<  18X,«X»,sX.*T».9X,  00003010 


LINES=LINES*8 

NPTS=OATAF(M2)*0.5 


iinmiMO 


00003030 

00003040 


00  10  L=ltND3 
Jr(L-l)  •9««2 


00003060 
0000  3070 


KP1=R*1 

KP2=X»2 

IF (CLINES*!) .LE. 50)  GO  TO  8 

WRITE(6,102)  (TITLESCMPfl)  ,PMH  = 1 ,8  0) 
RRITE(6,101) 

URITECb,! 13)  I 


OOC03O9O 

00003100 

00003110 

0C0O3120 

00003130 

00003140 


113  FORHAT(1HO,5X, ••FENCE  N0.»,I3/1H  , 7X ,• BOUND AR T POINTS  — BOOT  AXES00003150 
X (STATIONLINE, BUTTLINE, WATERLINE)  • 0000  3160 

I /IH  .31 1 RX. 'X* .9X . •¥ • .9X > '2*«2X ) ) 


LINES=7  00003180 

8 WRITE(6,114)  K,(DATAF  ( J*M)  ,N:1 ,3)  ,KP1,  < OAT  AFC  J*M)  ,M:4,6)  , 00003190 

1 KP2,(DATAF(J*H),m  = 7.9) 00003200 

114  FORHATdH  ,2X,3(6X,»(»,I2,»)»*2X,2(F7.2,3X),F7.2,2X))  00003210 

10  LINES=LINES*1  00003220 

IF(3.NQ3.EQ.NPTS)  GQ_j;Q_2a 0C00323Q 

J=ND3*9*N2  0000  3240 

K=3*N03*!  0000325G 

1F((3«ND3«1)  B T SI  _W  R LT  L(  6 , UJSl^AU  Q AT  AF  ( J U aOnaSELC. 

KP1=K*1  00003270 

IF((3*N03*2).EQ,NPTS)  WRITE(6,116)  K , ( OA TA F ( J*B ) , n=l , 3) , 00003280 

( 


LINES  = LINES*1  0000  3300 

115  FORMATCIM  ,«X  , • ( • , 1 2 ,•  ) • , 2 X ,2  ( F 7. 2 , 3X  ) ,F  7.  2 ) 00003310 

116  FQRMAKIH  . 2X  . 2 ( 6X  . » ( « . 1 2 . • 17  , 2Lf  .7  .2  t IXJ  ,fJ..7Li2XU 0000-332.0- 

20  CONTINUE  00003330 

25  CONTINUE  00003340 

WRITE(6.1Q21  (TITLES  (MHH)  .MHH=1.B01 00003350 

URITECb, 101)  00003360 

NPANEL=OATAC(23)*0.5  00003370 

tfiJJiI6dJJl_DA.TAC(23)_  DiaOJlBL 

117  FOBMAT(1HO,2X, •••REFLECTIVE  SURFACE  OATA^/IH  ,5X,^N0.  OF  SURF  ACE S • 000 0 3390 

I,5X,F3.0)  00003400 

LlNESsb gC00341C 

LTsO  OCC03420 

00  70  I=1,NPANEL  00003430 

LTsLT  *80 0000  34 4^ 

H2=PS(I)-20  0000  3450 

IF((LINES*6).LE.50)  GO  TO  27  00003460 

WRITE  (6.1 0 21  CTITLESCHWN)  .WHf«  = l,a01 attIlQ38JJl. 

URITE(6,101)  000034RO 


G LEVCL  21 


OUTPTl 


DATE 


76102 


10/A7/A8 


LINCS=3 

27  URITE«G<118I  I • I T I TLCS < LT« J) • J=1 *80 > 


1 BOAl) 

LlNeS=LINeS*2 

IF<ILINES*S)  .LE.5  0 )_  GO  T 0 30 


00003490 
00003500 
00003510 
00003520 
0000  3530 
0000  3540 


WRITE (6,1 02)  <TITLES(BRH>,«HN=1,80) 
WRITE(6,101) 

LINES=3 


00003550 

00003560 

00003570 


WRITE  (6, 118)  I,(TITLES(LT*J),J:l,8  0) 
LINES=LINES,2 

3 0 WR  1 TE  ( 6il  2 0 ) ( DA  T AS  ( H2*J  )j  -JM j(9_> 


00003580 
0000  3590 
00003600 


120  FORHAT(lhO,7X, ‘CURVE-FIT  INPUT  DATA’/IH  ,9X, "DESIRED  TERPS  FOR  CUR00CC3610 
IVE-FI TT ING*/1M  *11X."<1)  X TERM  • ,F4 . 0, 5X , • ( 2 ) Y TERM  "tFA.O,  0OCO362O 

2 5X*"(3 1 ^ I ER  ^ I*  1 -HI 4 . g . 5X-i-»  ( 51_XZ TERM  00003630 

3",F4.0/1H  ,11X,"(6)  YZ  TERM  • ,F  <t  . 0 , 5X  , ' ( 7 ) XSQ  TERM  "tFA.O.SX,  OG00  364  Q 
4 "(8)  YSO  TERM  ",64.0 ,5X, • (9)  2SQ  TERM  •,F4.0>  00003650 

i -LINES  "5 00003660 


IF((LINCS"6>.LE.50)  GO  TO  40 
WRITE(6,102)  (TITLES! MMM),MMM =1,80) 
_w_RXTE,iL6jL!LU 


00003670 

00003680 

opaa.369a 


00003700 
0000  3710 
00003720 


LINES=3 

WR  ITE (6,118)  I,{TI TLES(LT»J), J=1 ,80) 

L INESrLlNE S ♦ 2 

40  WRITE(6,121)  (DATAS(M2,J),J=10,1B)  00003730 

121  FORMAT(1HO,9X, "INPUT  VALUES  FOR  COEFFICIENTS  OF  SURFACE"/1H  ,11X,  00003740 

1 "(II  X COEF  ^jj:i2.4,5X,  "(2)  Y COEF  2_"El?."4"5Xr_^3)  Z COEF  "00003750 

X,E12.*/1H  ,11X,  00003760 

2 "(4)  XY  COEF  ",E12.4,5X,"(5)  XZ  COEF  " , E 1 2 . 4 , 5X,  " ( 6)  YZ  COEF  "00C03770 
3.E12.4/1H  ,I1X,"(7)  XSQ  COEF  " ,E  1 2. 4 , 5X  « " ( 8 > YSQ  COEF  ",E12.4,5X,  0000  3760 
4"  "(9)  YsQ  COEF  '",'e12.4)'  ’ - " 00003790 

WRITE(6,123)  0ATAS(M2*19)  00003800 

123  FORMATdH  , I Q X.,  MJ_0  )_.C  MSI*  (IT . 4J 00003810 


LINES=LINES"6 

IF((L1NES,4).LE.50)  GO  TO  50 

WR  ITE  (6.102  )_.XTiTLES.(I!MM)  ,MMM  =!,£  OJ 

WRITC(6,101) 

LINES=3 

JERllC  ttlii .«  _t it l tULI-rJ.)  t J.=L..e  IL 


LINES=LINCS"2 
50  WRITE(6,122) 

122  FORMATdH  0, 7X  ,_"1CUNDAR  Y POINTS  -- 


BODY  AXES 


00003850 
00003860 
00003870 
0000  3880 
00003890 

(STATIONLINE.BUTTLINEOOCO  390  0 


X, WATERLINE) "/IH  , 3 ( 18 X , "X • , 9X , 
1 "Y",9X,"Z",2X)) 


00003910 

00003920 


NPTS=DATAS(MJ,20)"0.5 

00003940 

ND3=NPTS/3 

0 00  0 3950 

DO  60  L=1.N03 

0000  3960 

J=(L-1),9*20*M2 

00003970 

8=(L-1)*3*1 

00003980 

XP1=I(*1 

KP2=K»2 

00004000 

6 level 


OUTPTl 


OATC 


78102 


10/47/*8 


IF((LINES*1>.LE«50)  go  to  52 
URlTE(b»102)  <TITLES«H«M>  *MHH=lt8  0) 
URITE(6f 101) 


URITC<b«llR)  I,(TITLES(LT*K),K=1«80) 

L1NES=5 

WR1TE(8»11»)  K. (DATAST J*H) . H= 1 .3 ) . KP 1 . < DATA S ( ) .M=A .6) . 


1 KP2, COATAS<J«n) ,n=7«9> 
GO  LINCS=LINES*1 


J=N05«9*20«H2 

K=N03*3*1 

lFt«3«H  Qi?.]X.  Efi  »N  p.LSJ_«.R  I TC(6.115)  UtLELilA  S(  J*H) 

KP1=K»1 

IF ( (3*ND3«2).EQ.NPTS1  WRITE(6,116)  K, (D ATA S ( J»M ) ,P= 1 , 3 > « 

LINES=LINES*1 
70  CONTINUE 

RETURN 

EnO 


ODOOAOlO 

OOOOA020 

00004030 


00004C40 

00004050 

00004060 


00004070 

00004080 


00004100 
0000411  0 

CL0QJH12.0 

00004130 

0O0C4140 

00004160 

00004170 

00004180 

00004190 


G LEVtL  21  NORM  DATE  = 78102  10/47/AB 


SUBROUTINE  NORN (DATA vISTART tNPTSt TNORN) 

DIMENSION  OATAIII 

TNORH=0. 

B0OOA41O 

OOOOAA20 

0000*430 

DO  20  1=1«NPTS 

00004440 

LsII-1>*3»ISTART-1 

00004450 

DO  20  J=1 t3 

000044GO 

20  TNORHsTNORN*  DATA ( L *>J  1*  *2 

00004470 

TNORR=TNORM/3./NPTS 

0000*4fl0 

TNORH=SQRT(TNORM> 

00004490 

RETUR  N 

00004500 

END 

0C004510 

170 


I 


G LEVEL 


CURFIT 


DATE 


78102 


10/A7/48 


21 


subroutine  CURFIT {FfXiN*MtCOEF  <OETtSIGtNH> 

00004520 

IMPLICIT  RE*L*8  <*-HtO-2) 

00004550 

-OIMENSIQN  A(9.91 tB(9)iX(l) iFI 1) iCOEFl  1) 

00004S4Q 

00  10  1=1. N 

00004550 

B( 1>:0.0 

00004580 

DO  10  J=1 .N 

nfinn4K7n 

10 

A( I,J>=0.0 

00004580 

DO  20  1:1. n 

00004540 

flO  iD  J=1.N 

QO&O4I5OO 

I>1 >.N 

00004610 

BlU>:P(J>.F(I).X(OJ) 

00004620 

..DO  20  L = 1.N 

00004620 

LLsL.«I-l)«N 

00004640 

20 

A(L.J>=A(L.J)*X<UJI*X(LL) 

00004650 

CALL  EQNSOLIA.B.N.COEF.OET.NH) 

00004660 

SIG=0.0 

00004670 

00  AC  I =1  .N 

00004680 

_SU.M=0.0_ 

00004A9n 

00  20  J=1.N 

00004700 

I-l ).N 

000Q4710 

2D  SUM=SU«»COEF(  J1  .X  f JJ) 

OQOO  4720 

SUH=SUM-OET 

00304720 

*0 

S 1G=SIG*SUM.SUM 

00004740 

.SIG=0SCRHS1G/M) 

OCG04750 

RETURN 

00004760 

ENO 

00004770 

G LEVEL 

21  SIGNF  DATE  = 78102 

10/97/48 

SUBROUTINE  S IGNF « COFB tNCOFS»DET ,*PTS, YPTS ,ZPTS»NPTS,I FI T 1 
REAL*8  VNORHC  3) VECT C 3> *X NAG • VP A GtOSQRT 

REALMS  COF0<  D.OET.XPTSll) tVPTSl 1 )»ZPTS< 1) 

00009780 

00009781 

00009790 

DIMENSION  IFIT(l) 

COSANG=a> 

DO  20  L=ltNPTS 

00009791 

00009792 

00009793 

XMAG=0SQRT(XPTS(L)*XPTS«L)*TPTS«L)»YPTS<L) ♦ZPTS(L)*ZPTS«L) I 
IFlXMAG.LT.l.  OD-1  0>  CO  TO  20 

XVECT I1)=XPTSCL1/XNAG 

00009799 

00009795 

00009796 

XVECT<2)=YPTS(L)/XNAG 

XVECT <3)=ZPTS<L»/XPAG 

DO  1 1=1.3 

00009797 

00009798 

0000*799 

1 

VNORN<I)=0. 

1 = 0 

IFIIF IT(l).FQ.P)  GO  TO  2 

00009POO 

0000*801 

00009802 

1 = 1*1 

VNORM(1)=VNORM«1)*COF0(I) 

IFII.EQ.NCOFSI  GO  TO  10 

0D004R03 

0000*804 

00009805 

2 

IF(IFIT(2>.EQ.O)  GO  TO  3 

1=1*1 

VN0RH*21=VN0RMI2)*CQFB( H 

0000*806 

00009807 

00004808 

3 

IFCl.EQ.NCOFS)  GO  TO  10 

IF<IFIT(3».EQ.O)  GO  TO  * 

1 = 1*1 

0000*809 
00004810 
0 n 0 0 48  1 1 

4 

VN0RNf3)=VN0RN(3)*C0Fe(  I) 

IF(l.EQ.NCOFS)  GO  TO  10 
lF<IFlT<4>.Ea.DI  GO  TO  5 

00009812 
00009813 
0000981 4 

1=1*1 

VNQRM( 1>=VN0RN( l)*COFB< I)*YPTS<L) 

VNORM(2  >=VNORM<.0)*COFB<  I ) - YPTSIL) 

00009815 

CCC04816 

OOOC9817 

5 

IFll.Ca.NCOFSY  GO  TO  10 

IFCIFITCSl.EG.C)  GO  TO  6 

1=1*1 

0030481R 

QC004819 

00004820 

VNORM ( 1 l = VNOBM( 1 ) *CCF  B< I) ‘EPTS (L ) 
VN0RM(3)=VN0RM(3)*C0FB(I)*XPTS(L) 

IFd.EO.NCOFS)  GO  TO  10 

00004821 

0000*822 

00004823 

6 

1F<1FIT(6).E6.0>  GO  TO  7 

I =1*1 

00004824 

00009825 

z_ 

VN0RN<3)=VN0RM<3)*C0Fb(I1*YPTS<L» 

IFd.EO.NCOFS)  GO  TO  10 

IFdFl  T(71  .EQ.fll  GO  TO  8 

00009827 
00009828 
noon  4P2Q 

1 = 1*1 

VNORM ll)  = VNORN  d > *2.*C0FBII)*XPTS<L) 

IFd.EO.NCOFS)  GO  TO  10 

00004830 

0000*831 

00009832 

A 

IFdF  IT(8).EQ.0)  go  to  9 

1 = 1*1 

VNORM12 )=VNORMt  2)*2 ..COFB  d )*YPTS  <L) 

00009833 

00009834 

00004835 

9 

IFd.EO.NCOFS)  GO  TO  10 

IFdFIT(9).EO.O)  GO  TO  10 

1=1*1 

00009836 

00009^37 

a00.fl9ft3B 

VN0BHI3)  = VN0RM«3)  ♦2.*C0Fe<n*ZPTS«L) 


1 
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10 

V(1*6=DSaRT(VN0RN<l)*VN(»N<ll*VN0RHC2>»WN0RM<2»*VN0Rm3)» 

1 VN0Rn(3>> 

IFtVMAC.EO.OI  GO  70  1000 

0000A8A0 

00004841 

00004842 

11 

DO  11  I=lt3 

VNORHC  I>  = VNORHU)/VnAG 

DO  12  1 = 1.3 

000048AS 
0000  4844 
00004845 

12 

20 

COSANG=COSANG«VNORniI>*XVECT( I) 

C0N7INUE 

IFCCOSANG.EQ.OI  GO  TO  2000 

00004846 
0000  484  7 
00004848 

28 

IF  (COSANG. GT.O)  RETURN 

DO  28  J=l.NCOFS 

COFB< JI=-COFB< J> 

00004849 

OOC048SO 

00004851 

1000 

0ET=-0ET 

retur  N 

ORITE<6.200) 

00004652 

00004853 

00004854 

2ca 

20C0 

FORNAKIHI .•BAGNI TUCE  OF  NORMAL  VECTOR  IS  ZERO  IN  SIGNF*) 

STOP 

URITE(G.20n 

00004855 

00004856 

00004857 

2C1 

FORMATdHO. 'WARNING  TEST  IN  SIGNF  SHOWS  VECTORS  PERPENOIC ULAR » > 

return 

END 

00004856 

00004850 

00004860 

G LEVEL 

21  EQNSOL  DATE  = 78102 

10/47/A8 

SUBROUTINE  ESNSOL ( A .8 « N.X tDET*NhOHO > 

00005000 

INPLICIT  REAL*B  ( A-H,0-S,U-Z1 

00005010 

REAL'S  TOL 

00005020 

DIMENSION  A(9«9>,B( 1) <X (1) 

00005030 

COMMON/TOLRS/TOLEQN«TOLSNT«TOLQPT,TOLINT,TOLCKl,TOLCK2t 

000050  40 

1 T0LC:K4 .TDELT  , T OLC IR LL ANGS^j TOLPLT 

00005050 

T0L=T0LE3N 

00005060 

NR  ANKzN 

00005070 

00  10  l=ltN 

10 

X tn^o.o 

00005090 

Iso 

OC005100 

I COLsO 

00005110 

12 

1 = 1*1 

00005120 

IC0L=IC0L»1 

00005130 

IP1=I*1 

00005140 

14 

AMAXI =X( I 1 ICOL) 

0C0C5150 

L:1  ! 

00005160 

IF(lPliGT.N)  GO  TO  150 

I1000S170 

DO  15  KrIPl.N 

00005180 

IF(OABS(A(K,ICOL)>>GT.OABSIAMAXn)  L=K 

00005190 

IS 

AMAXtrACL.lCOL) 

150 

I F(0ABS( AHAXI ) .OT. TOL)  GO  TO  16 

00005210 

X < ICOL ) = 1 .0 

00005220 

NRANKrNRANK-1 

00005230 

IF(ICOL.EO.N)  GO  TO  *5 

00005240 

ICOL  = tCOL*l 

30005250 

GO  TO  15 

00005260 

16 

IFIL.EQ.I  > GO  TO  25 

0000  5270 

DO  20  Jxl.N 

303052B0 

S VAP  = X 1 I . J> 

00005290 

A( I,J)=A(L«J> 

00005300 

20 

A(L.J>=SUAP 

0000531 0 

SWAP =8(1) 

00005320 

B ( I >=8(L> 

00005330 

S (L>  =SWAP 

003053SO 

25 

XDTV=AI7iIC0L) 

00005350 

00  30  J=IC0L<N 

OC005360 

50 

A( I,J)rA( I«J)/XDIV 

00005370 

00005380 

IFIIFI.GT.N)  GO  TO  *2 

00005390 

00  «0  KsIPl.N 

00005400 

XMUL=A(KtICOL> 

30005*10 

DO  35  JslCOL.N 

00005*20 

35 

A( Kf J)sA( K*J)-XPUL«A(I «J) 

00005*30 

SO 

B «K»=P«X)-XMUL*B« I ) 

00O05SS0 

IF ( ICOL.LT.N)  GO  TO  12 

00005*50 

«2 

IFTNRINK.EQ.N.OR.NHOMO.EQ.O)  X{ICOL):B(I)/A(I,ICOL> 

00005*60 

55 

00  50  I=2*N 

OOOQ5SRO 

K=N-I*1 

00005490 

1F(X(K).CQ.0.0)  GO  TO  AT 

00005500 

1 JUMP  = 1JUMP-1 

00005510 

174 
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♦ 7 

GO  TO  50 

IF(NRANK.Eo.N.OR-NHOHO.EQ«0)  x<k>=b«k> 

KP1  = K»1 

00005520 

00005530 

00005540 

98 

50 

00  48  J=KPltN 

X«K)=X(K)-A(K-IJUWPtJl**(J) 

CONTINUE 

00005550 

00005560 

00C05570 

03005580 

00005590 


EN 


SUBROUTINE  CFITS ( 0 A TAG »0ATAS« I STA RT tNP TStCOEFSB • ICOP • S IGtXCURF • 

1 FCURFl 

00005620 

00005630 

OIHENSION  SIGnA<8>«DETERN(e) 

REAL. 8 XCURF(  1)  ,FCURF(  1 ),C0FB(9).  SIG.DET 

REAL.a  XPTSI30).YPTS(30».2PTS(ia) 

00005650 

00005660 

00005661 

lOEFLTsOATACd  >.0.5 

00005670 

IPASS=0 

00005680 

CALL  NORMLOATAS.ISTART.NPTS.TNORH) 

00005690 

L = 1START-21 

00005700 

DO  5 1=1.9 

00005710 

5 

IFIT(n=OATAS(L*I)  *0.5 

00005720 

IFIIDEFLT.EG.O)  go  to  8 

00005721 

DO  6 1=1.3 

00005722 

6 

1 F I T H ) = 1 

00005723 

8 

CONTINUE 

00005724 

NCOEFS=0 

00005730 

DO  10  1=1.9  _ 

(UQ1157_4.Q_ 

10 

IFdFITd  l.EQ.l)  NCOEFS  = NCCEF  S.l 

00005750 

if(npts.lt.ncoefs)  go  to  5C 

00005760 

l5i__CaNTtNJUC  _ 

00005770 

IP*SS=IP*SS*1  00005780 


DO  35  Lrl.NPTS  00005790 


LU^=<L.-J  J ‘NCOEFS 

. . 0 0 0 0 5800 

1 J=«L-1 )*3*ISrART 

00005810 

X V=0AT AS ( IJ)/TNOR  M 
YVsOATASdJ^n/TNORH 

00005820 

00005830 

2V  = OATAS{  IJ.O/TNORf* 

00005840 

XPTSIL) =XV 

TP.lJLlLlsTV 

00005841 

..  _ _CiUlQ5a^2_ 

ZPTSIL) =ZV 

00005843 

I =0 

IF  ( IF  I T ( 1 > .FG.  C ) CQ 

T Q 3 1 0 

OOC05860 

nCQCSflTD 

1=1.1 

0G0Q5880 

XCURF(LJ.I)=XV 

310  IF«1FITI2>.EQ.0)  GO 

TO  311 

00005890 

00005900 

I =1.1 

00005910 

XCURF«LJ.I)=YV 

311  lFIIFITt3>.EQ.O)  GO 

TO  312 

00005920 

OCC05930 

I = I‘l 

00005940 

XCURF  <L J*I ) = ZV 
_JU. LF<  IF  ITI*  ) .ECi.O  ) GO 

TO  313 

000C5950 

00005960 

1 = 1.1 

00005970 

XCURr<LJ^I)2XV*YV 

SIS  IFdFircS)  «CQ>0)  GO 

TO  31* 

00305^80 

occe5®®o 

1 = 1.1 

XCURF(LJ.I)=XV.2V 

31*  IFllF  IKbl.EQ.C)  GO 

TO  315 

03C06C00 
G000601 0 
00CC6G20 

1=1.1 

00006030 

XCURF«LJ.I)=YV.ZV 

TO  316 

00006040 

00006050 

1 


I=I»1 


00006060 


6 LEVEL  21 


CFITS 


DATE 


78102 


10/A7/A8 


XCURF<LJ*I)=XV*XV 

90006070 

316 

1F(1FIT(8> .EQ.0>  GO  TO  317 

00006080 

1=1*1 

■HHBFrrnni 

XCURFCLJ*I>=YV*YV 

00006100 

317 

IF<IF1T«9).EQ.0)  GO  TO  35 

00006110 

1 = 1*1 

XCUHF(LJ*1)=ZV*2V 

00006130 

35 

FCU(<F(L)=1. 

00006190 

00006160 

CALL  SIGNF (COF 6 ,NC OEF S,DE T t XP TS» YPTS. 2P TS t NPTS , I F IT ) 

00006170 

IF(IOEFLT.EQ.O)  go  to  350 

00006180 

IFI  IPASS.rQ.9)  GO  TO  35  0 

DD0Q619Q 

SIGMA(IPASS)=SIG 

00006200 

DETERM( 1PASS>=CET 

00006210 

GO  TO  ( 318. 319. 320. 321. 322. 323»3r9. 325).  IPASS 

00006220 

318 

IF (NPTS.LE.3)  G3  TO  350 

00006230 

NC0EFS  = * 

00006290 

IFIT(7)=1 

00006250 

GO  TO  15 

00006260 

319 

IFITI  71=0 

00006270 

iFiT«e»=i 

00006280 

GO  TO  15 

0 00  0 629  0 

320 

I FIT«8)=0 

0000  6300 

IFIT<9>  = 1 

GO  TO  15 

00006320 

321 

IF<NPTS.LE.9)  GO  TO  325 

00006330 

NCOEFS  = 5 

OB3063-5J- 

IFIT(7)=1 

00006350 

IF1T(8)=1 

00006360 

• 

IF  IT(9)=0 

00006370 

GO  TO  15 

00006380 

322 

IFIT«8)=0 

00006390 

IFIT(9)=1 

■■■VFITRnil 

GO  TO  15 

00006910 

323 

IFIT(7)=a 

00006920 

IFIT(8)=1 

GO  TO  15 

00006990 

329 

IF«NPTS.LE.5>  GO  TO  325 

00006950 

NCOEFS=6 

00006960 

IFIT( 7)=1 

00006970 

GO  TO  15 

00006980 

■■HHTTITXEFI 

3250 

SlG=SIGnA (1) 

00006500 

L = 1 

C0C065! 0 

NCOEFS=a 

0CC06520 

IF(NPTS.EG.9)  NC0EFS  = <1 

00006530 

IF(NPTS.Ea.5)  NC0EFS=7 

0C006590 

DO  326  I=2.NC0EFS 

0C0C655C 

IF (AbS(0ETE8P( I) > .LT.O.S.ANO.IPATCH.CQ.O)  GO  TO  326 

00006560 

1F(SIGMA(I>.LT.SIG)  l=i 

00006570 

IF(L.EQ.I)  SIG  = SIGHA(n 

00006580 

326 

CONTINUE 

00006590 

G LCVEL  21  CFITS  DATE  = 78102  10/A7/«fl 


IF< INATCH.EQ.O. AND.SIG.lt. 1 .00-081  GO  TO  327 
IFdNATCH.EQ.l.ANO.SlC.LT. 1.00-061  GO  TO  327 
IFtINATCH.EQ.il  GO  TO  60 

00006600 

00006610 

00006620 

INATCH=1 

00006630 

GO  TO  3250 

000066*0 

^71 

TF  TT(  7)sQ 

QQ0066SQ 

IFIT(81=0 

00006660 

IFIT(9l=0 

00006670 

— Mimii 

IFtL.GT.l)  NCOEFSr* 

00006690 

IFtL.GT.*)  NC0EFS=5 

00006700 

IFtL.EQ.S)  NC0EFS=6 

00006710 

lPASS=a 

00006720 

GO  TO  ( 15. 328  *325.330.331,332 . 333»33«1  , L 

00006730 

328 

IFIT(  71=1 

000067*0 

CO  TO  15 

00006750 

32  9 

1FIT(8)=1 

00006760 

Gn  TO  

nnnftf.77n 

330 

IFIT(9)=1 

00006780 

GO  TO  15 

00006790 

1 

TF  TT I 7 1 = 1 

IF IT<81=1 

00006810 

GO  TO  15 

00006820 

332 

IFIT(7)=1 

00006830 

0000  68*0 
00006850 
00008880 
00006870 
00006880 
0000  6890 
00006900 
00006910 
Q0QQ6920 
0C0SC930 
000069*0 

00006960 
00006970 
00006980 
00006990 
000  0 7000 
0000  7010 
00007020 
00007030 
3C007;«0 
00007050 
0000  7060 
00QQ  7070 
00007080 
00007090 
0DQ071DD 


_li3_ 


rFIT<9)=l 
GO  TO  lb 
iFlT(B)=l 


IFIT(9»:l 
GO  TO  15 
33* IFIT<7)il 


IFITf8)=l 
IFIT<9)=1 
GO  TO  15 


350  CONTINUE 

L=( ICOF-1 )*10 
J=Q 


DO  *0  1 = 1 >10 
COEFSB(L*I )=0. 


IFdFIKI  l.EQ.O)  GO  TO  *0 
J = J*1 

36  C0EFSB<L*H=DATAS(1ST*RT>I-12) 
1F«C0£FSB<L*I>.NE.0.)  GO  TO  *0 
IFd.EQ.lO)  CO  TO  38 
C0EFSB(L*I)=C0F6< Jl/TNORN  


IF  ( I .LE.3 ) GO  TO  *0 
C0EFS8(L»I>=C0EFSB (L»I)/TN0R" 
GO  TO  *0 


38  COEFSB<L*I)=-OET 
*0  CONTINUE 
itiuAm 


5 


URITE(6.100>  ICOF 


00007110 


STOP ^ PIH07I*0 

60  yRlTE<6,110>  IcftF  00007150 

no  FORM/I  TflHl  .•SOUNDARY  POINTS  FOR  REFLECTIVE  SURFACE  NO.  ’.IS,  00007160 

1 » DO  NOT  FORW  A CURVE-FIT  WITHIN  ERROR  TOLERANCE*/ ^ 00007170 

2 IHO, ♦PROGRAM  WILL  CONTINUE  WITH  BEST  FIT  POSSIBLE*)  00007180 


GO  TO  327  00007190 

ENO  00007200 


179 


i 


G LCVEL 

21  0UTPT2  date  = 78102 

10/47/48 

SUBROUTINE  OUTP  T2  t COEF  B tS  I G .N  SE  TS  » TI  TEES, LOOP  ) 

DIMENSION  COEFB (1 ) .TITLES(1>,SIG( 1) 

•RITE <6.1 02)  <TITLES(MM).HH=l,80) 

00007210 

00007220 

00007230 

1C2 

103 

FORMAT! IHl.BOAl) 

WRITE(6.ia3) 

FORMAT!  lH0.<t6X,»*«*  RESULTS  OF  CURVE-FIT  ANALYSIS  ...♦) 

00007280 

00007250 

00007260 

10 

GO  TO  <10.50).  LOOP 

CONTINUE 

URITE<6.110) 

00007270 

00007280 

00007290 

no 

FORMAT! 1H0.2X.'.. FENCE  DATA*) 

LINeS=5 

DO  20  1=1  .NSETS 

00007300 

00007310 

00007320 

LL=!I-1)*A 

IF(!LINES*6).LE.50)  GO  TO  15 

VRITE!6.102)  !TITLESIMM).MM=l.flO) 

00007330 

00007380 

00007350 

WRITE!6.1 03) 

UR  ITE (6.110) 

L INrS=5 

0000  7360 
00007370 
0C00738Q 

15 

112 

URITE!6.112)  I 

FORMAT! 1H0.5X. ••FENCE  NO. *.13) 

LINES=L1NES*2 

00C0739C 

00007800 

00007810 

121 

URITE(6.121)  !C0EF9(LL*L).L=1 .A) 

FORMA T( IHO .OX . •VALUES  FOR  COEFFICIENTS  OF  SURFaCE^/IH  .IIX. 
1 •(!)  X COEF  •.EU  .A.5X,  • (2)  Y COEF  • .£  12  .<• . 5X  , • < 1)  2 

00007820 
00007830 
COEF  ^000  0 788  0 

123 

2.E12.A.5X.^(‘i)  CONSTANT  •.E12.*) 

URITE!6.123)  SIG<I) 

FORMAT!  1H0.9X. •ROOT  MEAN  SOUARE  ERROR  OF  FIT=^.E12.A) 

00007850 

00007860 

00007870 

20 

LINES=LINES*5 

CONTINUE 

RETURN 

00007480 
0000  7490 
00007^00 

50 

117 

CONTINUE 

WRITE!6.117) 

FOR  MA  T(1H0.2X.^««  RJLf  LC.C  TIVE  SUR  FAC£_  T A • ) 

000C7510 

00007520 

00007530 

LIN£S=5 

lt=o 

DO  60  1=1. NSETS 

00007580 

00007550 

00007560 

LT=LT.80 

LL=!I-1)*10 

IF!!LINCS*9).LE.50)  GO  TO  55 

00007570 

00007580 

00507590 

URITE(6.102)  (TITLES! MM >.MM=1.R0) 

URITE!6,103) 

UR  ITE  (6.117) 

00007600 

00007610 

00007620 

55 

UR 

LINES=5  00007630 
URITE!6.118)  I. !T ITLES !LT*MM) .MM  = 1 .80  ) 000076A0 
FORMAT! 1M0.5X. ••REFLECTIVE  SURFACE  NO . • . I 3 . 5X . • I DENT  IF IC A T I ON : •.  00007650 

1 80A1 ) 

L INES=LIN£S»2 

•RITE (6. 122)  (C0EFF(LL*L).L=1 .10) 

00007660 
000  0 7670 
000  0 7680 

122 

F0RMAT!lHa.9X. ‘values  FOR  COEFFICIENTS  OF  SURFACE*/1H  .IIX. 

1 COEF  •.E12.A.5X.^!2)  Y COEF  • .E 1 2.8 .5X . • ( 3)  Z COEF  •. 

2 IH  .11X.*(8)  XY  COEF  •.E12.8.5X.*!5)  ¥2  COEF  •.C12.8.5X. 

•!1)  X 00007690 
E12.8/  00007700 
0000771 0 

3 •!&)  YZ  COEF  •.E12.8/1H  .11X.^!7)  XSQ  COEF  •.E12.8.5X. 

00007720 

G LEVEL  21 


0UTPT2 


DATE 
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4 •iSi  TSO  COEF  •,E12.4*5X» ‘IS)  ZSO  COEF  ’.EIZ.^/IH  *10X* 

5 •ElO)  CONSTANT  ’.EIZ.A) 

WRITFIA>I23>SIG(II 

00007730 

000077*0 

00007750 

LINESsLINES*8 

00007760 

60  CONTINUE 

00007770 

R€ruf(N 

0000  7780 

END 

00007790 

G LCVCL  21 


GLINT 


GATE 


7«102 


10/A7/A8 


SUBROUTINE  6L I NTl DAT AC t OA TA St NS fC OEFSB .GAnmA t ALPHA2«BETa2 • INTREL*  00007800 
1 ISHAOUtINTERFtOWORKStOATAFtNF)  00007810 


DIMENSION  OATAC (1 ) tOATASC 1) tMS( 1) tCOEFSBC 1 > 
D IMENSION  DWORKSCl > tO A TAF ( 1 > t HF ( 1 > 


1 RFLTNB(S)tRFLTN( 31 

DATA  ONEt ZERO t rad/ 1. OEOOtO .aEaat57.2SS7Be00/ 

SUNM)g-COS<GAMM  A/RAD) 

SUN! 21=0. 

SUN( 3I=-SIN«CAMMA /RAO) 

C 02=j^A-T  A CJ  1 5 » *-C  A.T  AC  U 5 ) • 1 « * • _ . 

TMETAsDATAC  CA  ) 

PHI=0ATAC(5) 

CSTh=COS<THET  A/RaD) 

SNPH=SIN(PHI/RAO) 

_CiPM.=  COS<PH^l/RAC)  _ . _ 

SPSTsSNPH«SNTH 

SPCTrSNPHtCSTH 

CPST=CSPHtSNTH 

CPCT=CSPH*CSTH 

DIRCCS«1.3):-SNTh 

DIRCOS«?t  3 ) = - 

DIRCOS<  3.3)=CPCT 
OPSI=DATACni) 

(lEiI=DAIAC(?)  *0.5 

NFEMCE=DATAC<2r)*0.5 

NPANtL=0ATAC<23)»Q.5 

(iA.Sia - , . 

DO  ICO  I=1,NPANEL 
M10  = I I-1>  »10 


00007830 

00007840 


NPTS=0ATAS<KI )*0.5 
PSIKPiOATACI  10).-OPSI 

DO  bO  J=l.NPSl  . 

PSIKPrPSlKP»CPSI 

FSI=PSIKP 

IF<PS1KP  . L T . - 1 7T  . ) PS  I = -1 79  . 

IFCPSIKP.GT.179.)  PS1  = 179. 

IFfPSIKP.GT.-91.AND.PSIKP.LT.-89.)  PS  I = -89  , 

lECPSlKP.LT.Sl.ANO.PSIKP.GI.aS.)  PSI^eR. 

SNPS  = SINlPSI/RAO) 

CSPS=COS(PSI/RAD) 

DlfiCOStlti)=-CSTH*CSPS 

OIRCOS(1.2)=-CSTh.SNPS 
01 RCOSC  2. 1 )=SPST.CSPS-SNPS ‘CSPH 
DIRCOS  C 2.2) =SPS  T .SNPStCSPS .CSPm 
OIRCOS«3.1):-CPSTtCSPS-SNPStSNPH 
DIRCOS<3.2)=-CPST.SNPS  *SNPM.CSPS 

DO  20  M = 1.3 

SUNB(N)=0. 


00007860 

00007870 

00007880 

00007890 
00007900 
QjjQ  0 7910 
00007920 
00007930 
0QD07940 
00007950 
00007960 

0iJLDJ9-IQ 

0 00  0 798  0 
00007990 

00008000 

00008010 
00008020 
000  08030 

00008040 

occoeoEC 

atimeDiiii 

00008070 

00008080 

ODDCEOSC. 

00008100 

00008110 

30008120 

00008130 

00008140 

15.0  Oil  50. 

00008160 

00008170 

00008180 

00008190 

00008200 

0000 8210 

00008220 

00008230 

B00.B.B2.4Q 

300082E0 

00008260 

00!!i?8270 

00008280 

00006290 

Ofl-Q.Mm 

0000831 0 


182 
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DATE 
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10/47/A8 


r 


20 

00  20  K=1  ,3 

SUNB  (M)=SUNB(N>*OIRCOS(f1tK}*SUN(K) 

00  75  L=1.NPTS 

00008320 

00008330 

0a0083A0 

HAS=MAS*1 

00008350 

INTRFLCNAS)=0 

00008360 

1SHAOWIMASJ3C 

INTERF<MAS)=0 

00008380 

M3=(L-1)»3*KI 

00008390 

00  10  M=1.3 

0Q008A0Q 

10 

XB<NI  =0ATAS«N3»M) 

00008A10 

IF(NFENCE.GT.O)  CALL  S HADOU <0  AT AF , MF ,NFENC E tX 8 t SUNB 1 1 SHA OU * 

00008420 

1 MAS.DUQRKS) 

0Q00BA30 

VNORHB<1>=COEFS8CN10*1)*COEFSB(N10*A)«XB|2)*COEFS8(N10*5)*XB<3>* 

00008440 

1 2.«C0EFS8 «Nl 0*7) .XB { 1 ) 

00008450 

VNORHB(2)=COEFSB«N10»2) ♦COEFSB(N10«*>*xe«l  >*COEFSB(N10*6»»Xe<3)* 

00008A60 

1 2.*COEFSB<N10*8)*X8(2) 

00008A70 

VNORM8(3>=COEFSP(N10*3>*COEFSB(N10*5)»Xfi(l)*COEFS8«NlO»6)*XBC2)* 

00008480 

1 2.*COEFSB«N10*=)*XB<3) 

00008440 

ASN=0. 

00008500 

ABSNSO=0. 

00006510 

00  30  H=1 ,3 

ASN=ASN»SUNB(M) •VNORWe  <M> 

00008530 

30 

ABSNSC  = A0SNSG*-VNORr^8  CM)  **2 

000085AO 

ASN=ASN/SCRT( ABSNSO) 

Q000fl55C 

IF(ASN.GT.0.0017A5)  INTRFL(HAS)=1 

00008560 

00  AO  M=l,3 

00008570 

A 0 

RFLTNB«M>=SUNB<M>-2.*ASN*VN0RMB(f*)/SQRT(A8SNSQ) 

00008580 

IFCNFENCE.GT.O)  CALL  I N TF  ER  < 0 AT  A F t BF  . NFENC  E , X P.  , R FL  T NB  , INTER  F, 

00008590 

1 mas,duorks> 

00008800 

00  50 

OOOCR610 

RFLTM(M)=0. 

00008620 

X (B)=0. 

00008630 

0050  K=1 ,3 

000Q86A0 

RFLTN(M)=RFLTN<M)*0IRC0SIK.W)*RFLTNB«K) 

00008650 

50 

X(B)sX«M)*OIRCOS<K,f«)*XB<X) 

00008660 

CALL  OUAOP T< ONE. ONE.2FRO.2ER0.2ERO,2ERO,ZtRC, ZERO tZEROtOQ2t 

00008670 

1 RFLTNa).RFLTM<2>,RFLTN<3),X(l) .X(2),X<3>,XIl,TIl,ZIltXI2, 

00008680 

2 TI2»2I2. INTC 1 

00008690 

IF(INTC.EQ.l)  GO  TO  55 

■TTTmn 

CALL  SORNOT<RFLTN(l>,RFLTN(2>«RFLTN(3)«X(l>«X(2>tX(3>«XIl«YIltZIl«00006710 

1 ISHAO) 

00008720 

IF(lSHAD.EQ.l)  GC  TO  55 

X 1 1=X 12 

OOC087AO 

7I1=VI2 

00008750 

ZI1=ZI2 

00008760 

55 

CONTINUE 

00008770 

ALPHA2IMAS»=RA0«ATAN2(2Il.S(iRT<XIl*XIl*YIl*YIl)) 

00008780 

6ETA2(NAS>=RA0*ATAN2(YIltXIl) 

00008790 

75 

CONTINUE 

00008800 

60 

CONTINUE 

00008810 

■IIMIU!] 

RETURN 

00008830 

78102  10/47/48 


SUBR0U7INE  SH AOOy ( 0*7 AF .NF tNFENCC  «XBS« SUNB 1 1 SHA DWtMASfOUORKS)  00008BS0 

OmeNSION  DA7AF(1  >.HF(1>|ABS(1)  tSUNBdl  .ISHAOy  (1)  tOUORKSd)  «XBF13>000aBB60 


6PL=-SUNB( 2) 
CPL=-SUNB«  3) 


BlT*iWT^T*jyT 


7ME7A=0. 

IF< APL.NE.O.O.OR.CPL.NE.O.Ol  THE7A =A TAN2( APLt CPL > 


F.1^ 


SNTM  = SIN«  THETA) 
CSTH  = COS<  THETA) 
SNAL=SIN(ALPHA> 


CSAL  = COS(ALPHA) 

SAST=SNAL*SNTH 

SACT=SNAL  «CSTH 

CAST=CSAL ‘SNTH 
CACT  = CSAL«CSTH 

_iM£NCE 

KF=HF (I) 

NPTSF=OATAFIKF)*0.5 

TSF 

LI=«J-1)»3 

LK=(J-1).2 

QQ  10  K=1.3 

10  XBF(K)sOATAF(LI*K*l) 

CALL  PLANPT<iPLt0PL»CPL,OPL.SUNe<  l)»S<JNB<2)»SUNe<3)*0ISC,X«F«  1), 

1 XBFt2) .XBF(3) .XINTP.YINTP .ZldlPi 

CALL  S0RN0T(SUNR(1 ) • SUNB ( 2 ) t SUNB (3)«XBF(l)<XbF(2)fXBF(3><XlNTP, 

1 YINTP.ZINTP. ISH) 

IFIISM.ES.O)  GO  10.50 

20  OyORKS(LR»l)=XINTP«CSTH-ZINTP*SNTH 

OyORXS«LK*2)=-XINTP.SAST»YINTP*CSAL-ZINTP»SACT 


OOOOBBBO 

00008890 


90008910 

00008920 


000089A0 

00008950 

00008963 


00008970 

00008980 

00008990 

00009090 

00009010 

900090  30 
00009040 

00009060 
00009070 
00009080 
r .,09090 
00009100 
-■3BQa911tt 
0000P120 
00009130 
QQ3Q91AQ 
00009150 
00009160 


XC=XBS(1) •CSTH-XPS( 3)*SNTH 
YC=-XBS(1)*SAST*XBS«2)*CS4L-XBS<3)«SACT 

IF ( INCR .EQ.O)  GO  TO  50 
ISHAOy (HAS)=1 


00009180 

00009190 

00009200 

00009210 

00009220 
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subroutine  PLANPT(A*B«CtOtSL«  SN tS N«OI SCtXl »Y It 21 *X*T *2 > 
CONHON/TOLRS/TOLEQNtTOLSNTf T0LQPT*T0LINT«T0LCK1*T0LCK2« 

1 T0LCK»«0ELT.T0LC1R.ANGST2»T0L 

IF<AaS(SL>.GE.ABS(SN).ANO.ABS(SL> .GE.ABSCSnn  D 1 SC=A*B*SM/SL* 
1 CtSN/SL 

IFlABStSD.LT.ABSCSN) .ANO.ABSISM) .LE.ABS(SN>  > 01 SC  = A»SL/SNt 

1 a*sn/SN*c 

tF<ABS(SL>«LT.ABS(Sn>.ANO.ABS«SN>.LT.ABSCSM>>  0 ISC=A tSL/SN* 

1 BtC»SN/SH  

IFfABS(OISC>.LT.TOL)  RETURN 

IF (ABS«SL> .LT.ABS(SN) .AND.ABStSH) .LE.ABS(SN) > GO  TO  10 

IF  < ABS«SL  > .LT.  A B^S_<  SF  >_ . * NO  . A 8 S < S_N  )^.  L.I  .J.  B S«SM)  ) GO  TO  20 

X = <-0-B*ri»8tX 1*SH/SL-C ‘ZltC*  X 1 .SN/SL ) /0I SC 
Y = Y1*SN*«X-X1)/SL 

Z=Z1«  SN«( X-Xl ) /SL 

return 

10  2=<-p-A*Xl*A*ZltSL/SN-B*YltB*Zl*SM/SN)/DISC 

XsXl  «SL«<  2-21  )/SN 

Y=Y1»SN»<2-21)/SN 
RC  TtlR  N 

Q Y=(-D-A«X1*A«SL«T1/SP-C»21*C«SN«T1/SN)/(A«SL/SH«B*C«SN/SH) 
X=X1«SL*<  Y-Yl  )/SN 
Z=ZltSN*«Y-Yl)/SN 

RETURN 


00009270 

00009280 

00009290  ^ 

00009300 

00009310 

00009320 

00009330 

0000934C 

00009350 

00009360 

00009370 

00009380 

00009390 

OC009900 

00009910 

00009920 

00009930 

00009990 

00009950 

00009960 

00009970 

00009980 

00009990 

00009500 


END 

f 


00009510 


SUBROUTINE  SORNOT < SL t SN «SN*X1 «T 1 tZ 1 tX • Y «Z • ISHAO ) 
C0HH0N/T0LRS/T0LEQN,T0L»T0LQPT«T0LlNT,T0LCKl.T0LCK2t 
KA.DELT.T0LCIR.ANGST2.T0LPLT 


1SHA0=2 
SL1=X-X1 
Ni=Y-y 


SN1=Z-Z1 

IF<(ABS(SLl)»A8S(SNl)*ABS(SNl) I.LT.TOL)  RETURN 


IF(ABS(SL>.GT.TOLI  DIR=SL1/SL 
1F(A6S(SP>.GT.T0L>  DIR=DIR»SM1/SM 

IF(ABS<SN)  .GT  iT.C Ll_(U SN 

ISHaO  = 1 

IF(DIR.LT.Q.O)  ISHAC=0 


00009520 

00009530 


00009550 

00009560 

00009570 


00009580 

00009590 


00009610 

00009620 

00009630 

000096A0 

00009650 
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DATE  = 78102 
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SUBROUTINE  INTFER  < DATAF tHF t NF ENCE tXBS tRFLTNB • INTERF «n AS .OUORKS > 
OINENSION  OAT  AF(  1 1,I1F(1)  ,XBS(  1)  tRFLTNB  (1>  (INTERF  (1)  tOUORKSd)  « 

1 XBF(31 

APL=-RFLTNB( 1) 

BPL=-RFLTN8<2) 

CPL=-RFLTNB(3> 

0PL=-APL*XBS«l)-BPL*XBSI2)-CPL*xest5) 

THETA=0. 

IFIAPL.NE.OtO.OR.CPL.NE.O.O)  THETA=ATAN2« APLt CPL) 

ALPHA=ATAN2<BPLtSQRT(APL*APL*CFL*CPL) > 

SNTHrSINI THETA) 

CSTh=COS( THETA) 


00009680 

00009690 

00009700 

00009710 

00009720 

00009730 

00009790 

00009750 

00009760 

00009770 

00009780 

00009790 


► 


1 

f 


I 

I 


SNAL=SIN( ALPHA) 
CSAL=COS<ALPHA) 


SACTiSNAL  ‘CSTH 


cast=csal«snth 

iA  Cl^C  S A L*  C S T h 

DO  50  I=ltNFENCE 
KF=HF<I ) 

NPTSF=OATAF(KF)*0 .5 
00  30  J=ltNPTSF 
LI=( J-l)*3 
LKr( J-1 )«2 


10 


00  10  K=lt3 

XBF<K )=DATAF(LI*K*1  ) 

_tAl.L_ELAJ^RJLAPLjB.PLtCPLiDPL_iRFLINBiLJ_tfi.£LXfiaiJJJ.t_R 

1 XBF(l),X8F<2),XEF<3)iXlNTP,riNTPtZINTP) 


00O09B0O 

00009810 

00009830 

30009890 

A!UUllSa.5J!. 

00009860 
000  09870 
00009880 
00009890 
00009900 
0CC0991  0 
00009920 


00009930 

00009990 

00009950 


CALL  S0RN0T(RFLTN8»l)»RFLTNB<2).RFLTN9t3)tXBF*l)tXBF(2)(XBF(3>» 

l_XIjaP.IlNTPt2IfiTPtlSH)_ 

IFdSH.EO.l)  GO  TO  50 
20  DuORKS(LK»1)=XINTP*CSTH-ZINTP*SNTH 

OMORKS< LK*2>=-XI NTP^SAST ♦TINTP«CS AL-ZINTP« SAC  T 

30  CONTINUE 

XC  = XbS  < 1)*CSTH-XBS<3) (SNTH 

TC=-XBS<1  > (SAST ♦ X B S_(  2 )_lC S A LrX  B S_(  3 ) (SACT 

CALL  CIRCLE* DWORKStNPTSF.XCtTCt IN CR) 

IF tINCR.EQ.O)  GO  TO  50 

INTERF(HAS)=1 


0000  9960 
imjl“97Q 
00009980 
00009990 

00010010 

00010020 

0001003C 

00010090 

00010050 

00010060 


GO  TO  60 
50  CONTINUE 
60  RETURN  _ 
END 


00010070 
00010080 
JU!i!  1.0031)1 


0001 0100 
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QUAOPT 
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5 


subroutine  aUAOPTC  AtBtCtOf  E«F*6«HvCK«CLtSL.SH.SN,Xl,Tl,Zl,XIl,YIl  ,00010110 

1 211, XI2,YI2, 212. INTO  00010120 

COHWON/TOLRS/TQLEaN.TOLSNT.TOL.TOLINT.TOLCKl,TOLCK2. 08810130 

1 T0LCKA,0ELT,T0LCIR,ANGST2,T0LPLT  OOOIOIAO 

INTC  = 0 00010150 

IFiABSCSLl.LT.ABS(SN).AND.ABS(SM>.LE.ABS(SN>)  GO  TO  30 OCOIOIGO 

IF(A8S(SL> .LT.ABSCSN) .ANO.ABS(SN) .LT.ABS(Sn>>  GO  TO  50  00010170 

SA=A*B*<SN/SL) ••2*C*  <SN/SL) ••2»0*SN/SL*F»SM*SN/SL*«2»E*SN/SL  0001018  0 

SBr2.«B«Yl«SW/SL-2.«B«Xl«<SH/SL)««2*2.«C»21«SN/5L-2.«C«Xl« .00010190 

1 «SN/SL)«*2*0»Y1-P»X1*SM/SL*E*Z1-E*X1*SN/SL»F«Z1*SN/SL-F*X1*SN«  0001 0200 

2 SN/SL**2«F*Y1*SN/SL-F*x1*SM*SN/SL**2«G*H*  SM/SL«CK*SN/SL  0 0010210 

SC=B«Y1«Y 1«B«X1«X1 M SW/SL) ♦♦2*C«Z1«ZI-2.« 9«T 1«X1»SH/SL*C*X1«X1»  0001  0220 

1 (SN/SLS**2-2.«C*X1»Z1*SN/SL*F*Y1*21-F*X1*Y1*SN/SL-F*X1«Z1»SP/SL*  00010230 

2 F*X1  «X1  .SP^SN/SL  ••2*H.  Yl-M.X  1 * SM /SL ‘CK  ,21  -CK  *X  1 *SN/SL»CL  0001  02<t0 


1CQNTR.=  1 

OISC=SB«SB-»..SA.SC 
IFinlSC.LT.  SBSISAI/IO.)  RETURN 

IFiniSC.LT.  1.  OX_REJ  URPL 

IFCA8S«SA».6E.T0L)  GO  TO  10 
IF(ABS(S8).LT.T0L)  RETURN 

XI1=-SC/S8 

iNTCrl 

IF1IC0NTR.EC.2)  GO  TO  60 
IF<IC0NTR.EQ.3>  GO  TO  AO 


00010250 
00010260 
00010270 
0001 C2Hn 
00010290 
00010300 
00010310 
00010320 
00010330 
000103A0 


GO  TO  2 0 

10  XIl=-se/SA/2.»SQRT(OISC)/SA/2- 

XI2  = -SB/SA/2»-SQRT  1 Di^  )/iA/2._ 

INTC=2 

IF(IC0NTR.EG.2)  go  to  60 

IE  tL£.a!LIjL,I.(l,_31.  GO  J tL.A.Q 

20  Yll=Yl*SW*tXIl-Xl>/SL 
ZI1=21«SN*<XI1-X1>/SL 

TF fTNTC.FQ.l > RETURN 

Y I2rY 1»SM.<XI2-X1)/SL 
ZI2=Z1*SN*(XI2-X1)/SL 

30  SA=C*B*»SP/SN) .*2,A*«SL/SN),*2*F*SM/SN*E*SL/SN»D*SM*SL/SN**2 
Se  = 2.»e*Yl*SP/SN-2.«B ‘Zl*  «SP/SN> ••2*2.* A.X 1*SL/SN-2.»A*Z1* 

1 (SL/SN1««2*F.T1-F«Z1  -SPZS  NtEEXAHLLa  ♦SL/SN*Q  » X 1 . $M/SN-0  >2 1 »SH, 

2 SL/SN**2*0*Y1»SL/SN-D*21 • SM* SL/SN • *2 *0 K*H • SB/ SN*G • SL /SN 
SC=B«Y1*Y1-»B*Z1*Z1*(SP/SN)««2-2.*B*Y1*Z1«SH/SN*A.X1*X1-»A*Z1*Z1» 

L_J  SL  ZSN  T*  *2 -2  ..3  A • XI  *2.1  3 SL/SN»0«Y1  «X1- 0*  Y1  • 2 1»SL/SN-0«X1«21«SM/SN* 

2 0«Z1  •21*SP*SL/SN,*2*H*Yl-M',Zl*SP/SN*6«Xl-f.«Zl*SL/SN*CL 
ICONTR=3 

50  . T0_5 . 

AO  2Il:XIl 

YI1=Y1*SP«<ZI1-21 )/SN 

lLLL5Xi«aL.•J^ll,tZl)/SN 

IFdNTC.Ea.ll  RETURN 
Z12=XI2 


00010350 
00010360 
00010370 
00010380 
00010390 
OQC 1 OAnn 
0001 OAIO 
00010420 
0001 0430 
00010440 
00010450 
00010460 
00010470 
0301  0480 
00010498. 
00010500 
00010510 
00010520 
00010530 
00010540 
00010550 
00010560 
00010570 
8.0818.580. 
00010590 
00010600 


TI2=T1* 


0D01C610 


XI2=X1*SL*1ZI2-Z1>/SN 


00010620 


G LEVEL  21 
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DATE  = 78102 


10/A7/A8 


80 

RETURN 

SA=B»A«<SL/SN)**2»C»CSN/SN)*«2*0«SL/SH*E*SL»SN/SH*»2»F«SN/SH 

SB=2.»A»*1»SL/SH-2.*A»Y1»(SL/SN>»»2*2.»C»Z1»SN/SM-2.«C»Y1» 

00010630 

ooaio6Ao 

00010680 

1 |SN/S«)«*2»0**1-0*Y1*SL/SH*E**1*SN/SH*E*Z1»SL/SN-2.»E*Y1»SL*SN/ 

00010660 

2 SN**2*F*Z1-F«Y1*SN/SM*G*SL/SM»H*CK*SN/SM 

0001 0670 

SC=A*X1*X1*A*Y1*Y 1* (SL/SN> ••2-2.*A*ri*Xl*SL/S«»C*21 •21*C*Y1*Y1. 

000 10660 

1 fSN/SN)»*2-2.*C*Yl*Zl»SN/SM*E»Xl*21-E*Xl*Yl*SN/SM-E»Tl*Zl»SL/SH* 

00010690 

2 E»Y1*Y1*SL»SN/S««*2*G«X1-G«Y1*SL/SH»CK.Z1-CK»Y1*SN/SN*CL 

00010700 

1CONTRS2 

OOOlRTin 

GO  TO  8 

00010720 

60 

YIlsXIl 

0001 0730 

«I1=X1*5;i  •CYTl-Yl  »/SN 

00010790 

ZI1=Z1*SN*<YI1-Y1 >/SF 

00010750 

IFtINTC.EQ.l)  RETURN 

00010760 

Y12=XI2 . 

-00010770 

XI2=X1»SL«<YI2-Y1  )/SM 

OOC107BO 

ZI2=Z1*SN««YI2-Y1)/SH 

00010790 

RETURN 

000  10800 

EnO  0001  081C 
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0UTPT3 


DATE 


78102 


1O/A7/A0 


SUBROUTINE  0UTPT3 ( DA T AC tOA TAS*NS • TI TLES .G ANNA ,ALPHA2*BET*2 . 

1 INTRFL.ISHAOUtINTERF) 

DIMENSION  INTRFLd)  .ISHAOUCD.INTERFCl) 

0001 0820 
00010830 
00010840 

DIMENSION  DATAC(l).0ATAS<l>tMS(l)tALPHA2(l).BETA2<l).TITLES(l) 
LT  = 0 

MAS=0 

00010650 
000 1 0B60 
0001 0B70 

NPSI=DATAC (9)*0.5 

00010880 

0PSI  = 0ATAC<11) 

00010890 

NPANEL  = OATAC<  23)*0.5 

00010900 

DO  90  1=1 .NPANEL 

000  1 091  0 

K=ns< I ) 

0001 0920 

NPTS=DATAS(K)^0,e 

0001  0930 

VRITE<fatl02)  (TITLES(MM),MM=1,80) 

000  1 09«  0 

LT=LT*80 

OCO 1 0950 

WRITE (6,1 01)  I ,<TITLES<  LT*MN) ,NH  = 1 ,80) 

00010960 

1C2 

FORMAK  IHlfSOAl) 

00010970 

101 

FORMAT(1HO,51X,»,*,  SUN  GLINT  SIGNATURE  • • • • / IH 0 , 5X , • REF LECT I VE 

SUOOCl C9B0 

IRFACE  N0.*,I3,5X,»IDENT1FICAT ion:  dSOAl) 

00010990 

WR  ITE(6,1  03)  GAMMA,DATAC(  12)  ,OATAC<4)  ,D  A TA  C < 1 3 ) ,0  A TAC  (5)  , 

OCOllOOO 

1 OATAC(  14  ) ,OA  TAC(  15  ) 

00011010 

10  3 

FORMAT(1HO,5X.*SUN  ELE V AT  I ON* , F 1 2 . 2 , 76X , 'XROT • , F9 .2 / 

90011020 

1 IH  ,5X,*A7C  PITCH  ATTITUDE •,F7. 2 ,76X, 'TROT* ,F9 .2/lH  ,5X,'A/C  ROLL00011030 

2 attitude*, F8. 2, 7bX,»ZRCT«,F9.2/lH  , 1 06X , • D I STG • ,F 8 . 2 ) 

000  1 lOAO 

'WRITE<6,104) 

000 1 1050 

104 

FORMATdHO,  5X,*A/C  YAW  • ,20X, ‘BOUNDAR  Y POINTS  — BODY  AXES*,  7X 

, 00011060 

X 'STATUS  FLAGS  FOR  P0INTS*»1RX, 

00011070 

1 »GLINT*/1H  ,30X.*tSTATIONLINE,eUTTLINE, WATERLINE)*,  5X, 

9JJ.oiioe.o_ 

Y *INRFL*,4X,*I SHA  D* ,«X,*INTR4',11X, 

0C011090 

2 *AZIMUTH* ,3X , *ELEVATI0N*/ 1 H , 3 6X , * X • , 9X , • Y * , 9X , * Z * ) 

00011100 

L INFS:13 

OOMHIB 

PS1=DATAC  (10)  -DPSI 

00011120 

DO  60  J=1,NPSI 

00011130 

PSI=PSI*OPSI 

00011140 

DC  75  L=1,NPTS 

00011150 

MAS  = MAS»1 

000  1 1160 

N3=(L-1)*3'K 

00011170 

XR=DATAS(  N3'l ) *0  AT  AC  ( 12) 

QOOlllRO 

YB=DATAS(N3'2)*0ATAC(13) 

00011190 

ZB=DATAS(N3''3)40ATAC(14) 

00011203 

WRITE  (6,  10  5)  PSI,XB,YB,ZP,INTRFL(MAS),ISHADU(MAS)  , IN  TERF  (M  AS  ) » 

000 1 1210 

1 BETA2(MAS),ALPHA2(MAS) 

00011220 

105 

FORMAT(1HO,5X,F7.2,19X,F8.2,2F10.2,10X,I1,8X,I1,8X,I1,13X,F7.2, 

0001  1230 

1 4X,F6.2) 

000  1 1240 

LINES=LINES'2 

00011250 

75 

CONTINUE 

00011260 

60 

CONTINUE 

00011270 

90 

CONTINUE 

0001 1280 

RETURN 

000  11240 

END 

00011300 

AO-A061  322 


UNCLASSIFIED 


BOEING  VERTOL  CO  PHILADELPHIA  PA  F/G  1/3 

CANOPY  SUN  GLINT  EVALUATION  COMPUTER  PROGRAM. (U) 

SEP  78  F WHITE  DAAJ02-76-C-0061 

USARTL-TR-78-35  NL 


ri 

h , 

n 

fl 

1 

M 

LI 

1 

r^i 

fr-  1 

fl 

••  • 

L 

ri 

1 1 

1 

F 1 

n 

L 

1 

II 

mm\ 


f 


G LEVEL 

21  PINM*X  DATE  = 78103 

08/5A/AA 

SUBROUTINE  N INNAX  ( OATAC «0«TAS*MS«  ALPHA2 .BET«2*«PLOT  tBPLOTt 

INTRFL*  00011310 

1 ISHADUtlNTERF.IPLTOL) 

00011320 

nAfl  1 

1 APLOT(I)  ,BPL0T«1)  .INTRFL<l)«ISHAOW<l)tINTERFCl),IPLTOL<  1) 

00011340 

IRFLTN=OATAC<3)*0.5 

00011350 

MPAUri  -HATAr 

NPSI=0ATAC(9>*0.S 

00011370 

NPSIA=A«(NPSI-1)*1 

00011380 

TSTsl 

nnn i i xqh 

PASiO 

OOCllAOO 

LAS  = C 

OOOllAlO 

DO  70  in.NPANEL 

noci 1 44n 

KNrMSIU 

O0P11A70 

NPTSsOATASCKN>*G.5 

OOQllABO 

x-lil 

Kp=n 

00011500 

DO  5G  L=1,NPTS 

00011510 

HA<:=MA5:«1 

nnfti  1 «)9n 

IF(INTRFL<MAS>.NE.1RFLTN)  GO  TO  50 

00011530 

IFdSHADWCNASl.EQ.l)  GC  TO  5C 

OCOllSAO 

NP=NP*1 

00011560 

IFINP.GT.I)  GO  TO  10 

0001 1570 

91  0UrRrTA3£J*ACi 

onfl 1 i^®n 

kbnin=nas 

00011590 

eHIGH=BLOU 

00011600 

-JtfiiUJliBAS 

nnn i i4i n 

AL0U=ALPHA2( PAS) 

90C11620 

xanin=mas 

00011630 

AHIGHsALOW 

nn*?!  iR4n 

KAPAX=PAS 

00011650 

10 

IF<eETA2("AS>.GE.8L0U>  GO  TO  20 

O0C1166O 

Ri  nu=RrTA;(NA<;> 

nnn i i cin 

K6P1N=NAS 

00011680 

20 

IF(eETA2(  NAS)  .LE.BHIGH)  GO  TO  30 

00011690 

nHir.MsRrTA^tM  A^i 

nnn 1 i 7n  n 

KeNAXsNAS 

00011710 

JO 

1F(ALPHA2(NAS).GE.AL0U)  go  to  AO 

00011720 

41  ni4S4i  PM49fM4<;\ 

nnni i 7xn 

KANINsHAS 

000117AO 

*0 

if(alpha2(nas>.le.ahigh)  go  to  sc 

00011750 

AHTrM=ALPHA?C  PAS) 

nnni iTcn 

X APAX=PAS 

00011770 

50 

CONTINUE 

000117R0 

IF<NP.6T.01  GO  TO  55 

ftfini 1 7Pn 

aavg=o. 

00011800 

BAVG=0. 

POOIIRIC 

no  1 =1  .fciPK 

nnni \»on 

« = NAS-L*1 

0001 1830 

BAV6=BAVG*eETA2IPI 

0Q0118A0 

Si- 

AAWC=A*WC«A1  PMtailtl 

IU)0.t.l«.50 

AAVCsAAVG/NPTS 


00011860 


G LtVEL  21 


niNHAX 


DATE  > 78103 


08/SO/44 


8AVGSB4VC/NPTS 

00  54  L=l*4 
tPI  Tfll  tl  AS*(  izl 

00011870 
00011880 
aoBi i8«a 

APLOKLAS  *L)  = AAV6 

00011900 

54 

BPL07(LAS«L)=BAVG 

00011910 

fen  rn  fea 

55 

C0N71NUE 

00011930 

DO  56  LSI  *4 

00011940 

APLOT(LAS*1>sAHIGH«0.25 

00011960 

ePL0T(LAS*l)=BE7A2<KAWAXl 

00011970 

Rp^  nT#i  ru-n^?^ 

nnni  i«ar 

APL0T(LAS*2):ALPHA2(KPPINI 

00011990 

A plot  <LAS*3>=AL0W-C.25 

00012000 

Anni 9ni a 

PPLOT<LAS*4)rBhICH*0.25 

00012020 

»PL0T1LAS»4)  = ALPHA?*KB«AX  > 

00012030 

•■ft 

^ Afesi  A 

Aoo 1 9rftn 

bO 

CONTINUE 

OOC120SO 

IST=NPSI4*4 

00012080 

70 

continue 

00012100 

RETURN 

00012110 

G LEVEL  21 


PROBL 


DATE  = 78U3 


OS/5A/AA 


SUBROUTINE  PR OBL ( OATACtOA Tt S,HS «APLOT t BPLOT«XPLOT ,YPLOTf OUORK S*  00012130 
1DUPSP2*VPR0BL«IPLTQL>  000121*0 
DIPEW&IOW  DATACtH.DATASHI.HSll>«APLQTtll.ePLOTtl>.KPLQT<ll.  00012150 

1 VPLOT< 1)«OUORKS< 1>«DUKSP2I 1> tIPLTOLI 1)  00012160 
OIPENSION  RPS(200>,RPSN(200>«SPS(200>  00012170 
n I PENSION  BPSTiiiPt-Sniii.MTMPtjawBP«:niiTi*nni.MtmTt?in.i>p<;uKiiBai>i.  nfcigmn 


1 NUK1<201*RPSUK2(50C>»NUK2(20)«RPSWK3(500>  *NyK3(2  0> 
AL0U30ATAC(16> 

AMICMSPATACH71 


PETHrOATACi 19) 

PETL=DATAC(le) 

-J  A REAstAHlGM-ALQm.tPPTH-BrTLI 


NPANELsOATACI 23)*  0.5 
NPSI=0ATAC«9) *0.S 


00012190 

00012200 

00012210 


00012220 

00C12230 

000122*0 


00012250 

00012260 


Mf  liMjsr 

nnn i 997n 

NsuRri=o 

SWEEPsO. 

— LASPafl- 

00012260 

000122*0 

LASsO 

00  16  I:l*NPANEL 

NPSIP-NCST 

00012310 
00012320 
ft  nn  1 9SS  (1 

NPSISsNPSI ♦! 

DO  P J=1.NPS] 

IF  ( tPL  TCL  (L*SP*H.f  0.0.  ANfl.WPStS.r.T-NPSI)  NPSI  SrJ 

000123*0 

00012350 

nnm 

IFINPSIS.CT.NPSI  ) GO  TO  » 

IP(NPSIP.LT.NPSI)  GO  TO  P 

IF  (IPt  IDL  (L  ASP«1  I.FC.I  1 MPSIP=.l-t 

00012370 

09912360 
nnm  sspn 

>* 

LASP  = LASP»* 

00  12  J=1»NPSI 

- ()0  TO  la 

00012*00 

00912*10 

IFIJ.GT.NPSIP 1 GO  TO  12 

IFKNPSIP'NPSISI.EG.C)  00  TO  9 

IP  f.i.cT.uoctti  tn  Tn  1 

00012*30 

09912**0 

BAVGl>nPLOT(L AS** )*PPL0T(LAS*2) 

BAVC2sPPLOT(LAS*P  >*I>PL0T(LAS*6I 

IPIOTIDzAPlCTd  IS**) 

00012*60 
00012*70 
non  1 9*«n 

a. 

tPL0T(l)s»PL0T(L*S**) 

IF(PAVGl.CE.eAVG2)  SO  TO  10 

■ -XPLOTtIl»BPLOTtLAS*2l 

00012**0 

00012500 

00012510 

ID 

VPLOTd  >sAPLOT(L*S*2) 

LPxJ-NPStS*2 

XPICTK-PISBPIQTII  *S«)) 

09912520 
00012530 
ftnni 95*n 

TPL0T(LP)**PL0T(L*P*1) 

IF(  IMPSIP-NPSISl.EQ.O)  GO  TC  ICO 

T.MPSfPi  cn  Tff  n 

00012550 
00012560 
ftftn  t ?*Tn 

fc  AVCl  sPPLOHLAS**)*PPLOT  (LAS *2) 

BAvC2zPPL0T(LAS>  *BPLCT(LAS-2> 

L60 Pt  » -t.at  ut  ii  

00012f.Ra 

99P12595 

TPLOTaP*l)*APLOT(LAS**) 

IP((NPSIP-NPStS>.Ca.O>  60  TO  11 

|P4>*uc.i.er.*Auc*i  en  rn  ii 

00012610 
00012620 
ft  nil  1 9*  sn 

»PL0T(LP*l)a“PL0T<LAS*2) 

009126*0 

194 


G LEVEL  21 


PROBL 


DATE  = 78103 


08/54/4R 


YPL0T<LP*1)=*PL0HL*S*2) 

00012650  1 

1 

11 

LPs2*(NPSlP-flPSIS^l)^3-<J-NPSlS4ll 

VPI  aTfLPlsRPLOTfL 

00012660  1 

nftAi9fi7n  1 

1 

TPL0T(LP)=»PL0T(L*S*3) 

000126RO  1 

- 

12 

L*S=L*S** 

IFINPSIS.KT.NPSTl  GO  TO  1R 

00012690  1 

aaoiPTaa  I 

1 

NCFS=2«<NPSIP-NPSIS*l)«-2 

00012710  1 

1 

1 

• 

CALL  SORT(YPLOT,XPLOT,NCFSlt3ETH.PETL.OWKSP2) 

00012760  i 

lF(NCFSl.tQ*0)  GO  TC  14 

rAii  ^ABTrvDi  rT-.vDinT-.wrP«i-AuTi:u^iinu.nL‘i^cDo> 

00012770  } 

00012*^^0  1 

1 

IFlNCFSl.EQ.O)  GO  TO  1* 

00012790  ‘1 

CALL  AREAlXPLOTtYPLOT^NCFSlfYINTG) 

^urFP=4urrp*YTh.Tr. 

30012600  -1 

R0012030  \ 

NSURFI=NSURF1*1 

000128*0  1 

NINP(NSURFl)sNCFSl 

NPr9#NrF^1 

00012650  t 

nnniPA^n  | 

DO  13  Lsl«NCFSl 

00012870  1 

J2atL-l)«2*l 

RPSINP(NCUPI»J2)=XPLOTCL> 

00012880  1 
0001 2a90 

13 

RPSINP (NCUMI*J2*1 )=YPLOTIL> 

00012*00 

1 4 

NCUHIsNCUPI^N2 
f nwiTwiir 

00012910 

ncoipqpn 

00  15  l<KKsl.NCFS2 

00012930  i 

LPL=(RKK-1>«2 

000129*0  '■ 

vpf  (tTciririri  -nunpi/di  mi  *i  i 

^0012950 

15 

YPL0T<KFK>=0W0RKS<LPL*2) 

00012*60  J 

1 L 

NCF$2NCFS2 

r Ai  1 T f YPi  nr  .YPi  r.T.Mrp  .pr  ti  .ruk^p  ? » 

00012970  * 

00012940 

IFINCFS.EQ.O  GO  TO  18 

00012990 

CALL  SORTlXPLOT«YPLOT«NCFSf AHl6HtALCU«DWKSP2> 
iFfMrF.4.rn.ni  fin  in  ia 

00013000 
nflni  :^ni  n 

CALL  AREA(XPL0T,YPL0T,NCFS»YINTC-  ) 

00013020 

SWEEP  sSWCtP^YINTG 

MSURF  IsAiSlJRFT^l 

00013050 

0001 30Aa 

NINP(NSURFI)=NCFS 

00013070  ) 

N2=2oNCFS 

00013080 

-L2- 


J2*CL-1>*2»1 

RPSINPCNC UMI*j2)s*PLOT«H 
-aPSINP<WrUHl*.12»l  ISYPLtlTf 


J-X. 


C0013100 

0C013110 

000X8120 


NCUHIsNCU8I*V2 
18  CONTINUE 
VPR08L«0. 


00013130 

000131*0 

30813150 


I FINSURFI  .EQ.O  RETURN 
NPTSsNINPIl) 

N?a?«WPTS 


000131FO 
00013230 
and  1 30* a 


00  30  Lsl«N2 
30  RPSOUT(L>sRPSINP(L> 

MnumuhPTs 


00013250 

0001328C 

aaaiRpTo 


NSURFOsl 


00013280 


e LCVCL  21 


PROBL 


0«TC  s 78103 


08/54/44 


I 

\ 


I 


i 


NCUn0sN2 

IFlNSURFI.CO.l)  GO  TO  85 

MCHMISM? 

DO  Cu  Is2,NSURFI 
NR=NINP(I> 



00  40  L=1»M2 

40  RPSUKl(L)=RPSlNP<L*NCUni> 

NUKlCHsMB 

NSUK1=1 

NCUni=NCUHI*N2 

tLSUHsJ 

DO  70  K=ltNSURFO 
NS=NOUT«K> 

4?s2.WS 

DO  50  L=1«R2 

50  SPSIL)=RPS0UT(NSUM*L) 

MSUHsMSUK*K3 

NCUP1=0 

NCUH3=0 

MSUH3=0 

DO  62  JslfNSUKl 
NR=KIUK1(J) 

|[l2-2«tUI 

00  55  L=1  .N2 

55  RPS(L)=RPSyKl(L*NCUPl> 

HCUMl  sMCUai.*W2 

CALL  0VRLAP(RPS«rjRtSPS«NStRPSWK2«NSWK2tNUK2«f(PSN) 
lP(NSWK2.Ca.O)  GO  TO  83 

DO  550  N1s1»JlSUK2 

NPTS=NUK2<  Ml ) 

1P1NPTS.6T.5C)  GO  TO  80 

iS4j C0IILH^^IE 

NCUH2=C 

00  6D  NI=ltNSUK2 

H4U*3sMSy>C3*-l 

NPTS=NUK2<NI) 

MyK3«NSMK3>=NPTS 

M2S2.MPTS 

00  56  L=1«N2 

56  RPSyR5lNCUWI*L)=RPSWK2(NCUH2*L> 

MCUP.B  aNCUr2««2 

NCUH3sNCU«3*N2 

60  continue 

.62 CONTlNUr 

NSyKI :NSyK3 
VCU«1=0 

oa-6i-  oaraisuKj 

NPTS=NUK3(J) 

N2=2«NPTS 

NyiftlJlsWPTS 

DO  65  L=ltN2 


00013290 
00013300 
QOni  3310 
00013320 
00013330 
000«3A0 
000133S0 
0001 3360 
anni  .3.370 

00013380 
00013390 
00013400 
0001 3410 
00C1342C 

j.ccLa4?? 

00013440 
00013450 
0C0t3460 
00013470 
00013480 
00013490 
00013500 
0001  3510 
00013520 
00013530 
00013540 
<1001.357  0 

0001 3610 
00013620 
002.1263.11 
00013640 
rool 3650 
nnniJ660 
00013670 
00013680 
000.12690 
90013700 
0001 3710 
00012720 
00013730 
00013740 
00011750 
0001 3760 
00013770 
00013780 
00013790 
3001 3800 
00013B1C 
00013820 
00013630 
00013840 

00013850 


DATE  = 78103 


65  RPSWKllNCUHl«L>~RPSUK3(NCUni*L> 
68  NCUnisNCUni*N2 


4*1. 


NCUHlsO 

00  J=1«NSUK1 


NOUT(NSUftFO»J)=NPTS 

M2=2*NPTS 


72  RPSCUT(NSUM*L  >=RPSUK1(NCUNI*L) 
NCUPlrNCUWl»N2 


74  CONTINUE 

NSURF0=NSURF0«NSUK1 


08/54/44 


00013860 

00013870 


00913890 

00013900 


00013920 

00013930 


00013950 

00013960 

44ni  4a7n 


00013980 

00013990 


00  95  IsltNSORFO 
NPTS=N0UT« I) 


00  93  L=lfNPTS 
L2=»L-1)»2»1 


<»3  YPL0T«L)=RPS0UT<NCUP0*L2*1  > 

CALL  4REAcXPL0T.tPL0T»NPTS«YI NTG) 

suAiiSiiAitrum 

95  NCUMO=NCUrO*N2 

VPR0FL  = SWA2/TAREA 


ifivprobl.le.i.j  return 

WRITE  I6«200l 


1«  OVERLAP  PROBABILITY  GREATER  THAN  ONC«  ONE  ASSUNEO*) 
VPROBL=l. 


00014160 

00014170 


00014200 

00014210 


00014230 

00»’14250 

C0016280 

00014290 

00014330 


90014341 

00014342 


IllUift  ■ 


00014344 

00014345 


S LEVEL  21 


OVRLAP 


DATE  s TS102 


10/47/46 


1 


I 

I 

I 


r 


r 

1 

j 

[f 


! 

j 


SUWOUTINC  0«KLAPCRPS«NR*SPS«NS*I)PS1N*MRPLSN«NRPTSN*RPSN> 
OINCNSION  RPS(1>*SPSC1I«RPSN(1I«RPSIN<1)*NRP7SN<1> 

CALL  CWCCK4tRPS«NR«SPS«lt8l 

0001 4370 
00014380 
00014390 

\ 

i 

iS 

NRPLSNsO 

00014400 

i 

NSUHsO 

00014410 

j 

tRxl 

00014420 

• 

1R231 

00014430 

i 

16 

XCsRPS(lR2> 

00014440 

j 

ycsRPsriR2«i) 

00014450 

1 

NS2«2*NS 

00014460 

1 

CALL  ClRCLE(SPS«NS«XC*TCtINCR> 

00014470 

IFC1NCR.NC.2)  60  TO  IT 

00014480 

j 

IFdR.Ca.NRI  60  TO  120 

00014490 

IRsIR*l 

00014500 

i 

IR2«2«1R-1 

00014510 

60  TO  16 

00014520 

! 

17 

IPLA6SS 

00014530 

i 

IFIlNCR.EQ.ll  IFLA6S6 

0001454  0 

i 

NEXTsO 

00014550 

! 

IFIIFLA6.Ea.6>  NRNsO 

00014560 

IFIIFLAI.CQ.A)  60  TO  25 

10014570 

1 

NRPLSNxNRPLSN*! 

00014580 

3 

NRNxl 

00014590 

iipsaif  1 istc 

00014600 

RPSNI2)sTC 

00014610 

1 

IFLA6S1 

00014620 

1 

7% 

IRPlstR«l 

00014630 

1 

IFdR.EQ.NRI  IRPlsl 

00014640 

] 

IR2Pls2»IRPl-l 

00014650 

■ 

tOOUBLsO 

00014660 

30 

X1=RPS(  IR2I 

00014670 

i 

VlxRPS<IR2*l> 

00014680 

35 

X2>RPSdR2Pl) 

00014690 

J 

T2zRPSdR2Pl*l> 

00014700 

CALL  INTCRC(SPStNS*XltYl,X2»T2, 

X, V.KLtIFAULTtlOOUBLt 

00014710 

1 XALSOtTALSO) 

00014720 

IFdFLAC.EO.A.ANO.IFAULT.Ca.O) 

60  TO  70 

00014730 

IFdFLA6.C0.1.ANO.IFAULT.eO.O> 

60  TO  40 

00014740 

rF<tFL16.Ca.4.1NO.IFAULT.EO.Q) 

60  TO  40 

fi&OI 47^n 

IF(IFLA6.Ce.6>  GO  TO  38 

00014760 

? 

CALL  CHECKl(X2*Y2*RPSNtMRN«ICNK> 

00014770 

TFIirHK.ro. 11  GO  TO  90 

000147B0 

j 

XFlAGxI 

00014790 

NRN«NRN*1 

00014800 

NRN2S2*NRN-1 

00014810 

J 

RPSN(NRN2>sX2 

00014820 

! 

RPSN(NRN2*l)sV2 

00014830 

31 

IRxIRPl 

OOC14840 

IR2s2*IR-l 

00014850 

IFdRPl.E0.1.AN0.1FLA6.C0.6>  60  TO  90 

00014860 

i 

ltPl«lR«l 

nnnt4A7n 

IFdR.EO.NR)  IRPlsl 

00014880 

1 

■ 

198 
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IR2Pls2*IRPl-l 

00010890 

IDOU8L=0 

00010900 

00 

IFUCS2 

00010920 

NCXTslRPl 

00010930 

XSEPsX 

00010900 

tSEPst 

00010950 

CALL  CHECK1(X»Y*RPSN*NRN*ICHK> 

00010960 

IFllCHK.EO.l > 60  TO  90 

00010970 

NRN=NRN«1 

00010980 

NRN2s2*NRN-l 

00010990 

RPSN(NRN2>:X 

00015000 

RPSN(NRN2*1>=Y 

00015010 

IS=KL 

00015020 

IS2=2»IS-1 

00015030 

ISP1=IS*1 

OOOISOOO 

IF(IS.CQ.NS)  tSPl=l 

00015050 

IS2P1=2*ISP1-1 

00015060 

iooueL=i 

00015070 

XALSOsX 

00015080 

YALS0=Y 

00015090 

KDIR=0 

00015100 

so 

X1=SPS(IS2) 

00015110 

YlsSPS<IS2*l) 

00015120 

X2:SPS(IS2P1) 

00015130 

Y2=SPS(  IS2P1*1> 

00015100 

CALL  INTCRCIRPSfNR,Xl«Yl*X2tY2t 

X«Y*KL«IFAULT«tOOUBL«XALSO«TALSO) 

00015150 

IF(IFLAG.CQ.2.ANO.IFAULT.CQ.O) 

GO  TO  75 

00015160 

IFCIFLAG.EO.S.ANO.IFAULT.EO.O) 

GO  TO  75 

00015170 

IFIIFLAG. £0.3. ANO.KDtR.CO.il  60 

TO  60 

00015180 

lF(lFLAG.Ea.3)  GO  TO  56 

00015190 

IFdFLAG.EO.T.AND.IFAULT.CO.O) 

GO  TO  75 

00015200 

IFClFLAG.CO.T.ANO.KDIR.EQ.l)  60 

TO  50 

00015210 

CALL  ClRCLC(RPS.NR«Xl«Yl«tNCR I 

00015220 

IF(INCR.CQ.l)  GO  TO  56 

00015230 

IFCtFLAG.EQ.ll  CO  TO  30 

00015200 

50 

CALL  CtRCLE(RPStNR«X2,Y2.INCR> 

00015250 

IFdNCR.CQ.ll  GO  TO  60 

00015260 

IFIIFLAG.EQ.7)  IFLA6=1 

00015270 

IFdFLAG.EO.l)  GO  TO  30 

00015280 

IFdFLAG.CO.T)  IFLAGsl 

00015290 

IFLA6=7 

00015300 

tOOUBLsl 

00015310 

XALSOsX 

00015320 

YALSOsY 

00015330 

IFCINCR.Ca.O)  60  TO  58 

00015300 

GO  TO  68 

00015350 

58 

NRNsNRN*! 

00015360 

NRN2s2»NRN-l 

00015370 

RPSN(NRN2)3X1 

00015380 

RPSNINRNX^llsTl 

00015390 

tOOUBLsO 

00015000 

\ 


tiil-J-.'/V. 
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1FLA6S3 

08015410 

58 

ISPlslS 

00015420 

IS=1S-1 

00015440 

IFlIS.EQ.OI  ISsNS 

00015450 

KO  IR*-1 

^ 00015470 

GO  TO  50 

00015480 

60 

NRN2S2.NRN-1 

00015500 

RPSN(NRN2>sX2 

00015510 

RPSN(NRN2*1>=Y2 

30015520 

lOOUBLsO 

00015530 

IFLAGS3 

00015540 

68 

IS2=2*IS-1 

00015560 

1SP1*IS*1 

00015570 

IS2P1=2»ISP1-1 

00015590 

K0IR=1 

00015600 

GO  TO  50  00015610 


70  IFCNRPLSN.EQ.O)  GO  TO  74 

CALL  CHCCK2(XtY«RPSlNtNRPLSN«NRPTSN*ICHK2> 
lFHCHK2.EQ.il  GO  TO  120 


74 


76 


NRPLSN=NRPLSN«1 
GO  TO  76 


IM.lCMKl 


IFdCHK.EO.  1)  GO  TO  90 
NRNzNRN«l 

HRN2  = 2»NRN-1 


RPSNINRN2>=X 

RPSNCN8«l2*l)rr 

JLUlt 


YlsT 

lOOUBLsl 

XALSQsX 


YALSOsY 

IF(IFLAG.CQ.6I  IFLAG=1 
IFdFLAG.EQ.ll  GO  TO  35 


IRP1=KL*1 

IFtKL.EQ.NRl  IRP1=1 
IR2P1=2«IRP1-1 


_2JL 


1FLAG=4 
GO  TO  35 

IFtHRW.EQ.O)  GO  TO  Ufi. 


IF(NRN.LE.2>  NRPLSNsNRPLSN- 
IF<NRN.LC.2I  GO  TO  120 
IFtNEXT.EQ.O)  GO  TO  100 


NRPTSNCNRPLSN>=NRN 

N2s2*MRN 

DO  95  1=1.1112 


95  RPSINd«NSUR)=RPSN(I> 


00015620 
00015630 
■ 00015640 


00015650 

00015660 

nOfllKATO 


00015680 

00015690 

00015700 


00015710 

00015720 

noaisTao 


00015740 

00015750 

00015760 


00015770 

00015780 

00015790 


00015800 

00015810 

00015820 


00015830 

00015840 


00015860 

00015870 

00015880 


00015890 

00015900 

00015910 


00015920 
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NSUnsNSUn*N2 

00015930 

NRNsO 

0001 S9A0 

XlsXSEP 

Y1=YSEP 

00015960 

IFLAG=6 

00015970 

IR2P1=2«IRP1-1 

00015990 

lOOUBLsl 

00016000 

YALSO=YSEP 

00016020 

IFINRPLSN.EQ.5)  GO  TO  120 

00016030 

GO  TO  35 

00C160A0 

100 

X1=SPSI 11 

00016050 

Y1=SPS<2) 

00016060 

CALL  CIRCLE(RPSN,NRN.Xl>YltINCR) 

00016070 

IF< INCR.EQ.O)  GO  TO  110 

00016080 

X 1=RPSN(NRN2> 

00016090 

Y1=RPSN<NRN2*1) 

00016100 

CALL  CL0SC(SPS«NS<X1«Y1  «KN> 

00016110 

X1=RPSN(1) 

00016120 

Yl=RPSNf2> 

00016130 

CALL  CLOSECSPSfNStXl.YltKl) 

00016190 

00  103  K0IR=1«2 

00016150 

K=KM 

IFACT=<-1)*»»K0IR*1» 

00016170 

DO  102  IsltNS 

00016180 

L=K*IFACT»<I-1) 

00016190 

IFIL.GT.NS)  L=1 

00016200 

if<l-lt.i>  l=ns 

00016210 

IF(L.EQ.1.AN0.KDIR.CQ.1>  K=2-I 

00016220 

IF(L.E0.NS.AN0.KDIR.ea«2>  KsN 5*1-1 

00016230 

L2=2*L-1 

00016290 

«=NRN*I 

N2S2.P-1 

00016260 

RPSN(M2>=SPSCL2> 

000 16270 

RPSNfn2*l)=SPS(L2*l) 

noRi  f^9r^a 

IFfL.EO.Kl)  GO  TO  1020 

00016290 

102 

CONTINUE 

00016300 

KK2=2*KK-1 

00016320 

X1=SPS(KK2> 

00016330 

Y 1=SPS<KK2*1> 

00016390 

CALL  CIRCLE(RPSNtNtXl«Yl*tNCR> 

00016350 

IF( INCR.EQ.I)  GO  TO  103 

00016360 

102S 

CONTINUE 

OC016370 

GO  TO  lOA 

00016380 

103 

CONTINUE 

00016390 

WRITE (6<200> 

200 

FORNaT<1H1,*SOPCTHING  wrong  OVRLAP  SUBROUTINE'I 

00016910 

STOP 

0001692C 

NRPTSNINRPLSNlxN 

N2x2*n 

00016990 
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00  1040  lsl*N2 

000164SO 

1Q4« 

RPSIN(I*NSUN>=RPSN(1) 

OOOIGAGO 

OOOIGATO 

RPSN(l)sflPS(NRN2> 

00016480 

RPSN(2)=RPS(NRN2*1) 

00016490 

K=KN 

00  107  l3l«NS 

00016510 

L=K-IPACT*« I-l) 

00016520 

IFIL.ST.NTJ  LSI 

IFlL.LT.ll  L=NS 

00016540 

tF(L.EQ.l.ANO.IFACT.EQ.-l>  K=2-I 

00016550 

L2=2«L-1 

00016570 

M=I*1 

00016580 

RPSN(n2>=SPS(L2> 

00016600 

RPSN<P2M)=SPSCL2*H 

00016610 

TFlL.EQ.Rtl>  SO  TO  108 

00016620 

107 

CONTINUE  ■ 

00016630 

108 

NaN*l 

00016640 

N2=2*n-1 

RPSN(N2>=RPS(1> 

00016660 

RPSNTN2*1»=RPSC2) 

00016670 

NRPLSN=NRPLSN«1 

00016680 

GO  TO  115 

00016690 

no 

NsNRN 

00016700 

IIS 

NRPTSN(NRPLSN>=N 

00016710 

N2s2«M 

00016720 

00  no  Isl«N2 

00016730 

1 1 A 

RP<;lN  1 I*NSUH)=RPSNf  I > 

00016740 

NSUNsNSUn*N2 

00016750 

120 

CONTINUE 

00016760 

RETURN 

00016770 

ENO 

00016780 

6 LEVEL 
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SUBROUTINE  I NTERC ( POINTS*NPTS tX 1* Y1 *X2« Y2t X « Y «KL«I FAULT« 

00016790 

1 100UBL«XALS0«Y*LS0) 

00016800 

■■iHvnTrFm 

OIRENSION  IPTS<200),XINT(200) ,YINT(200) 

00016820 

C0RM0N/T0LRS/T0LEQN,T0LSNT*T0LQPT*T0L«T0LCKlfT0LCK2» 

00016830 

1 T0LCK4.DELT.T0LCIR .ANGST2.T0LPLT 

NINTC=0 

000168S0 

IFAULT=1 

00016860 

DO  7A  Ixl.NPTS 

00016870 

IP1=I*1 

00016880 

IFII.EQ.NPTS)  IPlrl 

00016890 

I2=2«I-1 

00016900 

I2P1=2*IP1-1 

00016910 

X3=POINTS(I2) 

00016920 

Y3=P0INTS«I2*1) 

X4=P0INTS(I2P1> 

00016940 

YA=P0INTS(I2P1*1> 

00016950 

XKAX=X3 

00016970 

YMIN=Y3 

00016980 

YMAXSY3 

00016990 

IF(X*.LT.X3)  xpin=xa 

00017000 

IF(X4.GT.X3)  XMAX=X« 

00017010 

IFCY4.lt. Y3>  YMIN=Y4 

IF<Y4.GT.Y3)  YMAX=Y4 

00017030 

XMIN=XMIN-TOL/2. 

00017040 

XMAX=XR4X»T0L/2. 

000  17060 

TMAX=YRAX»T0L/2 . 

00017070 

nri  1 

0EL2=X4-X3 

00017090 

IFCABSCOELD.LT.TOL)  go  TO  10 

00017100 

V9-Y1  l/nri  1 

non  171  IQ 

SB1=Y1-SA1*X1 

000  17120 

10 

IF(ABS(0EL2).LT.T0L)  GO  TO  20 

00017130 

SA2  = (Y4-Y3)/0EL2 

0001  7140 

SB2=Y3-SA2»X3 

00017150 

20 

IFCABSCOELD.GE.TOL)  GO  TO  40 

00017160 

ifcabsidcl21.6E.tol)  go  to  30 

00017170 

GO  TO  74 

0C017180 

30 

N1NTC=N1NTC»1 

0001 7190 

XINT(NINTC>=X1 

0001  72CO 

YINTCNIMTC>=SA2*X 1+SB2 

00017210 

IPTSCNINTC)=I 

00017220 

GO  TO  70 

00017230 

XO 

IFCABSCOCL2).GC.TOL)  GO  TO  50 

00017240 

NINTC=NINTC*1 

0001 7250 

XTNTlNTNTr)sX3 

0001 7260 

YINTCNINTC)=S  A1*X3»SB1 

00017270 

IPTSCN1NTC)=I 

00017280 

SO  TO  70 

00017290 

50  01FFkxS*l>SA2 


OOOlTJOO 


G LEVEL 


INTERC 


DATE 


78102 


10/47/48 


21 


IF(ABS(D1FFA>.LT.T0L)  go  to  74 

NIMTC=NINTC*1 

XINT<NINTC)=<SB2-SB1)/DIFF4 

00017310 

00017320 

00017330 

TINT  CNINTC)=SA1»XINT1MINTC)*SB1 

00017340 

IPTS<  NINTC)  = I 

00017350 

70 

CONTINUE 

00017360 

CALL  CHECKS!  X«IN»XM  AX  ♦TMINtTN  AX  tX  INTI  HI  NTC)»T  INTI  NINTC»»IBETU) 

00017370 

IFIIBET  J.EQ.n  IPTSINIMTOsO 

00017380 

74 

CONTINUE 

00017390 

IF ININTC.GT.O > IFAULTsO 

00017400 

IFI iFAULT.EQ. 1»  RETURN 

00017410 

XMIN=X 1 

000  1 7420 

XXAXsXl 

0001 7430 

YMIN=Y1 

00017440 

YMfcX=Yl 

IF(X2.LT.X1>  XMIN=X2 

0001  7460 

IFIX2.GT.X1)  XPAX=X2 

00017470 

IFIT2.LT.T1)  TMIN=V2 

00017480 

IFIT2.GT.T1)  TNAX=T2 

00017490 

XMlN=XMIN-T0L/2. 

00017500 

THIN=TN1N-T0L/2. 

00017510 

XNAX=XMAX«T0L/2. 

00017520 

TMAX=TMAX.T0L/2. 

00017530 

DO  80  Isl.NINTC 

0CC17540 

IFIIDOUBL.EQ.O)  GO  TO  75 

00017550 

IFIABSIXINTI I)-XALS01 .LT. TOL. AND . ABS I T INTI  I)-TALSO) .LT. 

000  17560 

t TOL)  IPTSII)=0 

0 00 1 7570 

75 

IF(XINTII).LT.XMIN)  IPTS(I)=0 

OOC175RO 

IFIXINTID.GT.XMAX)  IPTSU)=0 

00017590 

IFITINTII ).LT.TMIN)  IPTSII)=C 

0001 7600 

IFITINTlI ) .GT.THAX ) IPTSII)=0 

0001  7610 

80 

CONTINUE 

00017620 

IFAULT=1 

00017630 

DO  90  I=1.NINTC 

000  17640 

IFIIPTSID.EQ.O)  GO  TO  90 

00017650 

IFIIFAULT.EQ.O)  GO  TO  85 

IfAULTsO 

00017670 

DIST=IXINTII)-X1>..2*ITINTII)-T1)  **2 

00017680 

N = I 

00017690 

GO  TO  90 

00017700 

85 

DIST2  = IXINTI I )-Xl )•• 2*1  TINT  I I )-Tl )..2 

00017710 

IFI0IST2.lt. DIST)  N = I 

00017720 

IFIN.EQ.I)  OIST=OIST2 

00017730 

90 

continue 

00017740 

IFIIFAULT.EQ.1)  RETURN 

00017750 

KL=IPTSIH) 

00017760 

X=X1NTIN) 

00017770 

TSTINTIM) 

HTITIZlll 

RETURN 

00017790 

END 

00017800 
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subroutine  CLOSE! ARRAr*NPTS»Xl*YlfK) 

DINENSION  ARRAVdl 

K=1 

00017810 

00017820 

00017830 

01ST=(ARRAY!1>-X1)»»2*I  ARRAY<2)-Yl>*o2 

00017840 

IF(NPTS.EQ.l)  RETURN 

00017850 

00  10  1S2.NPTS 

00017860 

I2s2»I-1 

0 00  1 78  70 

01ST2=!ARRAYII2)-X1)»»2*!ARRAYII2*1>-Y1>»»2 

00017880 

IF(0IST2.LT.0IST)  K=I 

— — 1 IlillJI 

10  CONTINUE 

00017900 

return 

00017910 

END  OC017920 


205 


6 LCVCL  21 


CHCCKl 


o«re  s TBili 


ia/«7/«« 


SUBROUTINC  CMeCKl(Xl*Tl*«RK*V*N»ICORP> 
OINCNSION  *RR«T(1> 


1 T0LCK««0eLT«T0LCIR*ANCST2«T0LPLT 
ICORPsO 
DO  IB  Is 


I2s2*I-l 
XsARRAY(I2l 
sARR 


IFIABSfX-Xlt.LT.TOL.ANO.ABStY-Yll.LT.TOL)  60  TO  15 
10  CONTINUC 

RETURN 

15  ICORPsl 

omiB  H 


BBB1T9SB 

0BB179«t 


00017BA0 

OeOlTBTO 


I I 


aoaiTBM 

oeeiaooo 


0001S020 
00010030 
ooomoAO 
00010050 
BOO  1 BD6a 
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SUBROUTINE  CHECK2IX 1«V1 »MRAV «NPLS«NPTAR« ICONP > 

OlHENSION  ARR8V(1)*NPT«R(1> 

CONNQN/TOLRS/TOLEQNiTOLSNT.TOLOPT.TOLINT.TOLCKl.TOL. 

00018000 

08018100 

1 T0LCK4«0ELTtT0LCIR*ANCST2«T0LPLT 

00010110 

ICONRxO 

00018120 

NSUNsO 

00018130 

00  20  LslfNPLS 

00018140 

NPTSzNPTARCL) 

00018150 

00  10  I=1»MPTS 

000181G0 

K2=NSUn*2>l-l 

00018170 

X=ARRAY(K2) 

00018180 

T=ARRAY(K2*1> 

00018190 

IF(ABS(X1>X>.GE.TQLI  GO  TO  10 

00018200 

IF(A8S<Y1-Y).GE.T0L)  go  to  10 

00018210 

ICOHPsl 

00018220 

RETURN 

00018230 

10 

CONTINUE 

00018240 

NSUN=NSUH*2*NPTS 

20 

CONTINUE 

00018260 

RETURN 

00018270 
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CHECKS 
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iU7usa 

IF(X.L7.MH1N>  tBC7U3l 

lETUl  08018290 

00010300 

00018310 

IFtX.CT.XHAXt 

IBETysl 

00018320 

IFIT.LT.TNINI 

IBE TU31 

00010330 

IBETHsl 

RCTUMN 

00010390 

END 

00010380 
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CHCCK4 


OATC  « T8102 


10/%T/«8 


SU8R0UT1NC  CNCCK4CRPS«NR«$PS«NS> 

OINCNSION  RPS(l>t$PS«ll 

COWtOM/TOLRS/TOLEaW.TM.S>lT.TOLQPT.TOLINT»TOLCKl.TOLCK2. 
1 T0Lt0ELT«T0LClRtRN6ST2«T0LPLT 
DO  30  lsl«NR 

I2g<I-n«2«l 

12P2sI2*2 

IFCI.CO.NRI  I2P2S1 

XRl»RPSa2> 

TR1<RPS(I2*1> 

XR2sRPS(I2P2> 

TR2sRPSCl2P2*I) 

ISUTCHsi 

IF(«8SIXR2-XR1).LC.T0L)  ISUTCH=2 

IF<«BS(TR2-TR1).LE,T0L>  ISMTCMS3 

IFdSUTCH.NC.l)  GO  TO  IS 
*SL0PErfTR2-YRl>/(XR2-XRl> 

BlWTsYRi-«SLOPC«XRl 

AINSLs-1. /ASLOPE 
IS  00  29  Jzl*NS 

J2»<J-I>«2*1 

X$1=SPSCJ2> 

Y$l3SPS(J2*l) 

to  TO  (U.IT.IB),  ISWTCH 

Ifc  XS1L*«YR1-YS1*A1HSL»XS1-ASL0P£»XR1J/IAINSL-ASL0PE> 

TS ILSTRI* ASLOPE* (XSIL-XRII 

DlST«SQRTtlXSlL-XSl)««2*<TSlL-TSl>««2l 

IF(OIST.GT.TOL)  GO  TO  29 
ALPHAs«TAN(AINSL> 

XS 1 P»X S IL* DEL T»COS< ALPHA} 

TSIPsTStL*OELT*SlN(  ALPHA) 

CALL  CIRCLEIRPS.NRtXSlP.TSlPtlNCR) 

IFdWCR.EB.O)  to  TO  ItO  

XSlPsXSlL-OELT (COS (ALPHA) 

TS1P«YS1L*0ELT (SIN (ALPHA) 

lAO  SPS(J2)«XS1P 

SPS(j2«t) XYSIP 
to  TO  20 

17  1F(ABS(XS1-XR2I.6T.T0L)  tO  TO  20 

XSlPsXSl*OELT 

rSlPxYSl 

CALL  CIRCLC(RPS.HR.XSIP.TSIP.INCR) 

tF(INCR.CO.O)  CO  TO  lAO 
XStPxxSl-OELT 

to  TO  ItO 

18  IF(ABS(YS1-TR2I.6T.T0L>  CO  TO  20 
YS1P*YS1»0ELT 

AtlFxKRl 

CALL  CIRCLE (RPStNRtX SIP tTS IP* INCR) 

IF(INCR.Ea.O)  GO  TO  IGO 

YtlP«T$l»OCLT 

to  TO  IGO 


00018370 
00010380 
00818390 
000  1 84  0 0 
00018410 
00018420 
00018430 
00018440 
00018450 
00018460 
00018470 
09010480 
00018490 
00018SOO 
00010510 
00018S20 
0001SS30 
00018540 
0001855  0 
00018560 
00018570 
00018580 
00018590 
00018600 
00018610 
00018620 
00018630 
00018640 
00018650 
0001866  0 
00018670 
00018680 
00018690 
00018700 
00018710 
09018720 
00018730 
00018740 
00018750 
90018760 
00018770 
PflDiaTBQ 
00018790 
OOOIPBOO 
00918810 
90018820 
00018830 
00018840 
00018850 
00018860 
00018870 
00018880 
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2P  CONTINUE 
30  CONTINUE 


DATE  s 78102 


10/«7/4a 


00018890 

00018900 


G LCVeL  21 


CIRCLE 


DATE  = 78102 


10/A7/48 


SUBROUTINE  CIRCLE tPOINTS«NPTS«XC«VC*lNCR> 

00018930 

00018940 

00818990 

1 T0LCK4«0CLT«T0L1«ANGST2«T0LPLT 

00018960 

NPCSO 

NNCsO 

00018970 

00018980 

DO  10  IslvNPTS 

00018990 

IP1»I*1 

IFll.EQ.NPTSI  IPlsl 

00019000 

00019010 

L=II-1)*2 

00019020 

LP1=«  1P1-1M2 

01SXl=P0INTStL*ll-XC 

000190 3C 
00019040 

0ISYl=P0INTSfL»2>-YC 

00019050 

OISX2=POINTS(LP1*1>-XC 

D1SY2*P.0INTS<LP1»2>-YC 

00019060 

00019070 

IF(ABSCOISXl>.LT.TOLl.ANO.ABS(OISYl).LT.TOLl> 

60 

TO 

20 

oooioooo 

IF(  A8S(OISX2)  .LT.  TOLl .A NO . ABS  (0 IS Y2) . LT .TQL 1 > 
if(abs(OISxi>.gc.toli)  go  to  3 

GO 

TO 

20 

00019090 

00019100 

IF(ABS(0ISX2>.GE.T0L1>  GO  TO  1 

00019110 

IFttOISVl/DlSY2).6T.O.O)  GO  TO  10 

GO  TO  20 

00019120 

00019130 

1 

IF(ABS(0ISV2>.GE.T0L1>  GO  TO  2 

00019140 

X3D1SX2 

GO  TO  a 

000191S0 

00019160 

2 

IF(ABS(0ISY1>0ISY2>.LT.T0L1I  go  to  10 

00019170 

a 

GO  TO  5 

IFIABSCOISYD.GE.TOLII  60  TO  32 

00019100 

00019190 

IF(ABS<0ISX2).LT.T0L1)  go  to  10 

00019200 

IF(ABS(0ISY2>.6C.T0L1)  GO  TO  10 
IF((DISX1/DISX2I.GT.O.O>  60  TO  10 

06019210 

00019220 

GO  TO  20 

00019230 

32 

IF(ABS(0ISX2>.GE.T0L1>  60  TO  34 

60  TO  36 

00019240 

00019290 

34 

IFIA8S(01SY2>.GC.T0L1>  GO  TO  36 

00019260 

XS0ISX2 

60  TO  8 

00019270 

38 

lF(ABSCOISY2-OISTl).LT.TOLli  GO  TO  10 

00019290 

IFCABS(0ISX2-0ISX ll.GE.TOLl)  60  TO  S 

00019300 

X«0ISX2 

00019320 

5 

60  TO  8 

SAsCDIST2-DlSYl)/IDISX2-DISXl 1 

00019330 

00019340 

SBs01SYl-SA*0ISXl 

00019350 

Xa>SB/SA 

1FIX.LT.OISX1.ANO.X.LT.OISX2)  GO  TO  10 

0G019S60 

00019370 

IF(X.6T.0ISX1.AN0.X.GT.0ISX2>  60  TO  10 

00019380 

8 

IF<ABS(XI .LT.TOLl)  GO  TO  20 

IFIX.LT.O.O)  NNC=NNC*1 

00019390 

00019400 

IF(X.6T.0.0I  NPC=NPC*1 

00019410 

10 

CONTINUE 

00019420 

IEf2*(NPC/2>.Ea.NPC)  RETURN 


0001 944  0 
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CIRCLE 


DATE  s T8102 


lO/47/«« 


1F(2*IHNC/2>.E0.NNC>  RETURN  00019450 

|NCR=1  00019460 

RETURN Q8QHA78 

20  INCRS2  00019480 

RETURN  00019490 
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AR  E« 


DATE  = 78102 


10/47/48 


SUBROUTINE  AREA<X«V*NPTS«YINT6) 

DINENSION  XdltTCl) 

XLOUsX(l) 

00019910 

00019920 

00019530 

XHI6HSXL0U 

00019540 

KHINsl 

00019550 

KNAXsl 

00  10  Is2«NPTS 

00019970 

IF(XI1).6E.XL0U>  60  TO  9 

00019980 

XLOUsXCI) 

00019600 

5 

lF(X(Ii«LE«XHI6H>  GO  TO  10 

00019610 

KNAXsI 

XHIGHsXd  > 

00019630 

10 

CONTINUE 

00019640 

IFIKniN.Ca.KNAX)  GO  TO  70 

00019690 

L1=KHIN 

00019660 

L2sKnAX 

00019670 

IFIKNAX.LT.XNIN)  L2  = NPTS*KMAX 

0001961)0 

NTA0S=L2>L1 

00019690 

riNTGsO. 

00019700 

K1=L1*L-1 

00019720 

IFIX1.GT.NPTS>  K1=K1-NPTS 

00019730 

K2=L1*L 

IF<K2.6T.NPTS>  K2=K2-NPTS 

00019750 

0XSX(K2I-X(K1) 

00019760 

20 

VtNT6sTINT6*(0X/2.)*(YtKl >♦Y<K2) ) 

00019770 

L2SKNAX 

00019780 

IFIKHAX.GT.KNIN)  L2=KNAX-NPTS 

00019790 

NTAaS=Ll-L2 

00019800 

DO  SO  L=1*NTABS 

00019810 

K1=L1-L*1 

00019820 

1FIK1.lt. 1)  R1=K1*NPTS 

K2=L1-L 

00019840 

IF(K2.LT.1>  K2=K2*NPTS 

00019890 

DX:X(K2>-X<Kl> 

0001 9860 

30 

YINTG=YINTG-I0X/2.>*IY(K1>*T(K2>) 

00019870 

YINTGsABSI YINTG) 

00019880 

RETURN 

70 

YL0M=YU> 

00019900 

YHIGH=YLOU 

00019910 

DO  80  I=2.NPTS 

00019920 

IFtTin.LT.YLOU)  YL0U=YII> 

00019930 

IF(Yd>.GT.YMIGH)  YHIGHsT(I) 

00019940 

80 

CONTINUE 

00019950 

YINTC=YHIGH-YLOW 

00019960 

RETURN 

00019970 

^Mi)  aaoi3i8JL 


213 


r 


1 


6 LEVEL 

21  PLOTOL  DATE  s 78102  10/47/48 

SUBROUTINE  PLOTOL(D«TAC»OATAS«MS«TITLES«6AnN**8PLOT*BPLOT» 

00019990 

1 VPROBLvlPLTOL) 

00020000 

DINENSION  IPLTOL(l) 

00020010 

DIMENSION  OATACdl.DATASID.PSd)  tTI TLES<1  > t APLOTCl > *BPLOTf  1 > 

00020020 

DIMENSION  XPR(  10) <OUT< 91)* SYMBOLIC) 

00020030 

DATA  BLANK*PERI00«FI/1H  «1H.*1HI/ 

00020040 

DATA  SVHB0L/1HA,1HB*1)«C*1H0*1HE«1HF»1H6*1HK«1HH/ 

00020050 

DIMENSION  HEA0(9) 

00020010 

URITEIA.IOI)  (TITLESCD.Isl.BO) 

00020070 

ICI 

FORNATdHltSOAl) 

00O200R0 

UR1TE(8«209) 

00020090 

MM  I TE  18.102)  GAMMA.0ATAC(4)  .OATAC(S) 

00020100 

ALOUsOATACdfi) 

00020110 

AHI6HsOATAC(17) 

00020120 

BHICN=CATAC(19) 

00020130 

BLOWsOATACdS) 

00020140 

0ELX=(BHIGH-BLOU)/90. 

00020150 

0ELY=(AHtGH-AL0M)/20. 

00020160 

KPMINTsS 

00020170 

NPANELsOATACt23)*O.S 

000201BO 

NPSI=0ATAC(9)*0.S 

00020190 

DPSIsOATACdl) 

00020200 

PSI=OATACdO)-OPSI 

00020210 

DO  20  I=1.NPS1 

00020220 

PSI=PSI*DPS1 

00020230 

20 

HEAD!  I )=PS1 

00020240 

TY2=AHI6H*0ELT/2. 

00020250 

DO  140  1=1*21 

000202  60 

TT1=TY2-0ELY 

00020270 

DO  30  JS1.91 

00020280 

30 

OUT<J)=BLANK 

00020290 

IFCKPMINT.EQ.S)  OOT<1)=FI 

00020300 

IF(KPRINT.NE.S)  OUTlDsPERlOO 

00020310 

IFCKPRINT.EQ.S)  0UT(91)=FI 

00020320 

IFIKPMINT.NE.S)  OUT (91 ) =PeR lOD 

00020330 

IFd.NE.l.AND.I.NE.21>  GO  TO  SO 

00020340 

1*9 

00020350 

DO  40  Jsl.91 

00020360 

.LSL*1 

00020370 

QUT(J)=PERIOO 

00020380 

IF(L.EQ.IO)  OUT(J)=FI 

00020390 

40 

IF(L.EQ.IO)  L=0 

00020400 

SO 

CONTINUE 

00020410 

HASsO 

00020420 

DO  90  K=1*NPANEL 

00020430 

DO  60  J*1*NPSI 

00020440 

OO  7S  Lsl*4 

00020453 

MAS*MAS*1 

00O20460 

IFIAPLOTCHAS)  .GT.TY2>  GO  TO  7S 

00020470 

IF(APLOT(MAS)  .LE.TYl)  GO  TO  7S 

00020480 

.XF(1PLT0L(NAS).E0.1)_.60  TO  75 

ODD  2 0490 

TX2s8L0W*0CLX/2. 

00020500 

214 


e LEVEL  21  PLOTOL  DATE  = 78102  10/A7/A8 

00  55  Hrl,91 

TX1*TX2-0ELX 

ifcbpldtcnas).lt.tx2.ano.bplot(masi.6f.txii  go  to  SB 

00820510 

00020520 

00020530 

55 

58 

75 

TX2=TX2*0EL* 

OUTinlsSYnBOLtJI 

CONTINUr 

00020540 

00020550 

00020560 

60 

90 

CONTIMUe 

CONTINUE 

TPLOT=TT2-OELT/2.  . 

00020570 

00020580 

00020590 

IF(KPRINT.EQ.5>  URITE(6t201)  YPLOT<OUT«YPLOT 

IF(KPRINT.NC.5>  URITE(6,205>  OUT 

IF<KPRINT.Ea.5>  KPR1NT=0 

00020600 

00020610 

00020620 

lAO 

KPRINTsKPRINT*! 

TY2=TY2-0ELY 

XPRCllsBLOU 

00020630 

00020640 

00020650 

150 

00  150  I=2tia 

XPRII )=XPR(1-11»10.»OELX 

URITEC6.2a4) 

00020660 

00020670 

00020680 

URITE(6t202>  XPR 

URITE(6«207) 

URITC(6.217}  VPROBL 

00020690 

00020700 

00020710 

217 

FORHATtlHO  tl05Xt 'PROB  ability: 'tFG.S  I 

WRITE(6>210) 

DO  160  L:1.4 

00020720 

00020730 

00020740 

K=9-NPSI*L 

IF(K.LE.A)  URITC(6*211)  SYPBOL (L> f HE  AO ( L) tSYMBOL (L*5) (HE AO (L*5) 
IF(K.GT.4.AN0.K.LE.9)  WRITE  16.21 2)  STHPOLIL ) .HEADIL) .SYPBOL{L«S} 

00020750 

00020760 

00020770 

160 

IF«K.GT.9)  WRITE(6«213)  SYH80L(L)  «SYPB0L(L«5) 

CONTINUE 

IFtNPSI.LT.Sl  UR1TEI6.215)  SYHP0LC5> 

00020780 

00020790 

0002D8Ca 

210 

IF(KPSI.GT.5>  WRITE(6«216>  SY NPOL ( 51 «HCAO( 5> 

WRITE (6,214) 

FORHATl/Z/IH  .22X.T6(**t)/1H  . 22X . •• • .74X. •* •/ IH  .22X.***.2X. 

00020810 

00020820 

0 002  08  30 

1 ‘KEY  TO  PLOT  SYM60LS»,53X  , •/IH  , 22X  , ♦ ,74X  ♦ • » ’/IH  ,22X, 

2 •••,4x, 'SYMBOL*. 5x**A/C  HEAOInG  (OEG > • »10X . * SYMflOL • . 5X, 

3 *kJC  HEADING  (DEG)  •.4X.***) 

00020840 

00020850 

00020860 

211 

212 

213 

FORMAT! IH  . 22X . •* '.TX *A 1 . 12X *F7.2. I8X *A 1 ,1 2X «F7.2 .9X« > 

FORMA T(1H  «22X.*.*t7X«Al*12X.F7.2.18X.A1.28X, >.*) 

FORMAT! IH  .22x.*.* .7X.A1.37X.A1.2ax.***> 

00020870 
00020880 
0002  0890 

2M 

215 

216 

FORMATdH  .22X»*«*»74X.*.*/1H  » 22X » • • • . 74X  . ••  •/ IH  .22X,76(  •••)  ) 
FORMAT!  IH  «22X.*'*«7X«  A1.66Xt  •••> 

FORMATdH  .22X.*»*.7X.A1.12X.F7.2.47X.».*) 

00020900 

00020910 

00020920 

1C2 

FORMATdHO.*SUN  ELEW AT  ION* , FI  2 .2/ IH  .•A/C  PITCH  ATTITUDE*. 

1 F7.2/1H  .'A/C  ROLL  ATTITUDE* .FS . 2//1H  .•OBSERVER  ANGLE*'93X. 

2 'OBSERVER  ANGLE'/IH  . 3X. 'DEGREES* .100 X. 'DEGREES*/ ) 

00020930 

00020940 

00020950 

20A 
20  2 
207 

FOR  MAT!/) 

FORMATdH  .9X.10F10.3) 

FORMAT!/. 38X. 'GLINT  AZIMUTH  WITH  RESPECT  TO  SUN  - DEGREES') 

00020960 

00020970 

00020980 

201 

205 

209 

FORMATdH  .F10.3.5X.91A1.1X.F10.3) 

FORMATdH  .15X.91A1) 

FORMATIIHO.SDX.'SUN  GLINT  SIGNATURE') 

00020990 

00021400 

00021010 

RETURN 

00021020 

215 


I 
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21  ORPLOT  DATE  = 78102  10/47/48 

SUBROUTINE  ORPLOT (OATACtORTAS »HS*  GANNA vALPHAZtBETAZ t 000210A0 

1 lNTRFL,ISHADU«INTERF,xPLOT*YPLOT*IBUF«JCONT*OUORKS*OWKSP2tVPROBL) 00021050 
DINCMSION  DATAC(l).OATAS<l>.NS(l)  «T  1 TLESd  ) .ALPHA2  < 1 > *BETA2  «1  > . 00021060 

1 *PL0T<1),YPL0TU)*IBUF(1)»INTRFL(  D.ISHAOUCD.INTERFd) 

OINENSION  0w0RKS(l>t0UKSP2(l) 

IZPTSsO 

00021070 

00021080 

00021090 

ICONTsl 

JC0NT=JC0NT*1 

IRFLTN=0ATAC(3>«0.5 

00021100 

00021110 

00021120 

theta =DAT*C«4  1 

PHI=DATAC«5) 

NPANEL=OATAC(23>*0.5 

00021130 

000211*0 

NPSI=0ATAC(9) «0.5 

IALPL=0ATAC(16I*0.5 

IF»OATAC(lb».LT.O  .)  I ALPL=0  AT  AC  ( 1 6 > -0  .5 

00021160 

00021170 

I0AL  = (DATAC  (1  7)-04TAC<  16)  )/2.  ♦O.S 
ieET=0ATAC(18>«0.5 

IF(OATAC( 18) .LT.Q. > I BETrOA TA C( 1 8 ) - 0 . 5 

00021190 

00021200 

00021210 

IDBET=C0ATAC( 19)-0ATAC( 18) )/6.*0.5 

8ETL=I8ET 

OBTrlOBET 

00021220 

00021230 

00021240 

ALPL=1ALPL 

DAL=IOAL 

NAS  = 0 

00021250 

00021260 

00021270 

0PSI=0ATAC( 11) 

00  90  I=1»MPANEL 

K=NS(  I) 

000212RO 

00021290 

ooojuao. 

NPTS=0ATAS(K) *0.5 

PSI=OATACC10)-OPSI 

00  60  J=1  .UPSI 

00021310 

00021320 

00021330 

PSI=PSI*DPSI 

NPTS=0 

00  75  L=1.NPTS 

00021340 

00021350 

00021360 

NAS=MAS*1 

IF(IRFLTN.NE.INTRFL(HAS))  go  to  75 

NPTS=NPTS»1 

0C021370 

00021380 

00021390 

75 

XPL0T(NPTS)=8ETA2CNAS) 

YPL0T(NPTS)=ALPHA2(NAS) 

CONTINUE 

00021400 

00021410 

00021420 

IFINPTS.EQ.OI  GO  TO  60 

CALL  SORT«XPLOT,YPLOT,NPTStOATAC« l7),DATAC ( 16)»0WKSP2) 

IF(NPTS.EQ.  0)  GO  TO  60 

00021430 

00021440 

_!Ii0114^ 

CALL  S0RT2<XPL0T,YFL0T.NPTStNPTSl tNPTS2»DU0RK S tO WKSP2 ) 
1F(HPTS2.EQ.0)  go  TO  77 

CALL  SORT<TPLOT.XPLOT,NPTSltDATAC( 19).DATAC( 18)«CUKSP2) 

00021460 

00021470 

00021480 

IFIHPTSI .EQ.O ) GO  TO  20 

CALL  CPLOT(ALPL»DAL.GANMAtTHETA,PHl,PSI«N»TSl«ICONT*XPLOT.YPLOT, 

1 IBUF.JCONT .aETL.DBT.VPROBL.IZPTS) 

00021490 

00021500 

00021510 

20 

ICONT=0 

CONTINUE 

DO  76  KKK:l.NPTS2 

00021520 

00021530 

00021540 

LNL=«KRK-1).2 

00021550 

217 
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OnPLOT 


DATE  = 78102 


10/A7/A8 


XPL0T<KKK>50y0RKS(LNL*l> 

76  VPL0T(KKK)s0y0RKS(LnL'»2> 

MPTSSWPTS2 

77  continue 

call  SORT(YPLOT«XPLOTtHPTS«OATAC( 19>f DATACI 18>,0yKSP21 

IFtHPTS.EQ.ni  GO  TO  60 

CALL  6PL0T(ALPL»0AL*GANHA*THETA*PHI,PSI,NPTS*IC0NT«XPL0T.YPL0T* 
1 IBUF«JCONT*BETL«DBTtVPROBL»IZPTS) 

TCONTsn 

60  CONTINUE 
90  CONTINUE 

IFdCONT.EQ.Q)  RETURN 

I2PTS=1 

CALL  GPLO T ( ALPL <0 AL*G ANNA t theta «PHItPSI*NPTStlCONT«XPLOT*YPLOT* 

I IBUPtUCQNJ.BCTLtOBT.VPROBLtlZPTS) 

RETURN 

END 


00021S60 

00021570 

00021580 

00021590 

00021600 

00021610 

00021620 

00021630 

Q002169B 

00021650 

00021660 

00021670 

000216RO 

00021690 

00021700 

00021710 

00021720 


£ LCVtL  21  SORT  OATC  = T8102  ia/«7/«8 


SUBROUT  INC  S0RT(XPL0T*YPL0T«NPTS«AH16H«AL0U«0U0RKS> 

00021730 

OinCNSlON  XPLOTd ) ,YPLOT(l) *0 WORKS! 1) 

000217AO 

INTCsO 

00021750 

LPTS=0 

00021760 

00  S I=1,MPTS 

00021770 

KPTslI-11.2 

00021780 

OUORKSIKPT*1>=XPLOT(I) 

00021700 

5 

0W0RKS<KPT*2>=TPL0T!I ) 

00021000 

DO  70  1=1.BPTS 

00021810 

IP1=I*1 

00021820 

IF!  I.CO.MPTS)  IPlsl 

00021830 

IFL1=2 

00021860 

IFL2=2 

00021850 

KPTs!  1-1)  .2 

00021860 

KPT1=!IP1-1>«2 

00021870 

X1=DU0RKS!KPT*1) 

00021880 

Y1=OWORKS!KPT*2) 

00021800 

X2=DU0RKS!KPT1»1) 

90021900 

Y2=DW0RKS!KPT1*2) 

00021910 

IFIYI.GT.AHICHI  1FL1=1 

00021920 

1F!Y2.6T.AHICM)  1FL2=1 

00021930 

IF!Y1.LT.AL0W»  IFL1=3 

000219A0 

IF!Y2.I.T.AL0U)  1FL2  = 3 

00021950 

IF1Y1.NE.Y2)  RATI0=IX2-X1 I/IY2-Y1 ) 

00021960 

GO  TO  !1,2*3>*  IFLl 

00021970 

1 

GO  TO  !lltl2«13>t  IFL2 

00021980 

? 

GO  TO  121.22.23).  IFL2 

00021990 

3 

GO  TO  131.32.33).  IF(.2 

00022000 

11 

IFIINTC.CQ.O)  GO  TO  60 

00022010 

LPTS=LPTS.l 

0D022020 

XPL0T!LPTS)=XSAVE 

00022030 

YPL0T!LPTS)=TSAVE 

000220A0 

INTC=0 

00022050 

GO  TO  60 

00022060 

12 

IFlINTC.ea.O)  60  TO  120 

00022070 

LPTSsLPTS.l 

00022080 

XPL0T!LPTS)=XSAVE 

00022090 

YPLOT!LPTS)sYSaVE 

00022100 

1NTC  = 0 

00022110 

12  0 

LPTS=LPTS*1 

00022120 

XPL0T!LPTS)sX1.RATI0»!AHI6H-Y1) 

00022130 

YPLQTILPTSISAHIGH 

nnQ9?i40 

Go  TO  60 

00022150 

13 

IFIINTC.CQ.O)  GO  TO  130 

00022160 

LPTSxLPTS.l 

00022170 

XPLOT!LPTS)sXSAVE 

00022180 

YPLOT!LPTS)=YSAVE 

00C22190 

130 

LPTS=LPTS*1 

00022200 

XPL0T!LPTS)SX1*RATI0*!AHI6H-T1> 

00022210 

YPLOT!LPTS)SAMIGH 

00022220 

XSAVC»X1*RATIO«!ALOW-Y1) 

00022230 

TSAVCsALOW 

000222AO 

G LCVtL  21 


SORT 


DATE 


T8102 


tO/*7/A8 


1NTC=1 

•0022250 

GO  TO  60 

00022260 

21 

■ 1 A 

LPTS=LPTS*1 

00022280 

XPLOTILPTSIsXSAVE 

00022290 

210 

LPTSsLPTS*! 

00022310 

XPLOT«LPTS)=Xl 

00022320 

XSAVE=X1*RATI0*<AMICH-Y1) 

00022340 

TSAVEsAHIGH 

00022350 

1NTC=1 

00022360 

GO  TO  60 

00022370 

22 

IFCINTC>EQ.O)  GO  TO  220 

00022380 

LPTS=LPTS*1 

00022390 

xplotclpts>=xsave 

00022400 

YPLOT(LPTS>=YSAVE 

00022410 

INTC=0 

30022420 

220 

LPTSaLPTS*! 

00022430 

XPL0T(LPTS>=X1 

00022440 

YPLOTftPTS)=.ri 

00022450 

GO  TO  60 

00022460 

23 

IF(INTC.CQ.O)  go  TO  230 

00022470 

LPTS=LPTS»1 

00022480 

XPLOT<LPTS)=XSAVE 

00022490 

YPL0T(LPTS)=YSAVE 

00022500 

230 

LPTS=LPTS»1 

00022510 

XPLOT<LPTS>aXl 

00022520 

YPLOT  (LPTS)=Y1 

00022530 

XSAVEsX1»RATIO*(ALOW-Y1) 

30Jfi2.540 

YSAVEsALOW 

00022550 

INTC«1 

00022560 

GO  TO  60 

00022570 

31 

IFdNTC.EQ.O)  GO  TO  310 

00022580 

LPTS=tPTS*l 

00022590 

XPLOT<LPTS)=XSAVE 

0002260Q 

YPC0T<LPTS)*YSAVE 

00022610 

310 

LPTS*LPTS»1 

00022620 

XPLOT(LPTS)sX 1«RATI0«<AL0W-T1 1 

00022630 

YPLOTll.PTS>=ALOG 

00022640 

XSAVESX1*RATI0*I  AHIGH-Y1> 

00022650 

YSAVCsAHlGH 

00022660 

INTCsl 

00022670 

60  TO  60 

00022680 

32 

IF(INTC.Ca.O)  GO  TO  320 

00322690 

LPTS=LPTS*1 

00022700 

XPLOT(LPTS) rXSAVr 

00022710 

TPLOTtLPTS)sYSAVE 

00022720 

INTC«0 

00022730 

320 

LPTS»LPTS»1 

00022740 

XPLOT<LPTS)sXl«ffATIO*IALOy*Vl > 

00022750 

YPLOTlLPTS) 3AL0W 

00022760 

220 
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SORT 


DRTf 
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60  TO  60 

33  IFCINTC.EO.O)  GO  TO  60 

LPTS=LPTS*1 

XPLOTILPTSIsXSAVE 

TPLOTtLPTS)=rSAVE 

ifac.=(i- 

60  IFtlNTC.EQ.O)  GO  TO  70 

IFd.LT.RPTSl  GO  TO  70 

l>PTS=l,gJS«l 

*PLOT(LPTS)=XSAVE 

TPL0T<LPTS)=TSAVE 

. 7 IL_  CONTINUE 

HPTS=LPTS 

RETURN 


00022770 

00022780 

00022800 

00022810 

00022820 

00022830 

00022840 

00022850 

00022860 

00022870 

j}o;22Bao 

00022890 
00022900 
JID  022210 


6 LEVEL  21 


SORTE 


DATE  s 78102 


10/«7/«a 


SUBROUTINE  SORT2(XPLOT*YPLOT*NRTSfHPTSl*NPTS2*OUORKS*OIIKSP2) 
OINENSION  XPLOT(1>»TPLOTC1I*0u0RKSC1I»0UNSP2<1) 
CORMOli/TDUIS/TOLEaM.TOLSNT.TOLOPT.TOLINT.TaLCKl.TOLCKS. 

Baa22«2B 

BB022930 

1 TOLCK4*OELT«TOLCIR<RNGST2*TOLPLT 

00022950 

RPTS2X0 

00022960 

l/.V 

BHINsXPLOT(l) 

80022900 

BNAXsBHIN 

00022990 

Mini  1 ■ 

1F(XPL0T<1).LT.BNIN1  BHlNxXPLOTd  > 

00023010 

IFtXPLOTII I.GT.BNAX ) BNAXxXPLOTd  > 

00023020 

Li. 

CONTINUE 

00023030 

IFIBNIN.GC.O.O)  RETURN 

00023040 

IFtBHAX.LE.O.O)  RETURN 

00023050 

DO  84  Ixl.MPTS 

80023080 

if<absixploT(I)).lc«45.>  return 

00023070 

8* 

CONTINUE 

00023080 

00023091) 

IFIBAVC.lt. ANGST2)  RETURN 

00023100 

00  30  Ixl.NPTS 

00023110 

LPSX(1-1).2 

00023120 

OMORKS(LPS*llxxPLOT(l> 

00023130 

CU0RKS<LPS*2)xrPL0TII> 

00023140 

Mill  11 

30 

CONTINUE 

00023160 

NPTSlxpPTS 

00023170 

CALL  SORT(yPLOTfXPLOT,HPTSl.lBO..O..OWKSP2) 

0002  3180 

00  35  I=1.HPTS1 

00023190 

LPS=«I-ll*2 

00023200 

CURSP7CLPS.ll:XPLQT(Il 

30023210 

35 

0UK$P2<LPS*2|sTPL0T< I> 

00023220 

00  40  ixl ,NPTS 

00023230 

LPS:(1-11  .2 

XPLOT(I)=OUORKS(LPS*ll 

00023250 

TPLOTI I >=OUORKS(LPS.2> 

0002  3280 

IFtXPLOTCI)  .GT.O.  > XPLOTd  IzXPLOTd  >-360. 

00023270 

«0 

CONTINUE 

00023280 

00  45  Ixl.HPTSl 

00023290 

Mil  l! 

0M0RKS<LPS*1)x0WKSP2(LPS*1> 

00023310 

45 

0U0RKS(LPS*2>=DWKSP2(LPS*2> 

00023320 

NPTS2XPPTS1 

Mill!  ■ 

NPTSIXPPTS 

00023340 

CALL  SORT  tTPLOT.XPLOT.PPTSl .0 . «-l 80 . tOUKSP2 1 

00023350 

RETURN  oaa^asfcfi 


END 


00023370 


6 LCVtI.  21 


SPLOT 


DATE  : 78102 


10/47/48 


SUB80U7INC  6PL07(ALPL«0AL«SUNtLV*7H«C«PH«C*HC*0*NPTS*lC0N7» 

1 XPL07«TPL0T«  IBUF«JC0NT»BCTL*0BTtVP80BL*IZPTS> 

BIHEHSIOW  XPL07CH.YPLOTtH.lBUf til 

CALL  FACT0K(.TB57> 

IF(ICONT.ca.O>60  TO  10 

BP  Tniaa«2a«2Bi8at2tti2ni8ai2C.i2fl»8tt«20«2flliU£QtH 

so  call  PL07<2.*1.»-S> 

CALL  PL0T(-2.*-l.*3) 

CALL  PL0Tta.B.-1..2l 

CALL  PL0T(8.8*13.t2) 

CALL  PL0T(-2.«13.«2) 

CALL  PL0T<-2.«-l.t2> 

BO  TO  30 

29  CALL  PLOT(0.f4.«.3l 

SO  TO  30 

80  CALL  PL0T<12.«-e.«-3> 

CALL  PL0T(-2.«-l.*3) 

CALL  fLOTia.Bi-l.iZl 

CALL  PL0T(8.8«13.*2> 

CALL  PL0T(-2.«13.,2) 

CAU^  PLaH^2.«>l«t21 

39  CONTIKOE 

C********  OPAW  AXIS  •••••* 

CALL  A21HUTM  t-18 18.  . 0 . .BETL.  DBT» 

CALL  AXIS(0.t0.*21HeLEVATION  OF  OBSCA VCR«2 1 *2. «90.* ALPL«OAL> 
CALL  PLOT(0.*1.«3> 

CALL  ELfll (A«tl«i2> 

call  PLOT(6««0.*3> 
call  PL0T<6.*2.«2I 

CALL  PLOT<0..2.,2I 

CALL  BRIOi.01*.01»l.«l.*B»2> 

CALL  STMBOL(5.*2.3«C.10«12HPROBA8ILITr  *0.tlll 

CALL  MmiBERI9B9..SBB..B.lB.VPR0BL.B..3l 

IFdZPTS.CQ.l)  60  TO  19 
10  CONTINUE 

XPLQTtNPTS«ll«XPLQTH> 

YPLOT<NPTS»U«TPLOTn) 

XPL0TINPTS«2)3BETL 

XPLnTINPTS*31»DBT 

TPLOT  < NP  T S*  2 » « ALP  L 
rPLOT<NPTS*3>sOAL 

HPxWPTS*! 

CALL  LtNC(XPLOT,rPLOT*NO»l»0«4> 

IFIHEAO.LT.O.)  CO  TO  IS 

AHAXxTPLOTIl) 

LNAXsl 

00  12  ls2«NPTS 

lF<TPLOT<H.CT.AWAXI  LWAX»I 

IFILNAX.CO.I)  APAXsVPLOTCLNAX) 

12  continue 

AIIAXBCYPLnTtLWAX»..ALPLl/PAL»a.aB 

BRAXsCXPLOT(LNAXI>BETL>/OBT 


88023380 
80823390 
JM«9Qa 
00023410 
0002  3420 
88823430 
8082  3440 
00023450 
80828460 
00023470 
000234RO 
00023490 

00023SOO 

00023510 

J0-023S28 

00023530 

00023540 


00023560 
0002  3570 
0 0 0 2 35  40 
00023590 
00023600 
00023610 
0002  3620 
00023630 
001123648 
00023650 
00023660 
00023670 
00023680 


00023690 
28823788 
00023710 
00023720 
0002  3730 
0002  374  0 
00023750 
8082  3760 


00023770 

00023780 


88623790 


00023800 

00023810 

28ft2.M?0 


00023830 


00023840 

80023850 

00023860 

00023870 

88823888 

00023890 


1 


I 


( 


1 
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DATE  = 78102 


10/A7/48 


CALL  NUMBER (8NAX«ANAX«0.07«HeAOt0.«-l) 

00023900 

GO  TO  19 

00023910 

15 

AHIN3TPL0T(1> 

Mllllllll 

LMINsl 

00023930 

00  17  I=2,NPTS 

IFlYPLOTdl.LT.AHIN)  LHINsI 

0002  3940 
000  2 395  0 

IFtLMIN.EO.Il  AHINsYPLOTILHINl 

00023960 

17 

continue 

APIN=(YPLOT<LHIN>-ALPL>/DAL-0.05 

00023970 

000239BO 

BMINsCXPLOTILMIN)  -BETD/OBT 

00023990 

19 

call  NUMBER(eMIN««RINt0.07tHEAD«0.t-l > 

CONTINUE 

00024000 

00024010 

IFIICONT.EO.O  GO  TO  *0 

0002*020 

CALL  SYMbOL<1.6*2.1«0.10<25HHEA0ING  OF  AIRCRAFT  (OEG>f0.t2S> 
CALL  SYHeOL(0.»3.0i.l0i22HSUN  ELEVATION  ANG LE  = .n..22> 

0002*030 

00024040 

CALL  NUMBER (999.t999.».10tSUNELVt0.tD) 

00024050 

CALL  SYMBOLIO.«2.8«. 10« 16HA/C  PITCH  ATT  = «0.«16> 

CALL  NUMBER (999.«999.» .10  «THAC.a.tl> 

00024060 

0002*070 

CALL  SYMB0L(0.*2.6..10«15HA/C  ROLL  ATT  = «0.«15> 

000240BO 

CALL  NUMBER (999. «999.,.10«PHACt0.*l> 

1FI(( JCONT-l>/5)*3.EO.<JCONT-l)l  CALL  SYMBOLi 3.4.3 .♦« 

00024090 

00024100 

1 0.14,19HSUN  glint  signature. 0., 19) 

00024110 

«0  RETURN  0002A120 

END  0002A110 


G LEVEL 

21  GRID  0«TE  = 78102 

10/A7/A8 

SUBROUTINE  OR  10 IX ORG * YORGtDELVLS t DELHLS«NVSPSt NHSPS > 

00024140 

OYN«XsNHSPS*OELHLS 

00024150 

DX«AX=NVSPS»OELVLS 

00024160 

XnAX=XORG*OXnAX 

00024170 

yhaxsyorg^otbax 

00024180 

NWLSsNVSPS*! 

00024190 

NHLS=NHSPS*1 

00024200 

XPRESsXORG 

00024210 

YPRES=YORG 

000g*2?0 

DO  10  I=1«NVLS 

00024230 

CALL  PLOT (XPRES*YPRESt3) 

00024240 

FACT=<-1) **(I*1) 

00024250 

YPRES=YPRES*F  ACT*OTNAX 

00024260 

CALL  PLOT<XPRES.YPRES,2) 

00024270 

xpres=xpres*oelvls 

00024280 

CALL  PL0T<XPR£StYPRES*3) 

00024290 

10 

CONTINUE 

00024300 

XPRES=XORG 

0002431J 

YPRES=YORG 

00024320 

DO  20  I=ltNHLS 

00024330 

CALL  PLOTIXPRES.YPRES.3) 

afi524.318. 

FACT=<-1) •*«I*1) 

00024350 

XPRES=XPRES*FACT*DXMAX 

00024360 

CALL  PLOTIXPRES.YPRES.2) 

00024370 

20 

YPRES=YPRES*CELHLS 

00024380 

RETURN 

00024390 

END 

03024400 

